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SYNOPSIS 

Reverse micelles are the aggregates of amphipathic 
molecules such as synthetic surfactants and phospholipids in 
nonpolar media. These aggregates are capable of solubilizing 
large amounts of water in their polar interior, which is then 
called a 'water pool'. The surfactants are located only at the 
interface, with their polar head-groups directed towards the water 
pool and their tails sticking out into the nonpolar continuous 
phase. Such structures have received much attention in recent 
years as they constitute suitable working models for exploring 
physico-chemical properties of cell membranes, and also because of 
the potential use of the water pool as a microreactor for many 
enzymatic and nonenzymatic reactions. 

The surfactant chosen for the work embodied in this 

thesis is lecithin (phosphatidylcholine or PC) - one of the most 

commonly found glycerophospholipids which constitutes the bulk of 

the phospholipid content in animal cell and bacterial membranes. 

Reverse micelles of PC have been formed with cyclohexane and 

carbon tetrachloride (CCl^) as two typical nonpolar media. In the 

first part of this work, attention has been mainly focused on a 

variety of NMR studies of the pool water and the polar head-groups 
+ 

(-N(CH2)3 and phosphate) of lecithin. In the second part, the 
design and construction of a versatile magnetic field gradient 
controller (MFGC) is reported for making NMR measurements. 

In Chapter I, the present state-of-the-art pertaining to 
research in phospholipid reverse micelles is presented. 


In Chapter II, we report a new lecithin-based reverse 
micelle with dicyclohexyl carbodiimide (DCC) , a popular 'peptide 
condensing agent', as co-surfactant. Systematic water uptake of 
this reverse micelle at the chosen molar ratio of 5 ; 1 of PC to 
DCC in cyclohexane shows that a maximum of 15.8 moles of water 

can be solubilized per mole of PC. Based on and ^^P NMR 

chemical shift studies, a suitable model is proposed for the 
interaction of water with the polar heads of lecithin and DCC. 
Water proton T^'s recorded as a function of water addition are 
analyzed theoretically with the help of the well-known 

Bloembergen- Purcell- Pound model, and a rationalization is made 
of the rotational correlation time (x ) and microviscosity (tj) of 
water calculated from our experiments. The water proton 

measured from the experimental spectral linewidth is analyzed 
theoretically assuming that only dipole-dipole relaxation is 

dominant. A suitable explanation is also given for the 
observation of a contrasting trend in the P T^values observed 
for this reverse micelle. 

Chapter III describes our studies of intramolecular 
water dynamics and exchange of water between two different types 
of sites ('bound' and 'free') through deuteritim magnetic resonance 
(DMR) spectroscopy. For these studies we have used the 
PC/DCC/cyclohexane with D 2 O in place of H 2 O, as well as its 
non-DCC counterpart for comparison, and have performed our 
experiments at two different temperatures (263 K and 296 K) . Also 
presented in this chapter are (i) a general DMR theory, (ii) a 
very brief quantum mechanical background for calculating nuclear 
spin systems' response to radio-frequency (rf) pulse (s) and to the 
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quadrupolar echo pulse sequence used for recording DMR spectra, 
and (iii) a computational method used for simulating experimental 
spectra. The two-line DMR spectrum observed in both the reverse 
micellar systems at all D 2 O concentrations is correctly assigned. 
Our computed lineshape fittings show that the exchange rates of 
water between the 'bound' and 'free' sites are nearly 10^ times 
less at 263 K than at 296 K, the T 2 values of 'bound' water being 
10 times less than those of 'free' water. Interesting results 
obtained from both the reverse micellar systems at 263 K are also 
discussed. A tentative conclusion is drawn at the end of this 
chapter about water 'freezing over' at one of these sites at 
263 K, in comparison with analogous hydration studies of lipid 
head-groups . 

Chapter IV briefly reviews a recent technique for using 
the one-dimensional analog of a 2-D proton multiple quantum (MQ) 
filtered correlation spectroscopy for observing 'forbidden' cross 
peaks and extracting therefrom some useful cross-relaxation 
effects among coupled nuclear spins in a molecular aggregate. We 
have shown that, by fitting the MQ-filtered profiles of forbidden 
peak- intensities versus evolution time to our theoretical 
equations, valuable structural and dynamical parameters can be 
obtained for a specific molecular group in a macromolecule. 
However, calculations of the theoretical response for a group as 
large as the - N(CH 2)2 (which is a coupled case) turn out to be 
too complicated, as shown in this chapter. This has constrained 
us to use the simpler, well-studied sodium deoxycholate (DOC) as a 
surfactant for forming reverse micelles. The results of our 
zero- and double-quantum filtered 1-D experiments of this reverse 
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micelle have been interpreted using a suitable motional model for 
the group of the DOC molecule, which provides a good 
illustrative support for the theory developed in this chapter. 

Chapter V comprises work on instrumentation development. 

In the first part of this chapter, we review the different types 
of the magnetic field gradients (MFC) and the methods for 
generating them. The design of a versatile magnetic field 
gradient controller, adapted from the work of Fitzsimmons 
[Rev.Sci. Instrum. 1338 (1982)], with detailed circuit 
diagrams and performance results for two types of applied MFG, are 
given. The improved features of our magnetic field gradient 
controller over those of Fitzsimmons' are summarized. In the 
second part of this chapter, the results of our experimental 
measurements of the self-dif fusion coefficient (D) of 'pool' water 
in two reverse micellar systems ( PC/cycloh e xane/H 2 O and 
PC/CC1^/H20 ) using a commercial NMR equipment are presented. The 
D~value trends for the 'pool' water trapped in these reverse 
micelles substantiate some earlier findings of our laboratory. 

In Chapter VI we summarize the results of this thesis, 
and suggest some further work which should be done on 
phospholipid-based reverse micelles for an enhanced understanding 
of these fascinating molecular assemblies. 

The Appendix contains a listing of our FORTRAN 77 
program which simulates the experimental DMR lineshape discussed 
in Chapter III. 
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CHAPTER 1 

GENERAL INTRODUCTION 

I.l Reverse Micelles 

Reverse micelles are the aggregates of amphipathic 
molecules such as synthetic surfactants and phospholipids in 

. M— AN 

nonpolar media. ' ' These aggregates are capable of 

solubilizing a large amount of water and several other polar 
solvents, e.g., glycerol and formamide, forming in the case of 
water, a so-called 'water pool^.^^^ While reverse micelles that 
possess water cores are water-in-oil (w/o) microemulsions, they are 
still called reverse micelles in analogy with micelles. In 
reverse micelles, the polar 'head' groups of surfactant molecules 
are concentrated in the interior of the aggregate and directed 
towards the water pool while their hydrophobic 'tail' groups 
protrude into, and are surrounded by, the bulk nonpolar solvent. 
The surfactant dissolved in organic phase is found to reduce the 
interfacial tension between the water and the organic phase, which 
is the prerequisite for the formation of reverse micelles. These 
aggregates are thermodynamically stable, homogeneous and optically 
transparent . 

Reverse micelles have been the subject of recent 
increasing interest not only because of the water core's catalytic 
effects and their accommodative nature to host proteins and 

small molecules^ , but also because these multicomponent 

. . . (7) 

structures provide a model for processes in biomembranes' ' 

including transbilayer transport and fusion. While there 



is a considerable accumulation of experimental results about these 

aggregated systems, the information is incomplete. An unambiguous 

quantitative interpretation of micellar behaviour is still 

missing. In particular, much work is needed on the dynamic 

behaviour of reverse micelles, investigating the kinetics of their 

formation, their catalytic effects, and the solubilization of 

other molecules by the aggregates. The experimental and 

theoretical results accumulated on reverse micelles till 

f 3 1""3 8 ) 

recently are summarized in published monographs' ^ and review 

articles. The references given here are merely meant to be 

representative, and give by no means a complete view of the many 
studies which have been published. 

In this chapter, we restrict ourselves to reviewing the 
work pertaining to phospholipid-based reverse micelles only. In 
the following sections, we shall give a brief summary of the 
general structure of phospholipids and their different aggregated 
forms, the possible biological relevance of reverse micelles, 
their thermodynamics and phase diagrams, and the wide variety of 
methods (physical and chemical) used to characterize/probe reverse 
micellar structures. Finally a table containing important 
phospholipid— based reverse micelles studied till recently is 
presented . 


1.2 Phospholipids ; Source and Structure 

The biological membranes play a key role in cell 
biology. They control the cell life by modulating transport and 
exchange processes between endocellular (inside) and exocellular 
(outside) compartments and are involved in biosynthesis, energy 
transduction, information 


transmission. 


and 


cell 
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•a.- (41-44) 

recognition. The basic structure of a cell membrane is a 

bilayer lamellar sheet, approximately 4-5 nm in thickness, 
comprising polar and neutral lipids and proteins that are embedded 
in the bilayer (integral) and bound at the surface (peripheral) . 
The Singer and Nicholson fluid mosaic model of a membrane 
structure is given in Fig. I.i. 

Phospholipids are the primary structural components of 
all biological membranes. These are polar lipids containing 
phosphoric acid as mono- or di-ester. The natural 

phospholipids are classified into two main groups , namely 

(i) Sphingophospho lipids : These are derived from sphingosine. 

Its general formula is given in Fig. 1.2(a), where R 

represents a hydrophilic group and R' an acyl group. The 

most common sphingophospholipid is sphingomyelin, where the R 

+ 

group is -CH 2 CH 2 N(CH 2 ) 3 . 

(ii) Glycerophospholipids ; These have glycerol as a backbone with 
one or two 0-acyl or 0-alkyl, or O- ( 1-alkenyl ) chains 
attached at positions one and two and phosphate at position 
three. Its general chemical structure is given in Fig. 
1.2(b). The R' and R" groups are mostly derived from fatty 
acids, generally unbranched and with an even number of carbon 
atoms. Fig. 1.2(c) shows the main R groups found in natural 
membranes and their corresponding common names. 

Although they are different in terms of structure and 
function, both are amphipaths. Further details are given 
below about glycerophospholipids only , since these are the most 
commonly found membrane lipids and also since these appear 
exclusively in phospholipid-based reverse micelles. 
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Fig. 1.1 


The Singer and Nicholson 'fluid mosaic' model of a 
biomembrane, (adapted from ref. 45) . 
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Fig. 1.2 


(a) 


(c) 
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General Chemical structure of sphingophospho- 
lipids (a) glycerophospholipids (b) & main clases of 
glycerophospholipids with their common names (c) . 
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Lipid Polymorphism 

Phospholipids are seldom molecularly dispersed in 
solutions. As amphipathic molecule, phospholipids have a great 
tendency to aggregate and, like fatty acids, to form more than one 
structurally different phase in the crystalline (solid) state as 
well as in aqueous dispersions. This property is referred to as 
polymorphism. Various structures have been mainly established 

by X-ray diffraction , electron microscopy and magnetic 

(33,52) 
resonance. ' ' ' 

The most common phases of phospholipids in 
f40 49 53^ 

water' ' ’ ' are given in Table I.l. In the Table, L indicates 

a one-dimensional lamellar lattice, H, P and Q represent a 
hexagonal two dimensional, rectangular two-dimensional and cubic 
lattices, respectively. The Greek letters and a/3 indicate 

the conformation of the hydrocarbon chains. 

Reverse micelles, the subject of this thesis, belong to 
^ (reversed two-dimensional hexagonal type), P^^_^ (reversed 
two-dimensional tetragonal type) and a (reversed 

three-dimensional cubic type) containing a water core in an 
organic environment. 

1.4 Existence of Reverse Micellar Structures in Biological 

Membranes 

Although the existence of non-bilayer structures in 

biological membranes is by no means established conclusively, the 

existence of reverse micelles of membrane proteins and 

. r54) 

phospholipids in organic solvents is established. ' ^ The 

. . f55-57) 

possibility of reverse micellar dynamic intermediates' ' 
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their roles have been 


analyzed. 


These and the studies on 


Table I.l The most common phases of phospholipids in water 


Type of Packing 


Structure 


t 


Coherent double-layers 
of molecules separated 
by water; lamellar one- 
dimensional array 





mmmm 

mtmmm 

'-fi' 


Rod-like aggregates with 
organic core and 

aqueous core ; 

hexagonal two-dimensional 
array 


Rod-like aggregates with 
organic core and 

aqueous core '• 

tetragonal two-dimensional 
array 





B ••• *^55 







iaSK tPilwr affy * 

Pl-a 



Pll-a: 


Rod-like aggregates with 
organic core in aqueous 
environment; rectangular 
two-dimensional array 


Rippled double layers 
of molecules separated 
by water; lamellar two- 
dimensional array 







Pin-o: 




Uiiii 


imn; 

JIUl 


ttnnuiii^iti 

ttiiiuiiimliwr 




Pc/J 


Rod-like aggregates with 
aqueous core (Q^) and 

organic core » cubic 

three-dimensional array 



w~m 


^adapted from ref. 40,49a, 53 and 68. 





lipid polymorphism led to a new fundamental insight into the 

physical properties and functional roles of lipids in membranes. 

A new understanding of biomembranes necessitated a revision of the 

well-known 'fluid mosaic' model. The new model proposed by Cullis 
r 62 63) 

et ^1* ' , called the 'metamorphic mosaic' model, reproduced 

in Fig. 1.3, includes possible processes which involve reverse 
micellar structures such as exo- and endocytosis membrane fusion, 
movement of lipids across the bilayer , transport of ions across 
the bilayer and perhaps many other processes yet undetected. 

I . 5 Thermodynamics and Phase Diagrams 

Amphipathic or amphiphilic molecules aggregate in 
nonpolar solvents predominantly due to dipole-dipole and ion pair 
interactions and also due to a decrease in the overall free 

energy of the system. The hydrophobic effect of these molecules 
provides the driving force for aggregation and the repulsion 
between head groups limits the size of the aggregated structure. 
The nature of solvent (its polarity, dielectric constant etc.) 
also plays a role in the aggregation behaviour of surfactants. 

The thermodynamic principles of ampliphatic self- 
aggregation and the geometric aspects of the phenomenon have been 
well discussed. Israelachvili and Coworkers 
discussed these in terms of a parameter called 'critical packing 
parameter' which is connected with the thermodynamic 
characteristics of the amphiphile. The enthalpy and entropy 
changes of aggregation of potassixam benzene sulfonate in 
heptane^®®^ (enthalpy (AH) = -79.5 IcJ mol”^; entropy(AS) = -62.8 
JK”^ mol”^) suggest that the association energy for reverse 
micelles to be mainly of enthalpic rather than entropic origin as 
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Fig. 1.3 
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Integral 

protein 


Th© 

^ mosaic model of biological 

(3) Show •” (1) and 

reverse involvement of reverse micellar or 

exocvtof- 1* structural intermediate in an 

®oleculai 4- process and transbilayer 

2 shows process, respectively. Region 

structure *4- of reverse cylinder 
re at interbilayer tight junctions. 
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The thermodynamics of amphiphile aggregation in polar 
and nonpolar solvents is usually discussed in defining a term 

f 4\ 

called 'critical micellar concentration (CMC) ' ' - a concentration 

of amphiphile at which aggregation occurs (it is to be noted that 

(71-73 ) 

the definition of CMC is not unique). Three basic models' 

have been used for the thermodynamic treatment of micellization. 

They are : pseudophase, single step equilibrium or mass action law 

and multistep equilibrium. The first model is frequently used in 

NMR studies of micelles. Giomini et al.^^^^ used the mass 

action law model for finding out the CMC and other thermodynamic 

parameters of lecithin/benzene/water reverse micelle. Although 

these theories can in principle be applied to phospholipid 

aggregates, much progress has not been made due to the following 

two reasons ; (i) The cooperativity of reverse micellar formation 

is much smaller compared with that of normal micelles. When 

the aggregation number is small, the change in physical properties 

with increasing concentration of the surfactant is gradual and, 

therefore, a well-defined CMC can seldom be identified. 

(ii) Since minute amounts of water are invariably present in the 

organic solvent, these system should be considered as ternary and 

the CMC would then lose its significance, as it changes with the 

(3) 

amount of water. ' ' 

Despite the lack of a clear-cut concept of CMC for 
reverse micelles as mentioned above, the "CMC” has been 
determined for many surfactants in nonpolar 

solvents ^(4/3,74,77,78) Another way of understanding the 

thermodynamic properties of these systems would be through the 
construction of ternary phase diagrams. To the best of our 
knowledge, only two reports have been ptiblished on lecithin based 
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reverse micelles^^^' Shervani et have studied 

lecithin-based reverse micelles in benzene, cyclohexane, n-octane 
and n-dodecane, in the presence and absence of cholesterol. 

The role of water in the formation of reverse micelles 
has also received considerable attention. Previous studies 
have concluded in general, that water is a prerequisite to promote 
the aggregation of amphiphilic surfactants to reverse micelles and 
that part of the water present serves as 'gluing' agent between 
the polar head groups of the surfactant monomers in the core of 
the reverse micellar aggregates. However, Yu et al. reported 
recently from their dynamic and static light scattering studies on 
sodium bis (2-ethylhexyl) phosphate in n-hexane that water is not a 
prerequisite for reversed micellization, in contrast with earlier 
reports. Instead, water can function as an antimicellization 
agent 


1.6 Reverse Micellar 'Water Pool' 

A lot of importance attaches to the physico-chemical 
studies of water in reverse micelles is because of its rather 
peculiar chemical and physical properties^ ' ' ' which differ 

from those of normal water, especially at low water-to 
-surfactant ratios. For instance, it has been found by various 
spectroscopic and other physical methods that these water 
molecules , exhibit much higher viscosities ' ' , low 

mobilities , varied polarity and dielectric 
constant , varied and hydrogen bonding and low 

freezing point. 

Considerable studies have also been done on the 
interaction of water with the polar head ‘ groups of 
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phospholipids . 110) However, the conclusions of some of 

these studies are often contradictory or 
tentative . 

Reverse micelles are not rigid but dynamic entities. They 

can collide freely in the hydrocarbon solution which leads to the 

exchange of water pools. The dynamics of the exchange of water 

in reverse micelles and biological systems can be described by 

qualitative and quantitative models. Qualitative models assume 

that water molecules exist in two or more environments, with one 

of these environments being that of pure water. The simplest case 

would be a two-state model classifying the water molecules as 

'bound' and 'bulk' . This was later on extended to a 

three-state model 116)^ three types of state being : 

Type I ; "bulk Water" or "Free Water" - the rotational and 

translational motions of these water molecules are not appreciably 

altered compared to normal water, with t , the rotational 

c 

. . . —12 
correlation time, being ~ 10 s. 

Type II : "bound water" - the rotational motions, freezing point, 

etc. of these water molecules are significantly hindered by 

interactions with the polar head group of surfactants, leading 

-9 

to values of ~ 10 s. 

c 

Type III : "trapped water" or "irrotationally bound water" (this 

term, however, does not exclude translational or exchange 

mechanisms) ; these water molecules are essentially "trapped" or 

"site bound" to the surfactant or to the macromolecules, with - 

c 

lO”^ to lo"^ s. 

Quantitative estimation of these types of water has also 
been made. Koening and Schillinger ^ estimated Type III 

water from the magnitude of NMR spin-lattice relaxation time (T^^) 
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dispersion. Type II water can be directly estimated from freezing 
experiments, and the remaining water is to be considered as Type I 
or bulk water. Although previous studies have led to a reasonably 
consistent picture of the mobility and organization of water in 
the 'water pool', there still exists the necessity for obtaining 
true (or at least reliable) values for the physical properties of 
the aqueous core of the reverse micelles by using a variety of 
techniques. 

1.7 Experimental Methods Used to Probe Reverse Micelles 

Various experimental methods have been used to 
characterize the reverse micelles. Most of these methods are used 
to probe the compartmentalized, polar microenvironment, the 'water 
pool ' , but some of these are used to get a macroscopic information 
such as bulk viscosity . 

We may broadly categorize the experimental methods into 
three groups : 

( i ) scattering techniques 

(ii) physical methods 

(iii) spectroscopic techniques 
(i) Scattering Techniques 

Scattering techniques give a macroscopic structural 
information about reverse micelles. The three most frequently used 
scattering techniques for studying reverse micelles are ; light 
scattering (LS) , small-angle neutron scattering (SANS) and 
small-angle X-ray scattering (SAXS) . 

Light scattering techniques essentially probe the 
Brownian motion of particles. Static and dynamic light 

scattering (also known as quasi-elastic light scattering or 
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photon correlation spectroscopy) have been used to determine size 
and translational diffusion coefficients of droplets, to measure 
intermicellar interactions of reverse micelles and to investigate 
protein encapsulation in reverse micelles. A good review of these 
studies is given in Ref. (37) . 

Light scattering studies of egg lecithin reverse 
micelles in benzene, carbon tetrachloride (CCl^) and cyclohexane 
were reported from our laboratory . These studies indicate the 
presence of homodisperse aggregates in benzene and CCl^ and 
anisotropic tubular structures in cylohexane. 

The small angle neutron scattering (SANS) technique 

gives more detailed information on size and size distributions of 

reverse micelles, due to the fact that the neutron scattering 

vector has a large range (0.02-0.2 A* ^) . This can also be used 

to study the effect of additives in the reverse micelles.' 

Translational motion of reverse micelles can be monitored by using 

a quasi— elastic, incoherent neutron scattering method^ 

Several comprehensive reports of SANS studies in micellar and 

(37,146,151,152) 

microemulsion systems have been published. 

A SAXS experiment determines the so-called 'radius of 
gyration', which is related to the polar core radius of reverse 
micelles. These can be used to probe the change (increase or 
decrease) in size of the polar core in the presence of additives 
or solubilized enzymes 
(ii) Physical Methods 

There are quite few physical methods that have been used 
to probe reverse micelles. They are : 

(a) Viscosity measurement 

(b) electrical conductivity 
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(c) electron microscopy 

These have been used on egg-lecithin reverse micelles in 
three different solvents in our laboratory earlier The major 
results of these studies are : the reverse micelles contain 
dispersed water in small spherical droplets surrounded by a 
phospholipid layer; the electrolyte (KCl solution) added for 
measuring specific conductivity is contained in discrete water 
pools that are formed in the interior of the reverse micelles. The 
cyclohexane system at water-to-surfactant mole ratio of 5-6 
exhibits an elongated tubular structure containing aqueous canals, 
(iii) Spectroscopic techniques 

A wide variety of spectroscopic techniques may be used 
to probe the water pool, organic bulk phase and surfactant to get 
dynamic or motional information of the system. 

The selection of a particular technique depends on the type 
of physical property in which one is interested. A sximmary of 
characteristic frequencies of molecular motion of lipids and the 
approximate frequency ranges over which various spectroscopic 
techniques are sensitive to molecular motion is given in Fig. 1.4. 
This figure illustrates that the motional frequencies of different 
parts of the lipids occur in the range of lo”^ - 10^^^ s 
i.e., spanning a time-scale of 20 orders of magnitude. We note 
that some techniques yield a static picture since the molecular 
motions are 'slow' relative to the measurement and some see a time 
averaged picture of molecules which move between different 
environments very rapidly relative to the measurement time. 
Fig. I. 4 also indicates that out of all these spectroscopic 
techniques, NMR can probe the widest range of molecular motions 
thereby becoming an important tool for motional studies. 
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Fig. 1.4 Approximate frequency ranges of various 
spectroscopic techniques and molecular motions of 

The time scale of a particular motion being studied and the 
time scale of the technique being used for the measurement worth a 
further emphasis. Consider a general case which could equally 
apply to an EPR, fluorescence or ^H-NMR measurement of molecule or 
probe. All three spectroscopic techniques are sensitive to the 
orientation of the molecules with respect to an axis system 
defined in the laboratory. For example, in ^H-NMR, the observed 
spectrum is sensitive to the orientation of C-D, 0-D or other bond 
in relation to the externally applied field. In EPR, the 
situation is similar, where the spectrum is sensitive to the 
orientation of the N-0 (nitroxide) bond present in most commonly 
used probes, with respect to the applied magnetic field. In 
fluorescence spectroscopy, the measured polarization of the 
emitted light is dependent on the orientation of the molecular 
transition dipole moment with respect to the direction defined by 
the polarizer used to make the measurements. Some further 
elaboration of these techniques is given below : 
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(a) Optical Spectroscopy (UV-visible, fluorescence, IR, Raman and 

circular dichroism) . ' 148) 

In most of these experiments, some probe molecules are 

entrapped in either the water pool or the organic phase of 

micelles, and conclusions are derived from the observed trend of 

properties such as polarity, hydrogen bonding, mobility, 

vibrational motion of a bond in a molecule, protein confirmation, 

activities of solubilized enzymes, etc. 

For example using the nitrate ion as a probe and 
* 

monitoring the n tt band in the UV range, in egg lecithin based 

reverse micelles in benzene, carbon tetrachloride and 

cyclohexane' ^it has been found that the effective polarity of 

the water pool is lower than that of bulk water. This result is 

substantiated in the fluorescence emission spectra of the same 

three systems by using 8-anilinonaphthalene sulfonic acid 

(ANSA) as a probe. In the cyclohexane system, this technique 

indicated a possible presence of liquid crystalline phase when the 

water to surfactant ratio was six. The near IR spectra of the 
( 97 ) 

same systems' * give two absorption bands : one at about 1890 nm 
corresponding to water-dispersed in the organic phase and the 
other band over a range of 1930-1960 nm corresponds to the water 
solubilized in the reverse micellar interior. From this study , it 
is also estimated that only about 0.05% (V/V) of total water is 
present in the bulk solvent. 

The following spectroscopic techniques have also been used, 
but to a lesser extent, to characterize the water-pools of reverse 
micelles, namely. 
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(a) time-domain dielectric spectroscopy , 

M 54) 

(b) positron annihilation spectroscopy' , 

. QSS) 

(c) photon correlation spectroscopy' ' , 

(d) perturbed angular correlation spectroscopy . 
However, these shall not be further described here. 

(b) Magnetic Resonance Techniques 

Electron Paramagnetic Resonace (EPR); 

Electron paramagnetic resonance (EPR) technique is very 
sensitive and extremely useful in the study of reverse micelles. 
It is usually carried out by incorporating either a paramagnetic 
ion^^^^^ or spin labels or probe in the water pool or organic 
phase. The properties that can be obtained from EPR are : the 
aggregation behaviour of different surfactants, microviscosities 
and local polarities around the probe molecule, the effect of 
different solubilizates on these properties, and the dynamics of 
micellization processes. 

A wide variety of spin labels have been designed and 
synthesized to have physico-chemical properties such that they 
either remain in the water pool (e.g., Fremy's salt 

(Na 2 (SO”) 2 NO) , CUCI 2 ) or organic phase (e.g., 16-doxyl stearic 
acid) , or they are exclusively present at the water/oil interface 
(e.g., long chain fatty acid spin label analogs; substituted 
semiquinones , 5-doxyl stearic acid etc.) or they exchange between 
different compartments (e.g., 5-doxyl hexylsulfate, sodium salt). 
Nuclear Magnetic Resonance (NMR) : 

Since all nuclei with nonzero nuclear spin quantum 
number are NMR 'active', useful information can be obtained from 
virtually all components of the reverse micelle. This essentially 
eliminates the need for using a molecule which may be 'invasive' 
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and change the physico-chemical properties of the system. 

The NMR parameters are of two kinds. The first kind 
consists of steady state parameters (the chemical shift (5) and 
scalar (J) , dipolar- and quadrupolar coupling constants) which are 
related to the molecular structure of the system under study; they 
give information about the chemical environment of the observed 
nucleus and about its stereochemical relationship with other 
nuclei in close spatial proximity. The second kind includes 
dynamic parameters (the longitudinal or spin-lattice relaxation 
time (T^) and the transverse or spin-spin relaxation time (T 2 )) 
which are related to the motional characteristics of the nuclear 
spin system. 

The following information can be obtained from NMR 
studies on reverse micelles : aggregation numbers and 
exchange rates of different components of the system, extent of 
hydration of the polar head group, properties and activities of 
the solute solubilised in different phases, macroscopic 
self-diffusion of a molecule and segmental motions, conformations 
and order parameters of surfactant molecules, etc. Much of the 
earlier NMR literature on reverse micelles has been somewhat 
contradictory (see Table 1.2) . We shall review very briefly the 
kind of information that has been obtained from NMR studies of 
egg-lecithin reverse micelles after 1980. 

Systematic studies have been carried out in our 

(92 94 97 110) 

laboratory on egg-lecithin based reverse micelles' / > > ' m 

131 13 

benzene, CCl^ and cyclohexane using extensively, H, P and C 

2 

NMR and to a limited extent H and water proton diffusion 
coefficient measurements. Based on these experiments several 
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guest molecules, peptides, enzymes etc. which are solubilised 
in the ‘water pool’ and bulk organic solvent (an updated and 
and integrated version of Table III, IV and VI of Ref. 40) 


Sol- R = 
vent 

[Water] Guest molecules, 

tEgg enzymes etc. 

Lecithin] 

Method (s) used Ref. No. (s) 

Benzene 

0-45 

Light scattering. 

1,18,92, 



Osmometry, viscosity. 

94,97, 



diffusion, electron 

101,103, 



microscopy, NMR, IR, 

105,108, 



Fluorescence , elec- 

110,119, 



trical conductivity 

120,138, 



and optical 

142,144 



biref ringence 



0-16.2 

KCL 

Electrical 

conductivity 

92 

0-30 

NaCl,KCl, MgCl^.CaCl^ 

NMR and IR 

127,141 



Iodine 


128 


Dibasic fatty acid 

Solublization 

129 


m- (p-ani 1 ino ) -phenylazo- 
benzene sulfonic acid, 
potassium and sodium 
salts 

Absorbance studies 

119 


8-Ani 1 inonapthalene 
sulfonic acid (ANSA) 

Fluorescence 

97 


Nitrate ions 

UV spectroscopy 

97 

0-30 

4- [ (2-hydroxy-l- 
naphthalenyl )azo] - 
benzenesulfonic acid 
monosodium salt 
(tropacolin 000, 

Orange II) 

Light scattering. 
Diffusion, 

Viscosity and 
solublization 

1 


bis(p-chlorophenyl) acetic 130 

acid (p,p’-DDA); 1,1 (4- 
chloropheny 1 ) -2 , 2 , 2- 
trichloroethanol 


Cyclo- 

hexsine 


Meisenheimer complex; 131 

sodium 1 , 1-dime thoxy-2, 

4, 6-trinitrocyclohexa- 
dienylide 

0-50 amino acids; di-tri- IR, NMR 140,18 

and penta peptides 

18 azobis (2-aminopropane) ^^P-NMR, Kinetics 139 

(ABAP) hydrochloride- 
an auto-oxidation 
initiator; antioxidants- 
6-hydroxy-2,5,7,8 tetra- 
methylchromane-2-carboxy- 
late (Trolox) (water- 
soluble) and a-Tocopherol 
(organic solvent-soluble) 


13-14 


Proton NMR for 

CMC determination 

74 


Buffer solution 
(KH^PO^ and Na^HPO^) 

31 

P-NMR 

143 

0-19 

Cholesterol 

Temperature- 
dependent water 
solubilization and 
ternary phase- 
diagram studies. 

118 

5-30 

Phospholipase 

from Vespa 

Oriental is 


132 

0-15 


Optical 
birefingence, 
viscosity, light 
scattering 

92 

0-15 


NMR (^^P and ^H) 

110,94 

0-15 

ANSA 

Fluorescence 

97 

0-15 

KCl 

Electrical 

conductivity 

92 


0-15 


Nitrate ions 


UV spectroscopy 


97 


U-16 


Lholesterol 


118 


Temperature- 
dependent water 
solubilization 

^^P-NMR, IR, Laser 
Raman spectroscopy 


Carbon 0-15.8 Optical birefrin- 92 

tetra- gence, viscosity, 

chloride light scattering 


0-15.8 


and ^^P-NMR 

110,94 

0-15.8 

KCl 

Electrical 

conductivity 

92 

0-15.8 

Nitrate ions 

UV spectroscopy 

97 

0-15.8 

ANSA 

Fluorescence 

97 

0-18 


^H-and ^^P-NMR 

109 


Meisenheimer Complex 

Absorbance studies 

131 

5.5 or 

11 

Indole ; 1 -methyl indole; 
skatole 

NMR 

133 


NMR 


102 


n-Octane 3 Chloesterol 

n- 2.5 Cholesterol 

Dodecsuie 


Temperature 118 

dependent water 
solubilisation, 
ternary phase- 
diagram studies, 

NMR and IR 


Chloroform 


NMR 

142 


COCI 2 . 6 H 2 O 

Visible 

spectroscopy 

137 


Meisenheimer complex 


131 

Chloro- 18 

benzene 

an autoxidation 
initiator - ABAP; 
antioxidants-Trolox 
(water-soluble ) 
and a-Tocopherol 
(Organic solvent- 
soluble) 

Kinetics 

139 


0“ 

dichloro 

benzene 

0-10 


Ultracentri- 

fugation 

109 

Heptane 


Iodine 

Kinetics 

128 

Toluene 

0 


Osmometry, light 
scattering 

122 




viscosity, 

diffusion 

123 

Butanol 

0 


Light scattering 

124 

Ethanol 

0 


Light scattering 

124 

Methanol/ 

10-25 

phospholipase 


134 

Diethyl 

ether 


from Crotalus 
adamanteus 



Methanol/ 

pentanol/ 

octane 

11.1 

Acid phosphatase 
from wheat germ 


135 

Methanol/ 

pentanol/ 

octane 

13 

Peroxidase from 
horseradish 


136 

Hexanol 

0 


Light scattering 

124 

Methanol 

in 

benzene 

0 


Light scattering, 
viscosity and 
diffusion 

125 

Diethyl 

ether 

0-46 


Ultracentri- 

fugation 

106 


0-44 

COCI2.6H2O 

visible 

spectroscopy 

107 






conclusions have emerged: (i) only one water molecule becomes 
highly bound to the polar head group ; (ii) the microviscosity of 
the water environment is very high and falls off rapidly as a 
function of added water; (iii) the added water is solubilized in 
the polar core of the reverse micelle; (iv) water addition 
increases the motion of the polar head group but has no effect on 
the hydrocarbon chains. Proton studies on these systems have 
indicated additional relaxation mechanisms in addition to 
intramolecular dipole relaxation. Interestingly, the cycohexane 
system has shown unusual physical properties in the 
water-to-surf actant molar ratio range of 5 to 6 . 

During the same time 

have extensively studied the and NMR of egg-lecithin based 
micelles in benzene, and have come to conclusions essentially 
similar to those set forth above. They have also studied the 
effects of a wide variety of small and medium sized guest 
molecules such as NaCl, KCl, MgCl 2 , CaCl 2 , amino 

acids, peptides etc. (see Table 1.2) and have concluded that these 
guest molecules alter the distribution and organization of water 

molecules. They have also suggested the possible use of 

31 . . . ri43) 

P-relaxation parameters for indicating pH of the water pool' ' 

using a buffer solution. By measuring proton NMR chemical shift 

of egg-lecithin in benzene , the same group demonstrated the 

use of this NMR parameter for measuring the CMC of the reverse 

micelle. 

Most of these studies have been done on the benzene 
system. Very little attention has been paid to egg-lecithin in 
cyclohexane reverse micelles . In addition to work from our 


25 


laboratory, only one other group (see Table 1.2) has done 
31 

some ( P-NMR and water-solubilization) studies recently. 

Considering the abovementioned unusual features of the 
cyclohexane-based lecithin reverse micelles, and also taking into 
account that relatively little work has been done regarding the 
effects of water addition in this system, we have taken up this 
particular solvent for preparing egg-lecithin based reverse 
micelles. Further, we also been interested to prepare egg-lecithin 
based reverse micelles with a cosurfactant which may have 
potential uses in synthesizing biochemically important biomembrane 
models. We have attempted to characterize such system as far as 
possible by NMR (^H, and only. Preparation of one such 
system, and its physical properties, have been presented in 
Chapters II and III. 

Although the main bulk of our experiments has been done 
on egg-lecithin in cyclohexane, due to reasons spelt out in 
Chapter IV, we have taken up a deoxycholate reverse micelle for 
some multiple quantum NMR experiments. 

In Chapter V, we describe, in detail, the construction 
of a versatile magnetic field gradient controller for doing NMR 
measurements and measurements of self-dif fusion coefficients. 
Measurements of the water self-diffusion coefficients in 
lecithin/cyclohexane/H 20 have been compared with similar 
measurements in lecithin/CCl^/H 20 system. 
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CHAPTER II 

PREPARATION AND AND ^^P NMR STUDIES OF A NEW REVERSE MICELLE, 

LEC I THIN/DCC /C YCLOHEXANE/H^O 

II. 1 Introduction 

It has been reported recently that peptide bond 

( — ^NH — CO — ) formation occurs at the interfacial region of an 

Aerosol-OT based reverse micellar system . This result has been 

( 2 ) 

achieved by taking a 'peptide condensing agent' ' ' as cosurfactant 
while forming the reverse micelle. 

Based on the above, it occurred to us that if we could 
prepare a phospholipid-based reverse micellar system containing 
a 'peptide condensing agent' as cosurfactant, it may lead to useful 
information regarding how the freely diffusing globular proteins 
are both embedded in the bilayer (integral) and bound at the 
surface (peripheral) of the biological (or cell) membranes. ( see 
Figs. I. 1 and 1.3). 

We were successful in preparing a lecithin based reverse 
micelle with the very common 'peptide condensing agent'. 
Dicyclohexyl carbodiimide (DCC) , as cosurfactant. 

Q_n=c=n^Q 


Dicyclohexyl carbodiimide (DCC) 

Such a reverse micellar system should provide a very 
appropriate model for cell membranes, as it can also facilitate 
peptide bond formation and hold peptides and proteins at the 
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interfacial region. Although we had originally intended to check 
the formation of peptides inside this reverse micelle, our 
attention had to be diverted to other important aspect of our 
further studies. 

In this chapter, we report the preparation of the novel 
reverse micellar system, lecithin/DCC/cyclohexane/H20, and the 
results of our and NMR investigation on this system. 

II. 2 Experimental Details 

II. 2.1 Materials Used 

The extraction of chromatographically homogeneous egg 
lecithin was carried out using the method of Singleton et. al. . 
The extracted lecithin was stored at — 20°C in a 9 : 1 v/v mixture 
of chloroform and methanol. A single spot on a Silica Gel G 
(Merck) thin layer chromatogram confirmed the purity of the 
extracted lipid. 

Prior to using the lipid, the solvent was evaporated 
under vacuiim and the residual solid lipid was left in vacuum for 
about 6 hours. This solid was then dissolved in dry cyclohexane 
for forming reverse micelles for further characterization. The 
molecular weight of egg lecithin reported in most studies is in 
the range of 750-770. We have assumed a value of 760 in our 
studies . 

The cyclohexane used was of AnalaR grade, and this was 
further dried by distilling and refluxing it in the presence of 
added sodium wire for about three hours and then redistilling (b.p. 
81°C) . The dry solvent was stored over sodium in a bottle with a 
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tight stopper. 

The water used for making reverse micelles was deionised 
and double glass-distilled. 

II. 2. 2 Maximum Water Uptake Measurement 

As a first step towards preparing reverse micelles, 
solutions of 2%, 5%, 8% and 10% w/v lecithin in cyclohexane were 
prepared. By using 'Sigma Chemical' recalibrated microlitre 

pipettes, different volumes of water were added to these 
solutions. Then the samples were shaken thoroughly to ensure 

homogeneous mixing, allowed to stand for fifteen minutes and 
centrifuged at 800 rpm in a Clay-Adams table-top centrifuge and 
then inspected for phase separation. The amount of water added 
was increased in fixed steps. Based on the results of water 
uptake studies, the 5% w/v solution of lecithin in cyclohexane was 
chosen for finding out a suitable (lecithin : DCC) ratio to 
prepare a lecithin/DCC/Cyclohexane/Water system. On 

systematically varying (lecithin : DCC) molar ratio in the 5% 
lecithin solution of cyclohexane and subsequently adding water to 
form reverse micelles, we have chosen a (lecithin : DCC) molar 
ratio of 5 : 1 for further investigation. A schematic 

representation of the structure of this reverse micelle is given 
in Fig. II. 1. 

Since our main aim is to study the changes in properties 
of reverse micelle as a function of added water, we have kept the 
concentration of lecithin constant at 0.067 M , i.e., an 
approximately 5% w/v solution of lecithin in cyclohexane. Water 
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0 

I— J Polar region of Lecithin C-0-P-0-CH2-CH2-N(CH3)3] 

0 “ 

Polar region of DCC (-N=C=N-) 


Hydrophobic region of Lecithin 



Fig. II. 1 


Schematic representation of the structure of 
lecithin/DCC/cyclohexane/ water reverse micelle. 
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is then added in fixed increments. Here we define a quantity R as 
the molar concentration ratio of water to lecithin, namely. 


R = 


[water] 


[II. 1] 


[lecithin] 

This definition of R shall be maintained throughout this thesis. 
Further, the terms phosphatidyl choline (PC) and Lecithin shall be 
used interchangeably. 


II. 2. 3 NMR Studies 

(i) Chemical Shift Studies 

To probe the electronic environment surrounding the 

protons of the water molecule and phosphorous of lecithin 

indirectly, we have measured the chemical shifts of water and 

+ 

the trimethyl amino group (-NCCH^)^) by proton NMR, and the 
phosphorous present in the phosphatic head group of lecithin by 
NMR for the reverse micelles at different R values. 


(ii) Relaxation Studies 

To elucidate further the role of water in this reverse 

micelle, the H relaxation times of water molecules and those of 

31 • . 

the head group of lecithin, as well as the P relaxation times of 
the head group, have been measured at various R values. 

The longitudinal nuclear relaxation rate (1/T^) , also 
known as spin-lattice relaxation rate, leads us to getting 
information concerning the motion of molecules in solution^ . A 
change in molecular motion usually reflects a change in the 
structure of the medium and/or the interactions between the 
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molecules observed and their surroundings. 

All the proton NMR spectra were recorded on a BRUKER 

WP-80 FT NMR spectrometer operating at a frequency of 80 MHz. 

As this spectrometer needs a deuterated solvent for establishing 

magnetic field - to - frequency locking, we inserted a sealed D^O 

capillary tube in the 5 mm NMR sample tubes for all our spectral 

recordings. The common inversion recovery (180 °-t- 90°) pulse 

sequence was used for spin-lattice relaxation time (T^) 

measurements. The 90° pulse width for protons in this 

spectrometer was 5.2 ^s. 

31 

All P NMR spectra were recorded on a BRUKER WM-400 

31 

FT NMR spectrometer operating at a P frequency of 161.984 MHz. A 
capillary tube containing 80% phosphoric acid was used as a 
reference. T^ measurements were carried out once again by the 
inversion recovery method. The 90° pulse width in this case was 
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values were 


< '^2 = *’ 1/2 


measured from experimental spectra 
is the linewidth at half height) . 


II. 3 Results and Discussion 

From our water uptake-studies, we could determine that 
the maximum amount of water that could be solubilized in 5:1 
lecithin : DCC in cyclohexane solvent was 0.37 g of water per g of 
PC, or 15.8 mole of water per mole of PC. This value is only 
marginally more as compared to the value reported earlier from our 
laboratory for PC in cyclohexane ^ , which is 0.35 g of water per 
g of PC or 15.0 mole of water per mole of PC. This indicates that 
the addition of DCC does not increase the size of the reverse 
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micelle substantially to accomodate more water in the aqueous core 
as compared to the non-DCC containing counterpart. 

The slightly turbid and viscous solution formed between 
R values of 5 and 7 in DCC added reverse micellar system, which 
could be observed visually, were very similar to those observed 
earlier in our laboratory in PC/cyclohexane/water system. In the 
latter system, our studies have conclusively proved the 
viscosity/turbidity changes to be due to the change in the shape 
of the reverse micelle from spherical to tubular aqueous 
canals . 

Fig. II. 2 shows a proton NMR spectrum of lecithin and 
DCC in Cyclohexane. The observed peaks are assigned according to 
their chemical shift (5) values (relative to TMS) . The absence of 
a peak around 5 ppm in this spectrum clearly shows that there is 
no water present in the solution before adding water for forming 
the reverse micelles. Fig. II. 3 shows a representative spectrum 
for the R = 5 reverse micelle. 

On systematically increasing the water content of 

lecithin-DCC solutions from R = 2 to R = 15 for forming the 

various reverse micelles, and by plotting the observed chemical 

+ 

shift (5) values for water and -N(CH 2)3 protons in each case, one 
gets the plots of Fig. II. 4 and II. 5. 

From Fig. II. 4, the water proton resonance is upfield at 
lower water content compared to that of normal bulk water (S ~ 5) . 
This indicates a minimal amount of hydrogen bonding. As the water 
concentration is increased, there is an upfield shift of the water 
proton resonance followed by a downfield shift towards the normal 
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Fig. II. 3 80 MHz proton NMR spectrum of lecithin/DCC/ 

cyclohexane/water at R = 5. 
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NMR chemical shifts of water protons in 
lecithin/DCC/cyclohexane/H-0 as a function of R. 


Fig. II. 4 




Fig. II. 5 


NMR chemical shifts of — NCCHj)^ protons in 
lecithin/DCC/cyclohexane/H20 as a function of R. 
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Hh 

bulk water S value of around 5. However, the •-N(CH 2)2 proton 
resonance moves continuously upfield as the water concentration (R 
value) increases (see Fig. II. 5). 

The above trend in the water proton resonance arises 
because the chemical shifts of the -OH protons are strongly 
influenced by the following two factors ; (i) changes in the 

hydrogen-bonded state of water and (ii) changes in the chemical 
environment of the water molecule, i.e., the ratio of water bound 
to the lecithin head group to the free, bulk water in the inner 
'' pool ' . 

The trend of the water proton chemical shift observed by 

us at 80MHz is in the same direction as that reported by Shaw et. 

al. at 250 MHz for lecithin/CCl^/water system, Klose and 

{n\ 

Stelzner' ' at lOOMHz for lecithin/benzene/water system, and Kumar 

and Raghunathan ^ ® ^ for lecithin/cyclohexane/water system. 

31 

In P NMR, a slight increase in 5 value is observed 
upto around R = 4, after which there is a jump in 6 value at 
around R = 6, the chemical shift further increasing slowly as R is 
increased (see Fig. II. 6). 

Given the above trends in the ^H NMR chemical shifts for 

31 

water and -N(CH 2)2 , and P NMR chemical shifts for the 

phosphorous present in the lecithin polar head groups, we consider 

a plausible arrangement of lecithin head group in the reverse 

micelle. This arrangement derives support from a couple of 

(9) 31 

earlier findings, namely, those of Henderson et. al.' ‘ from P 
NMR studies of phospholipids in dry organic solvents and of 
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Fig. II. 6 


NMR chemical shifts of the head group phosphate 
in lecithin/DCC/cyclohexane/H20 as a function of R. 
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Pullman and Berthed from their quantum mechanical studies of 
energetically prefered conformations of the polar heads of 
phospholipids. Both the above studies favour the fact that 
the polar heads of the phospholipids are likely to adopt highly 
folded structures, as depicted in Fig. II. 7. 

For a water molecule (or molecules) getting 'bound^ to 
the head group of lecithin, the electrostatic interaction model 

H" 

with positively charged terminal trimethyl -amino group (-N(CH 2 ) 2 ) 
and negatively charged phosphate oxygens suggested by the above 
studies becomes significant as we shall see presently. 

As the water content is increased the most probable 
effect will be disruption, at least in part, of the folded 
arrangement of the head group and the production of a more 
'swollen' or extended conformation. Indeed, such a phenomenon has 
been envisaged by Pullman and Berthod^^^^ in connection with the 
theoretical model of the zwitterionic gamma-aminobutyric acid. 

Water addition reduces the electrostatic interaction 
between the positively charged trimethyl -amino group and the 
negatively charged phosphate oxygen atoms. Water is more likely 
to bind to the phosphate portion of the head group. The protons 
of water will thus form a weak bond with the phosphate group; 
likewise, the oxygen atoms of water will be in weak electrostatic 
interaction with -NCCH^)^. The suggested interaction of water 
molecule with the polar head group of lecithin is shown in Fig. 
II. 8(a). 

The possible picture for the electrostatic interaction 
of water molecule with DCC is given in Fig. II. 8(b). 
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ig.II.7 


Structure of folded polar head group of Lecithin. 
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twl No. » 

In addition to the mechanisms suggested above, there 
also exists the possibility of electrostatic interaction 
between the polar head groups of lecithin and the polar region 
(— N=C=N— ) of DCC. At the moment, however, we do not have 
experimental evidence to show possible electrostatic interactions 
of water molecules and lecithin polar heads with DCC. In our 
suggested model (Fig. II. 8(a)), the water is weakly hydrogen 
bonded with itself, or equivalently, the oxygen atoms are further 
removed from the hydrogen atoms ' ' . Such a weak hydrogen bonding 
would, of course, lead to an upfield shift of the -OH protons 
because it is shielded more thoroughly from the external magnetic 


field. This explains the initial upfield shift of the -OH proton 

resonance in Fig. II. 4. In fact, this shift is essentially very 

similar to that observed whenever structure breakers (i.e., 

electrolytes) are added to water Furthermore, the above 

trend in the observed chemical shifts strongly resembles the trend 

(14) 

in the PC/cyclohexane/water system studied in our laboratory ' 

At higher values of R,i.e., when the water is in excess, 

free water pools are formed with an increase in the average 'level 

of hydrogen bonding of water due to water-water H-bonds (i.e., 

self-association) , leading to downfield shift in the -OH 

resonance. This is in accordance with the result of Fig. II. 4. 

+ 

From Fig. II. 5, it is evident that -N(CH 2)3 proton 
resonance continuously moves upfield as a function of added watea: 
upto a certain value and then levels off. This, again, is quite 
expected in view of the added water electrostatically shielding 
the trimethylamino group. 


52 


The arguments given above for rationalising the -OH and 

Hh 

-NCCH^)^ proton chemical shifts are nicely complemented by our 

NMR chemical shifts of the lecithin head group presented below. 

• 31 

Fig. II. 6 shows the variation in the P NMR chemical 

shift as a function of added water. As R increases, the ^^P 
chemical shifts move downfield and, at R values between 4 and 6, 
there is a big jump in downfield shifting of chemical shift. As 
above, the explanation, of course, is that there will be a strong 
electrostatic interaction between the hydroxyl proton of the added 
water and the oxygen atoms of phosphate of lecithin. Such an 
interaction would be expected to deshield phosphorous and hence 
result in a downfield shift. The jump in chemical shift in the 
R-value range of 4 and 6 may be due to the onset of a change in 
structure of the reverse micelle and other physical properties 
such as viscosity etc. Such changes have been observed in a 
similar system studied earlier in our laboratory . 

The variation of the spin-lattice relaxation time, T^, 
of the water protons as a function of added water is shown in 
Fig. II. 9. One notices in this figure that T^ of the water 
protons increases roughly linearly from R = 1 to R = 14, except in 
the range R = 4 and 8, when there is a sharp 'peaking' in T^ 
values . 

We shall now analyze quantitatively the water proton 
spin-lattice relaxation times presented above in terms of 
microviscosity (i.e., viscosity of the reverse micelle at 
different smaller regions) values. It is well known that the 
relaxation times are related to molecular motion by the rotational 
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-orrelation time, x , which gives a measure of the mean residence 
-ime of the two proton nuclei in a given relative orientation, 
he correlation time, in turn, depends on the microviscosity. The 
pin-lattice relaxation time may be considered to be the result 
f a combination of intramolecular and intermolecular relaxation 
-ffects for any small molecule such as water containing two 
dentical spins. These two relaxation effects can be evaluated 
■heoretically. The relevant equations for intermolecular and 
ntramolecular relaxation effects for the water molecule, as 
ieveloped in the classic work of Bloembergen et al. , are : 


1 


T 

1 ( intra ) 


2y^h^I(I+l) 

5r^ 


4T. 


1+W^T^ 

c 


+ 


1 + 4 ( 4 )' 




[II. 3] 


1 


T 

1 (inter) 


2 4 2 

97r^r Tin N 
' 0 

5kT 


[II. 4] 


^ [II. 5] 

1 1 ( intra ) 1 ( Inter ) 

n eqns. II. 3 to II. 5 T, • , and T, , , are the 

1( intra) 1 (inter) 

ntramolecular and intermolecular parts of the dipolar mode of 

elaxation, y is the gyromagnetic ratio, r is the proton-proton 

istance in the water molecule, is the rotational correlation 

c 

ime, 0 ) is the resonance frequency of the water molecule, t) is the 
icroviscosity and is the volume concentration of spins. 

Substituting the value of obtained from eqn. II. 3 in 
he well known Debye-Stokes ecjuation 
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3 

_ 4Tn7a 
c 3kT 


[ 11 , 6 ] 


where a is the radius of the water molecule (approximated to be a 
sphere) , one can get ti, the microviscosity of the water molecules 
whose is known. 

We have analyzed the data of water protons in two 
ways : (i) under the simplifying assumption that the 
intramolecular relaxation mode alone is dominant in equations 
[II. 3] and [II. 6] and (ii) under the assumption that both intra- 
and intermolecular relaxation mechanisms are operative in 
equations [II. 3] - [II. 6]. The results of such calculations are 
presented in Table II. 1. 

We notice from Table II. 1 that the microviscosity (tj) of 
water for this reverse micelle is initially very high and falls 
off rapidly with the addition of water. The explanation for the 
downward tj trend is that, at lower R values, all water molecules 
will be bound to the polar head groups and show higher ri; when R 
increases the amount of free 'pool' water also increases, leading 
to decreased tj values. It reaches a lower limiting value at 
around R = 5. These trends are quite in agreement with the 
microviscosities of water in other micellar systems. For example, 
Shinitzky et al^^®^ reported microviscosities of 17-50 cP for 
water in cetyltrimethyl ammonium bromide micelles and Pownall and 
Smith calculated a microviscosity of 151 cP for water in 
hexadecyl trimethyl ammonium bromide micelles. 
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Table II. 1 Microviscosity of water environment in the 

Lecithin/DCC/Cyclohexane/H O reverse micelle. 


R 

(10 ^s) 

T according to [II. 3] 
and T), [II. 6] 

T according to [II. 5] 
and T] according to [11.61 

T^(10”^°S) 

V(cP) 

T^(10"^°S) 

TJ(CP) 

2 

83 

1.7S7 

53.6 

1.044 

31.3 

3 

133 

1.097 

32.9 

0.650 

19.5 

4 

181 

0.802 

24.0 

0.477 

14.3 

5 

382 

0.378 

11.3 

0.226 

6.8 

6 

328 

0.441 

13.2 

0.263 

7.9 

7 

312 

0.464 

13.9 

0.276 

8.3 

8 

249 

0.582 

17.4 

0.346 

10.4 

9 

298 

0.486 

14.6 

0.289 

8.7 

10 

328 

0.441 

13.2 

0.263 

7.9 

12 

352 

0.411 

12.3 

0.245 

7.3 

14 

360 

0.402 

12.0 

0.240 

7.2 


Constants used ; a = 1.50 X 10 m ; r = 1.52 X 10 m; 

h = 6.626 X 10~^^ J s; N = 6.798 X 10^^ ; 

o 

k = 1.38 X 10“^^ J/K/molecule; 

T = 307 K 

From Table II. 1, it is also seen that of the proton 
nuclei of water molecules decreases as a function of added water 
upto an R value of ~6. This shows that the mobility of the water 

molecules bound to the polar head groups of lecithin and DCC is 

. . . fl7) 

significantly reduced from that of normal bulk water' ' 

(t = 3 X 10 s) and hence showing higher x values. There is a 
c c 
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sharp decrease in r around R values 5-6. A similar trend was 
observed earlier in our laboratory for lecithin/ cyclohexane/ 

water system. Based on our earlier results^^'®^ ,a plausible 
reasoning for this sharp reduction in x may be due to the change 
in shape of the reverse micelle from spherical to tubular aqueous 
canals. The decrease in x is expected as x depends on the size 
and shape of the water pool. x decreases as R is increased due 
to the fact that once the hydration of polar head groups of 
surfactants by water molecules is complete, the further addition 
of water results in a free 'pool', which causes x^ to move 
gradually towards its normal bulk water value. 

The water molecules are only in contact with the polar 
head groups. Addition of water to the system increases the 
fluidity of the central aqueous pool and the interfacial region. 

4* 

As a result, the relaxation times of the —N (€ 112)2 protons would 
be expected to increase. This is found to be roughly the case as 
shown in Fig. II. 10. 

Fig. II. 11 shows the trend of spin-spin relaxation 
times, for the water protons as a function of R. T 2 values 

increase till R = 6 and then decrease. This leads us to suggest 
that some of the molecular organizational rearrangements necessary 
for the phase transition take place in this region . 

The Tj^ and T 2 values of water protons roughly follow the 
same trend, and the noticeable difference is that T 2 values are 
lower than the T^^. This difference shows that in addition to 
intramolecular dipole-dipole mode of relaxation there may be other 
modes of relaxation available for the water protons. The 
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rationalization for this is given below by calculating 
theoretically the values assuming that the intramolecular 
dipole-dipole relaxation mechanism is dominant and that the 
intemmolecular relaxation mechanism can be ignored. The relevant 
equation is^^*^^: 


1 _ y^h^I(I+l) 

T 6 

2(intra) SlT 


5t 


4T 


3t + 
c 


2 2 
1+(J T 

C 


l+4ca^T^-l 

c 


[II. 7] 


The rotational correlation time, obtained from eqn. 
[II. 3] is substituted in equation [II. 7] to calculate the T^ at 
various water concentrations (i.e., R values). The calculated and 
experimental values at different water concentrations are given 
in Table II. 2. 

It is seen from the Table II. 2 that the experimental T 2 
values are consistently much lower than those calculated. 
Ignoring such unlikely causes as poor magnetic field homogeneity, 
the results of Table II. 2 may be interpreted to mean that other 
relaxation mechanisms should also be considered for calculating T 2 
values, in addition to intramolecular dipole-dipole relaxation. 
Likely additional relaxation mechanisms are exchange narrowing, 
spin-rotational interaction and chemical shift anisotropy. We 
mention here that our experimental T 2 values have been evaluated 
from the absorption linewidths at half-maximum height. Slight 
inhomogeneities in the magnetic field would cause additional line 
broadening and thereby result in a decrease in the apparent 
relaxation time. However, we are aware of the fact that this can 
be eliminated by measuring the signal line width at two different 
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magnetic field strengths, because the field independence of 
linewidth is a classical test for magnetic field homogeneity . 

Table II. 2 Calculated and Experimental values of water 

protons in Lecithin/DCC/Cyclohexane/H^O reverse 
micelle. 


calculated from obtained from spectra 

egn. [II. 7], [T = ]/(in seconds) 

(in seconds) ^ * 1/2 


2 

0.068 

0.037 

3 

0.110 

0.049 

4 

0.150 

0.063 

5 

0.318 

0.098 

6 

0.273 

0.083 

7 

0.260 

0.071 

8 

0.207 

0.062 

9 

0.248 

0.060 

10 

0.273 

0.062 

12 

0.293 

0.057 


31 

Fig. 11.12 and 11.13 show the variation of P linewidth 
and values, respectively, for the phosphorous in the choline 
head group of the reverse micellar system. 

From Fig. II. 13 it is clear that, when no water is 
added, lecithin molecules show a high value of about 1.4 s. 
This may be ascribed to free— tumbling lecithin molecule which 
occur in solutions free, non-aggregated monomers. When water is 
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added for forming reverse micelles, reduces drastically. 

This is due to the electrostatic interaction between the hydrogen 
of the water and the oxygen of the phosphate, which restricts the 
motion of phosphate head groups, thereby leading to a decrease in 
T^. Further addition of water enhances the "structuring” of more 
water molecules around polar head groups of lecithin, and 
consequent reduction in the motion of phosphate groups. After all 

the phosphate groups have been hydrated, further addition of water 

3 1 

does not change the of P sxibstantially, as is evident from 
the final 'plateau' value of -760 ms (Fig. 11.13). 

Interestingly, the decreasing trend of the phosphorous 
'T^ vs R' values observed for lecithin/DCC/Cyclohexane/water, 
which we have justified in terms of an electrostatic model, is 
quite contrary to the general trend reported for similar lecithin 
based reverse micelles . It is conceivable that the presence 
of the co-surfactant (DCC) in our system plays a distinct role in 
altering the head-group dynamics. However, we shall not pursue 
this point further here. 

. 31 

Fig. 11.12 shows a pronounced increase in the P 
linewidth as R increases. This is only to be expected, 
considering that there is a 'visible' increase in the bulk 
viscosity of the reverse micelle as R increases, thereby heralding 
an overall slow-down in micellar motion and/or possible 
structuring at the micellar interface due to hydration upto R 
value of -8. 


We have described in this chapter our extensive NMR work 
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to characterize the novel reverse micelles formed by 
lecithin/DCC/ cyclohexane/water . What is perhaps crucial to a more 
thorough understanding of the temperature and composition range of 
the 'reverse micellar phase (s)' of this system is the construction 
and study of the phase diagram of this multicomponent system. 
This we have not attempted in our work. 

II. 4 Conclusions 

(i) The novel reverse micellar system and lecithin/DCC/ 
cyclohexane/ water, which may be of much potential use in 
peptide chemistry, has been prepared and characterized by 
NMR . 

(ii) A -5% w/v solution of lecithin in cyclohexane containing 
DCC in 1:5 mole ratio with lecithin solubilises about 0.37 
g of water per 1 g of lecithin. During the course of water 
addition, a slightly turbid and highly viscous phase is 
observed visually. 

(iii) Chemical shift values of water protons indicate that the 

water molecules are weakly hydrogen bonded at lower R 
31 

values. P chemical shifts of phosphate head group show 
evidence for an electrostatic interaction between water 
protons and the oxygen of the phosphate head groups of 
lecithin. 

(vi) The rotational correlation time (t ) , calculated from 
proton NMR spin-lattice relaxation time for the reverse 
micellar water protons, shows that at lower R values it is 
higher (as compared to normal water) and decreases toward 
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the normal water value as R is increased. 

(v) The microviscosity of the water environment is very high at 

low water concentrations, and decreases rapidly as a 
function of added water. 

(vi) Analysis of proton values of water suggests contributions 

from other relaxation processes such as exchange narrowing, 

spin-rotational and chemical shift anisotropy, in addition 

to the dominant intramolecular relaxation process. 

. 3 1 

(vii) The P T^ values show a decreasing trend as R is 
increased; this is in contrast to reported trends observed 
in other lecithin based reverse micelles, and suggests that 
DCC may be playing a role in altering the head-group 
dynamics . 
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CHAPTER Hi 

DEUTERIUM NMR LINESHAPE ANALYSIS OF REVERSE MICELLAR POOL WATER 

III.l Introduction 

In recent years, deuterium magnetic resonance 
(DMR) spectroscopy has been increasingly widely exploited as a 

means of investigating molecular dynamics in lipid 

M— 10) ^11-17) 

membranes ' ‘ and surfactant aggregates ' ' . The great 

utility of DMR derives from the fact that deuterium lineshapes and 

relaxation times are determined by a single well-defined 

interaction, the quadrupole coupling, which permits one to study 

internal dynamics spanning a wide range of rates between about 
3 ""1 8 *"1 

10 s and 10 s in a relatively direct manner. 

With a view to employ DMR for probing the dynamics of 
reverse micellar 'pool water' , we have chosen two systems for our 
study, namely, (i) Lecithin/Cyclohexane/D 20 (ii) 

Lecithin/ Dicyclohexyl Carbodiimide /Cyclohexane 
different temperatures. In similar contemporary experiments in 
aqueous media, DMR studies have been extensive, showing that water 
near model membranes and biomembranes possesses significantly more 
structure than, and dynamically very different from, aqueous phase 
'bulk 'water. 

In this chapter, we shall first briefly present the 
general aspects of DMR spectroscopy. This will be followed by a 
description of the 'quadrupolar echo' method and the spin system's 
response to the quadrupolar echo pulse sequence. Our 

computational method for DMR lineshape simulation, and the 
interpretation of our experimental results on the two reverse 



micellar systems by means of these simulation procedures, shall 
constitute the rest of this chapter. 
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III. 2 DMR Spectroscopy 

2.1 General Theory 

In addition to its magnetic dipole moment, a nucleus 
with I a 1 possesses an electric quadrupole moment, eQ, which has 
its origin in a non-spherically symmetric nuclear charge 
distribution^^® 20) ^ This moment will interact with eq, the 
nonuniform electric field gradient (efg) , which is produced by the 
surrounding valence electrons and other nuclei, at each nucleus to 
give a series of quantized energy levels, even in the absence of 
an external magnetic field. 

Consequently, the total Hamiltonian for the nucleus in 
the presence of an applied magnetic field is 


total z 




[iii.i] 


where is the Zeeman Hamiltonian representing the interaction 
between the applied magnetic field and the nuclear spin (I), and 
Hq the quadmpolar Hamiltonian due to the interactions between 
quadrupole moment associated with spin I and efg existing at the 
nuclear site. In eqn. III.I, all other internal Hamiltonians 
except Hq have been momentarily suppressed for simplicity. 

The energy levels resulting from a first order 
perturbation solution^^® of the total Hamiltonian, eqn. 
III.I, are given (in angular frequency units) by 
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Then, 


E 


la 



+ E 


( 1 ) 

m 


[III. 2 ] 


E. 


m 


= - 0) iti 

o 


‘"q 

r 2 

3cos 0—1 

3 

2 


+ 5" 


sin^e 


[ 


cos2<p\ 3m -1(1+1) , [III. 3] 


]' 


where is the Larmor frequency, m is the magnetic quantum number 
in a representation where the z-component of I is diagonal, e and 
^ are the usual polar and azimuthal angles specifying the magnetic 
field direction in relation to the principal coordinate system of 
the efg and is given by 


3 e^qo 

Q 8 hl(21-l) * 


[III. 4] 


The term 173 in eqn. III. 3 is the asymmetry parameter, defined by 


V 


a 


V. 


XX 



zz 


/ 


[III. 5] 


where V. . are the components of efg tensor defined such that [V I 
- I^yyl “ * Since the electric field surrounding the nucleus 
is produced only by charges external to the nucleus, the Laplace 
equation, given below, holds : 


V 

XX 


+ V. 


yy 


+ V 


zz 


0 


Thus 1 ) ranges from 0 to 1. 

cl 


[III. 6 ] 


The term 


e^qQ 


in eqn. III. 4, where eq = V 


zz 


d^v 

dz^ 


is referred 


to as the nuclear quadrupole coupling constant. 



72 


2 . 2 DMR in D^O 

Since the deuterium nucleus has I = 1, m will assume 
values of 1, 0, and -l. The DMR spectrum arising from a 
particular deuteron will consist of two sharp peaks centred about 
0 }^ and separated by 

Ao) = Wq (3cos^e-l) + Tj^ Qq sin^e cos2(^ . [III. 7] 

One notices from eqn. III. 7 that, for cases where the 
efg at a deuteron site has cylindrical symmetry, the asymmetry 
parameter is practically zero, reducing eqn. III. 7 to the simpler 
form 

Aw = cJq (3cos^e--l) . [III. 8] 

For water deuterons in our systems, although the major 
contribution to the efg is from the intramolecular 0-D bond, there 
can be other contributions of intermolecular origin such as the 
charge distribution in the vicinity of polar head group and other 
interactions. Furthermore, hydrogen bonding can reduce the 
contribution of the 0-D bond from 312 KHz (as measured in the gas 
phase) to 213 KHz (22-24) ^ such reasons could be different 
for different sites. Also, chemical exchange can take place 
between nuclei in different environments. If the exchange rate is 
much faster than the quadrupolar splitting (Aw) , the observed 
splitting will be a weighted average over the different sites, 
given by (25-27) 
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Aw = 


I Pi 


1 ~Q 


[III. 9] 


where is the probability (defined by the fraction of nuclei at 

site i) of a nucleus at being site i with a characteristic 
The exchange rate will also vary with temperature. 

In the extreme limit, if isotropic tumbling motion of 
the quadrupolar nuclei is fast compared to the static splitting 
frequency, the angular term in eqn. III. 8 is averaged to zero and 
no splitting is observed. On the other hand, the motion in a 
lamellar or liquid crystalline phase would in general be expected 
to be anisotropic, and a splitting of the first order spectrum 
will usually be observed. 

The energy levels and spectra corresponding to single 
cyrstals and polycrystalline samples of deuterium are given in 
Fig. III.l. 


In summary, then, Aw data and exchange of deuterons 
between different environments can provide useful dynamic 
informaion through DMR spectra. 


III. 2. 3 System Response to the Quadrtipolar Echo Sequence : 

2.3.1 Theoretical Background 

In this subsection, we review the general procedure used 
to follow the time evolution of a spin system. 

(i) Density 0perat3.or and Time Evolution of the Spin System 
The dynamics of coupled spins are often understood 
elegantly with the help of a quantum mechanical formalism which 
characterises the state of the system by an operator called the 
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Hami Itonian Quad rupolar 

Hamiltonians 


(a) Energy levels 


tjQOCos^e-i) 



(c) When tumbling motion of (d) Anisotropic motion, first order spectrum 

Quadrupolar nuclei is fast 
compared to Aco 


Fig. III. 1 


Energy levels (a) and DMR spectra corresponding to 
single crystals (b) of a spin I - 1 system (and) at 
two different motional states ( (c) and (d) ) 
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density operator, 

For a multilevel spin system with a complete set of 

eigenstates |^j^(t)>, the state of the system can be written as 

|(;r(t)> = Y. C. |^,(t)>. Then, a mixture of independently prepared 
k 

states |^/^^(t)> (n = 1,2,...) with statistical weights can be 
completely described by the density operator p(t) as 

P(t) = I X [III. 10] 

n 

The attributes of this operator and its matrix elements 

are : 

(a) it is Hermit ian, that is, 

«^i|p(t)|^j> = <^j|p(t) 1^.>* , [III. 11] 

(b) the diagonal element p(t)^j^ of the density matrix is equal to 
the probability of finding the system in the eigenstate |^j^(t)>, 

(c) the off-diagonal element indicates a 'coherence 
superposition' of eigenstates, simply called 'coherence', and 

(d) the expectation value <Q> of an arbitrary observable 
Hermitian operator Q is given by 

<Q> = Tr {Q p(t)} [III. 12] 

The time development of the spin system under various 
interactions defined by the Hamiltonian, K, can then be studied 
using the equation of motion of the density matrix. 


- i [«, P(t)] 


[III. 13] 
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This differential equation is familiarly known as the 
Liouville~von Neumann equation (30 ) , and has the solution 

p(t) = U(t) p(o) U{t)“^, [III. 14] 
where U(t) = T exp - i J (t') dt' j- . [III. 15] 

U(t) is variously called the 'Liouville operator' and 
'Dyson Expression' ^ and T is the Dyson time-ordering operator 
(31) which has the property of ordering the exponential terms 
appearing in eqn. III. 15 with increasing time from right to left. 

In magnetic resonance, the time evolution of spin 
population can be thought of as rotation of spin vectors in a 
(21+1) dimensional space. Exponential operators of eqn. 111.15 
are used to explain such a time evolution. The mathematical 
physics of exponential operators may be found in standard 

(ii) Nuclear Spin Interaction Hamiltonians 

Various interactions of a nuclear spin system in NMR 
experiments can be described by a laboratory frame 
Hamiltonian ^ ^ 


«L = 


n . + 

ext 


n 


int 


[III. 16] 


where and represent the nuclear spin interactions 

between external fields (®o' ®1^ internal fields. 
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respectively. consists of two terms, namely, 


where 


«ext = «rf 

«z = 


= -20)^ (t) cos (cjt + <p) 


[111. 17] 

[111. 18] 

[111. 19] 


^rf' Hamiltonian, describes the interaction of 
the spin system with a linearly oscillating rf field perpendicular 
to the static field. w^(t) equals rB^(t)/2, where B^(t) is the rf 
magnetic field amplitude, o and 4> 3^re the frequency and phase of 
the rf. Typically, the rf field is applied in pulses that are 
ideally square, so that tj^(t) is piecewise-constant, taking on the 
values 0 and tj° only. will be referred to as the nominal rf 
amplitude. Also, typically, w is constant and 4> is 
piecewise-constant . 

In eqn. III. 16, describes interactions of the spin 
system with internal fields, which include dipole-dipole 
interaction (K^) , J-coupling interaction (^j) / spin-rotation 
interaction (^gj^) «ind quadrupolar interaction (M^) . Since our 
present focus is only on the secular part of (assximing that 
the field gradients have axial symmetry) may be written as 


= 


? “Qi [^li - (V3) 1 ?] 


and 




"Qi 


^ Ocos^fl.-l) 


[III. 20] 

[ 111 . 21 ] 


where, as defined earlier, eQ and eq are respectively the 
quadrupole moment and the efg along the unique directions for the 



spin. 6^ is the angle between the static field direction and 
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the unique field gradient direction. 

The dominant term in Eqns. III. 16 and 17 is the 

interaction with the static field. All the contributions to ?{. . 

mx. 

are usually at least 100 times smaller. Thus, to a good 
approximation, the eigenstates of when the rf is turned off 
are eigenstates of I . They occur in Zeeman manifolds 

Z 

characterized by the eigenvalue m of and separated in energy by 

increments of roughly When the rf field is turned on, it 

induces a coherent mixing of the eigenstates if w is about equal 

to (j) . 
o 

(iii) Rotating Frame Transformation 

The well-characterized, and therefore uninformative, 
interaction with the static field can be largely removed, the 
oscillating time dependence of the rf can be made to vanish, and 
the important parts of ^int made apparent by a 

transformation to a new frame of reference called the 'rotating 

(18 , 19 , 32a, 33 ) type of transformation from one frame 

f34 35) 

to another commonly known as interaction frame ' " ^ in quantum 

mechanical systems is described very briefly from a general point 
of view. 

A change of reference frame is accomplished by 
performing a unitary transformation on an equation of motion in a 
representation, in which time dependence exists in both operators 
and state vectors, either to remove or to take account of selected 
portions of the Hamiltonian from the description of motion of the 


system. 
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Let us consider a general case with the total 
Hamiltonian consisting of two terms and , one or both of 
which may be time dependent. The equation of motion of the system 
is 


id| 

dt 


+ H ^) |^> 


If we perform a unitary transformation, then 


[III. 22] 


|^/>(t)> = |$;(t) > 


[III. 23] 


where |iA(t)> is the state of the system at time t in the original 
frame and |^(t)> is the same state at time t in the new frame. 

Assuming that is determined by via the 

differential equation 


idU^ 

"dF" 


= KiUi 


[III. 24] 


,- 1 . 


with the boundary condition that (t =0) = (t=0) = 1 eqns. 

III. 22-24 give 


dt 


n.+n. («,+«,) u-|>> 

|^> H ^ 


[III. 25] 


dt 


idu- ^ n u |^> 


[III. 26] 


Differentiating egn. III. 23, 
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dt 




+ U 


d|iA> 

1 dt 


[III. 27] 


comparing III. 2 6 and III. 27 


1 dt I 


[III. 28] 


Multiplying III. 28 on the left by iU^^ yields the time dependent 
Schrodinger equation for 

i^^ = U"^ J<2 ^i\^> ± [III. 29] 

where = u“^ 

It is seen in eqn. III. 29 that |^i> no longer explicitly depends on 
the interaction used to perform the transformation. This is a 
transformation to the 'interaction frame' of 

In particular, the type of unitary transformation used 
to transform the spin system description from laboratory frame to 
the rotating frame of the Zeeman interaction is 


U(t)rot = ®xp (-i t) [III. 30] 

It is important to note that in this new frame, the Zeeman field 
is removed from our description of the system, while the rf field 
becomes stationary. 


RF Hamiltonian in the Rotating Frame 


The rf Hamiltonian in the laboratory frame is 
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= -2ca^I^cos (wt + (^) 


[III-31] 


The rotating frame RF Hamiltonian is obtained by performing the 

unitary transformation defined by eqn. III. 30, in which n = 

z 

-rhl^Bo = = rB^. Thus, the rotating frame rf 

Hamiltonian is 


*r£ - “^ot «r£ “rot 


Using the value, and neglecting high- frequency terms. 





I cos^ + I sin^ 

. ^ y . 


[III. 32] 


Quadrupolar Hamiltonian in the Rotating Frame 

The quadrupolar Hamiltonian (Kq) for the axially 
symmetric field gradient case is 



This does not change its form after a unitary transformation to 
the rotating frame as the operators present in the secular part of 
always commute with i.e., Therefore, 


% ' “roi «Q “rot 




[III. 33] 
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(iv) Density Matrix for 1=1 case 

The evolution of the density matrix in pulse experiments 
on spin I =1 can be calculated by expressing it in terms of an 
orthogonal basis set of matrices. Many such orthogonal basis sets 
are available in the literature . We have chosen the 
orthogonal basis set of matrices proposed by Vega and Luz^^^^, as 
their method is well-suited to our problem. Their basis set 
matrices are Hermitian and therefore give a real physical 
significance to the spin states they represent. Also, these 
matrices have convenient commutation relations with the operators 
that occur in the Hamiltonians of our interest. 


A system of non interacting 1 = 1 spins in an external 
magnetic field can be defined by the spin density matrix, p, 
having a dimension of (21+1) X (21+1) and a total number of 
(21+1) ^-1 traceless independent Hermitian operators, • This 

density matrix can be expressed as 

N 

Pit) = I a (t) A + a 1 [III. 34] 

n=i 

where 1 is a unit matrix and a^(t) time-dependent coefficient. 
These can be obtained by solving the equation of motion for p(t) . 

For our problem, the set of eight independent traceless 
operators required are ly, for spin 1 = 1 and the following 
additional operators defined as 





Vy ’^y^x 



[III. 35] 
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Alternative linear combinations of Q and K are : 

z z 



[III. 36] 


[III. 37] 


The nineth operator is the 3X3 unit matrix, U. 

The matrix representation of the operators defined above 
are given in Table III.l 

The transformation relation existing between the three sets of 
second rank operators (Eqns. III. 36, 37) are 


= I (-°K + V 

Oj = 5 (-«K - •ic) 

J'i = I <-='«K - ’St> 

•'j = 5 <='«K - 


[III. 38] 


where i,j,k = x,y,z or their cyclic permutation. 
Different commutation relations existing between the basis-set 
operators and the commonly encountered spin-interaction operators 
are given in Table III. 2. 



Table III.l 


Matrix representations of I =1 spin operators (from 
reference 38) : 
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Table III. 2 
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Commutation relation between 1=1 spin operators 
and commonly encountered spin-interaction operators 
(adopted from ref. 38) 



i[K^,l^] = 2J^ 

i[Q,,I^] = 0 


i[l^,Iy] = -I, 
i[J^,Iy] = J, 

= -2Ky 

i[Iy^V = 0 

= ^y 

= "^y 
= -2K^ 

it^x'Qz^ = ^z 

iCIy^Q^] = -Jy 
iEl^^Qz] = 0 

i[Q^,Q,] = 0 


iCJ,,ly] = 

i[Ky,ly] = 2Jy 

i[Qy,Iy] = 0 

iCIy.Ij.] = -Ix 

iCJy'^z^ = ""x 

i[K^,I^] = 2J^ 

itQz'^z^ = 0 
iCJx^Qz^ = -Ix 

i[Jy.Qz3 = ly 
iCJz.Q^] = 0 

iCKj^Q^] = 0 

(j = X, y, z) 


III 3. DMR Measurements 

3.1 Solid Echo or Quadrupolar Echo Pulse Sequence : 

In conventional FT NMR, spectra are obtained by applying 
a 90° pulse and then Fourier transforming the free induction decay 
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(FID) . Application of this method to lineshape investigation in 

deuterated viscous liquids and solids does not give a faithful or 

true lineshapes as the very early part of the FID, containing the 

information about broad components of the line, is lost due to 

'dead time^ effects of the rf receiver and due to pulse length, 

recovery times, magnetoacoustic ringing, etc. To circumvent dead 

(39) 

time problems, Powles and Mansfield have proposed a pulse 

sequence for a pair of spin -1/2 nuclei which has come to be known 
as the 'solid echo'. The sequence is 


Tr/2 1^/2 1 _ acquire 

(Jo ( J 90° 

where the subscripts denote the phase of the pulse. 

(40) 

Later, Davis et al. showed that for the deuteron 

the same sequence gave a perfect refocussing for the quadrupolar 
Hamiltonian. Hereafter, this sequence will be referred to as 
quadrupolar echo sequence throughtout this thesis. 

3.2 Quadrupolar Echo Pulse Sequence and the Spin System's 

Response : 

The quadrupolar echo pulse sequence is illustrated in 
Fig. III. 2. 

In this sequence, z is the variable interval between 
pulses, the superscripts - and + indicating 'just before' and 
'just after' the pulse, respectively. Signal (FID) acquisition 
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starts from 2x and recorded for time 't' . The response of the 
spin system to this sequence is described below. 

At thermal equilibrium, when the high-temperature, 
high-field approximation becomes valid, the density matrix in the 
laboratory frame can be written as 


P 


eg 


P(o) 



P 


eq 


- bl 


z 


[III. 39] 


where b is equal to z-component of the total 
spin operator. This density matrix, on transformation to rotating 
frame using eqn.III.30, becomes 


P(o) - p(o) 

p(o) = -bl^ [III. 40] 

Eqn.III.40 represents the spin system at time 0~ (see Fig- III. 2) ; 
thus 

p(0“) = p(0). [III. 41] 

At the instant 0~, rf pulse is applied to the spin system for the 
duration t^. The resulting rf pulse perturbation can be explained 
in terms of the unitary operator. Up, defined as 

Up (tp) = exp (-i tp). [III. 42] 

Density matrix of the system immediately after (^r/S)^ pulse (at 
the instant of O^) is 
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p(0‘^) = Up p(0~) Up^ 

Using Table III. 2, and taking the rf phase factor (0) arbitrarily 
as zero (as it is the first pulse) , the above eqn. may be written 
as 

PCo'*') = b [I^ cos 0 )^ tp - ly sin tp] [III. 43] 

We know that the pulse tip angle, 6, is e = yH.t = ti/2 and 

I P 

substituting o)^ = yH^, 

0 = tJ^tp = 90° [III. 44] 

Putting [III. 44] in [III. 43], we get 

PCo"^) = - bly [III. 45] 

After this pulse, the spin system evolves freely under the 
influence of Neglecting the relaxation effect for the time 

being and assuming 'on resonance^ condition, the density matrix, 
P(0^)f just before the (n/2)^ pulse becomes 

~ / ~ ,^+. i^«T 

p(T ) = e Q p(0 ) e Q 

Using Table III. 2, it can be shown that 

p(T~) = b ly cos WqT + Jy sin [III. 46] 

Just after the second pulse, the density matrix is 
P(t''’) = exp (“iji^^t^) p(r ) exp (tift^^tp) , 
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P(t'*’) = b [- ly cos WqT - jy sin [III. 47] 

Following the {n/2)^ pulse, the density matrix, p{z^), evolves 
under Hq again, and spins along y refocus at time 2x according to 

P(2t) = exp pCt"^) exp(i«QT) 

= -bly ^sin^QgT + cos^ 

p(2x) = -bly [III. 48] 

We notice the interesting fact that the density matrix 
at time 2t is the same as the density matrix immediately after the 
first pulse, thereby showing that the inherent detector 'dead 
time' problem is completely removed by employing this pulse 
sequence for the excitation of the spin system. 

3.3 Computational Method for Lineshape Calculation ; 

The response of the spin system immediately after the 

pulses, evolution under following pulses, and at 2x, i.e. eqns. 

^ (38) 

III. 45, 46, 47 and 48, prompted Vega and Luz to write the 

evolution of density matrix during the quadrupolar echo 
experiments in terms of just two components of the basis set as 

P(t) = y(t) ly + C(t) Jy [III. 49] 

This p(t) under can be calculated from the usual equation of 
motion. 
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= i [p(t), , [III. 50] 

that is, 

= i [y(t) ly + C(t) Jy, qJ . 

Using the commutation relations of Table III. 2, 

d[y(t) I + C(t) J ] 

" Jy + 

The expectation value of 1^ and is calculated as 


d Tr [y(t) ly. ly + C(t) Jy. ly] = 

Wq Tr [-y(t) Jy. ly + C(t) I y ly] [m.Sl] 

d Tr [y(t) ly. Jy + C(t) Jy. Jy] = 

Wq Tr [-y(t) Jy. Jy + C(t) ly. Jy] [lll.52] 

Using Table III.l, eqns. III. 51 and 52 become 

C(t) [III. 53] 

dt Q ' ' 

= -Ug y(t) [111.54] 

[III. 53] + i [III. 54] gives, 

d [y(t)-hiC( t) ^ . [yjt) + iC(t)] [III. 55] 

dt Q L -1 


(38) 

Vega and Luz define y(t) + iC(t) as an effective complex 

maanetization. Gft] . namelv. 
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where y(t) is the expectation value of the usual magnetization, 
MySlong y— axis and C(t) is an NMR unobservable coherence such as 
Jy etc. The effect of ( 71 / 2 )^, pulse on G(t) is proved to be simply 
its complex conjugate, i.e.. 


G(t) 




G*(t) 


[III. 57] 


Eqn. III. 55 may now be written as 


dG(t) 

dt 


-i Wq G(t) . 


[III. 58] 


It has been shown that G(t) is following a Bloch type of equation 
under guadrupolar Hamiltonian. 

Using the equations described above, a method is derived 
for simulating deuterixim spectra obtained using quadrupolar echo 
sequence. For this purpose, we assume that there are only two 
types of sites for D 2 O molecules in the reverse micelles; one is 
D 2 O bound to the polar head group in the interfacial region and 
commonly called as 'bound water' , and another is D 2 O which are 
present in the interior of the reverse micelles and called as 
'free water' as it is similar to normal bulk water (See Fig. 
III. 3) . 

As explained in Sect. 2.2, these two types of water will 
have different and T*. There is a continuous exchange of water 
molecules taking place between 'bound water' sites and 'free 
water' sites and also the two protons of a water molecule flip 
about its C axis. With this model in mind, we have derived the 
necessary equations for simulating experimental ^^ 2 ^ spectra 




water' molecules as site 1 water molecules and 'free water' 
molecules as site 2 water molecules. 

Incorporating exchange and relaxation tejrms for these 
two types of water sites, eqn. III. 58 becomes 


dG^(t) 

= - [“qi + "21 

1 1 

a 

(t) 

+ K^2S<'^> 

[III. 59] 

dt 

CN 

^1 

dG2 (t) 

= - [i“Q2 + Ki2 

1 

G 

(t) 

+ Kj^G^Ct) 

[III. 60] 

dt 

T2<1) . 


In the 

above eqns . , K ^2 

H 

<N 

tl 

= 

K, 

the exchange 

rate. 


subscripts 1 and 2 correspond to site 1 and site 2, respectively, 
and is already defined in eqn. III. 21. Eqns. III. 59 and 60 
can be combined as below : 


dG(t) 
dt 


= - RG(t) , 


[III. 61] 


where G(t) is a column vector having G^{t) and G 2 (t) as its 
components. R is a matrix whose elements are 

Rli = i + K + 5^, 

Rjj = -K, i j and i = j = 1/2 

Immediately after the first {n/2)^ pulse in the quadrupolar echo 
sequence, all magnetizations corresponding to both sites Ire along 
-y direction (ref. eqn. III. 45) with magnitudes proportional to 
their fractional populations 


GCO"^) = P , 


[III. 62] 


where P is a real column vector defined as 


in which p.'s are the fractional populations of site i 

P2J ^ 
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Following the first pulse, the system evolves accoriding 
to eqn. 1X1. 61, and just before the second pulse we have 

dG(T”) -n “v 

“dt = - R G(t: ) . [111.63] 

G(t ) is obtained by integrating the above equation as 
shown below : 


dG(T ) _ 
G(t") 


R dt 




Then 


G(t ) = exp (-Rt) G(0 ) = exp (~Rt) P 


[111.64] 


The spin system just after the (n/2)^ pulse, according to ecpi. 
111. 57, is 

G(t:''') = G(t")* 

= exp (-Rt)*P [111.65] 


Similarly, the complex magnetization at 2 t is 

G(2t) = exp (-Rt) G(t^) = exp (-Rt) exp (-Rt)*P [111.66] 

Here we define J as J = G (2 t) where J is the colvimn vector 
having components j^ = G^(2 t) and ~ G2(2t). 

The total echo amplitude, E, is then calculated by 
summing j^ and j 2 
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E = 1. J = 1. exp (-Rx) exp (-Rx)* P 


[III. 67] 


where 1 is a row vector with its two components equal to unity. 

In the actual quadrupolar echo experiment, FID is 
recorded starting from 2x and for a time 't'. Therefore, 


FID(t) = Re 
= Re 


-|l.G(2x+t) j 
-|l.exp(-Rx) 


.J 




The absorption spectrum is obtained 
transformation on FID(t) , 


by 


doing 


I(w) = 


00 

[ FID(t) cos £Jt dt 


o*' 


[III. 68] 
the cosine 


[III. 69] 


The exponential matrices appearing in the above 
equations can be very conveniently expressed ' ' for our 
calculation in terms of a diagonal matrix, containing complex 
eigenvalues of R, and a matrix S of its eigen vectors 


Rx = S e^^ 


e-^t = S e-^^ S-" 
where s“^RS = A 


Substituting [III. 70] in [III. 68] 


[III. 70] 


FID(t) = Re 
= Re 


{l • S J } 


} 


[III. 71] 
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and the eqn. III. 69 becomes. 


I(u) = Re 


{e — 

'^ijk 0) 




^ 




[III. 72] 


The I(tJ) of Eqn. III. 72 is calculated for all possible O-D bond 

orientations with respect to magnetic field, and the overall 

spectrum is obtained by summing I (a) over all the possible 

orientations. A fast Monte Carlo algorithm based on systematic 
. (42 43) 

sampling ' 'has been previously developed in our 
(44 45) 

laboratory' ' ' ; this algorithm has been used in the present 

work for the overall simulation. Eq. III. 72 is implemented by 
means of a programme written in FORTRAN 77 to simulate all our 
experimental deuterium magnetic resonance spectra. These results 
will be discussed further in the next section. The programme in 
its entirety is listed in Appendix I. This programme was executed 
on the HP 9000 computer system of I.I.T. Kanpur. 


III. 4 Results and Discussion 

All deuterium NMR spectra have been recorded at 61.25 
MHz on a JEOL JNM-GSX400 FT NMR spectrometer, by using the 
quadrupolar echo pulse sequence with full phase cycling and a 
pulse sequence recycle delay of 10s. Due to the low rf power 
transmitter of the spectrometer, 90° pulse width corresponding to 
deuteron nucleus was 20.2 us. The interpulse spacing (x) chosen 
was 60 us. Depending on signal -to-noise ratio, from 24 to 500 
fid's were averaged before doing FT and 32K data memory was used. 
Also, the lOmm sample tube was not spun and the field frequency 
lock not established. 
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The experimental spectra and their simulation shall now 
be discussed for the following two reverse micellar systems : (i) 
PC/Cyclohexane/D 20 (ii) PC/DCC/Cyclohexane/D 20 at room temperature 
(296 K) and at 263 K. 

(a) Spectra at 296 K 

The spectra for a number of R— values for both the 
reverse micellar systems are given on the left of Figs. 
III-4-III.10 and all these spectra show two lines, almost 
identical for the same R value of both micellar systems, namely, 
one broad and one sharp line. The sharp line is not shown in some 
of the experimental spectra as we have expanded the broad peak for 
obtaining an accurate linewidth measurement. 

At first sight at any rate, the appearance of broad and 
sharp lines leads us to the following tentative interpretation : 
the broad line may be due to the water molecules which are bound 
to the polar head groups (see Fig. II. 8) present in the 
interfacial region of reverse micelles and are hence restricted in 
their translational motion ; slow motion of D 2 O molecules does not 
average out the angular part of the quadrupolar Hamiltonian, and 
hence the line is broad. The sharp line may then be rationalized 
on the assiamption that, even at low R values the reverse micelles 
contain some amount of D 2 O as 'free water' either dispersed in the 
organic solvent or present in the inner core of reverse micelles. 
As this latter type of water is very much similar to normal water 
whose Hq is averaged out to zero due to fast motion, the 
corresponding line is narrow. 

If the above interpretation is correct, then one would 
expect that the integrated intensity of the sharp line should 
increase with increasing R value, because free water content 












Fig. III. 9 '^D-NMR spectra (left) and best-fitting simulations 

(right) for PC/DCC/Cydo^sxane/DjO at 296 K; (a) R 

s= 5 and (b) R = 6. 
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should increase as R is increased. Somewhat unexpectedly, 

however, its intensity starts decreasing as R is increased. Also, 
the sharp line is appearing in all spectra around 180 Hz instead 
of around zero Hz (this is due to chemical shift and small 
variations in the magnetic field which had been left 'unlocked'' 
with respect to deuterium frequency when recording our spectra) as 
all spectra have been recorded with respect to normal ^ 2 ® line 
(reference) frequency. This forced us to re-examine the origin of 
this sharp line. When we performed a NMR measurement on (i) 

cyclohexane and (ii) lecithin dissolved in cyclohexane solution 
(i.e., a 'blank' run, which should give no spectrum), a sharp 
line appeared at exactly the same frequency where we obseirve the 
'sharp' component in all our reverse micelles (see Fig. Ill 11 (a) 

and (b) ) . Even the distilled water (H 2 O) in just one FID give a 

2 

strong residual signal in the D NMR (see Fig. III. 11(d)). This 

has been a matter of considerable surprise to us, especially since 

2 —2 

the natural isotopic abundance of D is 1.56 x 10 % and the NMR 

sensitivity at constant field relative to equal number of protons 

-3 

should only be 9 . 65 x 10 

Having thus identified the origin of the sharp line as 
being due to naturally abundant D, we have simulated the broader 
component of all our experimental spectra. The input data used 
for simulations are given in Table III. 3 and Table III. 4, 
respectively, for the two reverse micellar systems. 

The T 2 values (Tables III. 3 and 4) for the 'bound' water 
are seen to be nearly 100 times less than the corresponding 'free 
water' values. This leads one to conclude qualitatively that the 
motion (rotational and translational) of the 'bound' water is 
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slower than that of the 'free' water or, in other words, the t of 
the 'bound' water is longer than that of 'free' water. This also 
substantiates our group's earlier result and the result 
presented in Chapter II, by assuming the Stokes— Einstein eguation 
[II. 6], that the microviscosity of the water near the interface is 
higher when compared with 'free' water microviscosity. 

As expected from our assumed model , the exchange rates 
in Table III. 3 are seen to increase as R is increased, but the 
initial rates of increase are rather slow. There is a substantial 
increase in exchange rate when R is increased above 3. This is in 
line with our earlier findings that the reverse micellar structure 
changes from spherical shape to rod-like or tubular channels . 
In our simulations, the quadrupole coupling constant chosen for 
the free water is the most commonly used value of 220 KHz (48) 
while the 'bound' water quadrupole coupling constant has been 
varied between 150-194 KHz for getting best-fitting simulations in 
each case. Lower coupling constant was also reported earlier by 
Quist and Halle (49) . These workers have reported 150 KHz for 
A0T/D20/2,2,4- trimethyl pentane system. The plausible reasons 
for the observations of lower coupling constants for the 'bound 
water' are : 

(i) change in 0-D bond length due to hydrogen bonding ^ 

(ii) electrostatic interaction between the polar head groups and 
'bound' water molecules, and 

(iii) partial averaging of quadrupole coupling due to the libration 
and single axis rotation (49) as the 'bound' water is loosely 
held by the polar head groups by electrostatic interaction. 

At least in the case of lecithin/ cyclohexane/D20, this would be 
the most reasonable mechanism in view of our earlier 
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findings 


(50,51) 


Table III. 3 


Parameters used in simulating NMR spectra for 
PC/Cyclohexane/D^O 


2 ^ ^ Exchange 
Rate . , 


(10 -"s) (10 ■'s) 4 -1 site 1 


Fractional population 'bound water 
of quadrupole 

site 1 site 2 coupling 

ori^ 0 ^ 2 ^ constant. 


1 

66.3 

44.1 

62.0 

3 

64.6 

60.0 

64.0 

5 

64.3 

96.7 

69.0 

6 

64.3 

72.9 

70.0 

8 

61.8 

52.0 

71.0 

10 

65.2 

23.0 

66.0 

12 

65.2 

21.1 

67.0 


0.99 

0.01 

190 

0.85 

0.15 

175 

0.70 

0.30 

170 

0.63 

0.37 

180 

0.55 

0.45 

165 

0.35 

0.65 

160 

0.15 

0.85 

150 


2 

Table III. 4 Parameters used in simulating D NMR spectra for 
PC/DCC/Cyclohexane/D20 

m 7T^\ T f 2 l ~~ Fractional population 'bound water' 

R -^2^ ^ 2^ ’ Exchange quadrupole 

(10-=s) (10‘^s) STEiT site 2 coupling 

' (10 s ) Dfl^ P(2) constant. 


constant , 
KHz 


1 

67.1 

74.1 

60.0 

3 

67.3 

75.3 

61.0 

5 

62.3 

98.6 

67.0 

6 

61.1 

87.2 

68.0 

8 

60.6 

85.5 

67.5 

12 

61.2 

77.0 

63.0 


0.99 

0.01 

194 

0.85 

0.15 

180 

0.70 

0.30 

170 

0.63 

0.37 

176 

0.55 

0.45 

165 

0.15 

0.85 

155 
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(b) Spectra at 263 K 

Most of the earlier studies of water dynamics in reverse 

micelles have been done at room temperature. In this work, it is 

of interest to investigate water dynamics at subzero temperature 

as the anomalous properties of water are known to occur at subzero 
(52 53) 

temperature ' . It is observed that the water present in 

microemulsion samples does not freeze at the thermodynamic melting 

point of bulk water. The lowering of temperature in this type of 

systems initiates crystallization (ice formation) through the 

mechanism of homogeneous nucleation. The temperature of this 

homogeneous nucleation process decreases as the water droplet size 

decreases. For example, at atmospheric pressure, for droplet size 

of a few hm, the homogeneous nucleation temperature is close to 

-40°C. In the case of our reverse micellar systems, the water 

pool size is about 350A° upto an R value of five^^^^. Very little 

information is available about the supercooling of these systems. 

In aqueous solutions of macromolecules, it has been found that a 

small amount of unfreezable water remains well below the 

(54 55) 

homogeneous nucleation temperature ' ’ . 

Before undertaking the present study, we have not come 
across any report of 'pool water' freezing point in phospholipid 
based reverse micelles. Of possible relevance is a very recent 
report by Wu et al.^^®^ on the freezing of interbilayer free 
water, phosphocholine head group and bound water of fully hydrated 
egg sphingomyelin bilayers. The freezing points of interbilayer 
free water, phosphocholine head group and bound water have 
been reported as — ll^C, — 34^C and around “70 C, respectively. In 
the present study, we have recorded all the spectra at —10 C. 
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since egg sphingomyelin and egg phosphatidylcholine have the same 
phosphocholine head group structure (see Fig. III. 12), we can 
reasonably assume that, at temperature of our study, only the 
'free' water is approaching freezing. 

Interpretation of low temperature deuterium NMR spectra ' 
becomes difficult due to the fact that the phase behaviour of 
these reverse micellar systems at subzero temperature has not been 
characterized to our knowledge. Our interpretation of these 
spectra is, accordingly, rather tentative. 

Spectra of both systems are given in Figs. III. 13 to 
III. 19, and the parameters used for simulation are given in Tables 
III. 5 and III. 6. 

Table III. 5 Parameters used in the 'best-fit' NMR simulation 

for the system PC/Cyclohexane/D 20 


r 

R 

fjd) 

(lo'^s) 

T2(2) 

(lO'^s) 

Exchange 
Rate , ^ 

Fractional population 
of 

site 1 site 2 

p(l) P(2) 

'bound water' 
quadmipole 
coupling 
constant , 
KHz 

1 

14.8 

3.5 

4.23 

0.99 

0.01 

192 

3 

13.8 

1.3 

4.40 

0.85 

0.15 

200 

5 

13.1 

1.1 

4.70 

0.70 

0.30 

200 

6 

12.3 

1.0 

4.90 

0.63 

0.37 

186 

8 

13.0 

0.75 

4.80 

0.55 

0.45 

185 

10 

13.2 

1.0 

7.0 

0.35 

0.65 

182 
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(a) Phosphatidylcholine 


CH 2 CH-CH2-0-P-0-CH2-CH2-N(CH3)3 

0 0 Q- 

1 I 

CO CO 

I I 

(CH2)n(CH2)n (n = 12-24) 

CH3 CH3 


(b) Sphingomyelin ^ 

HO — CH — CH -CH2 — 0-P-0-CH2-CH2-f5(CH3)3 

CH NH 0" 

CH CO 

(CH2),3(CH2)i4 

CH3 CH3 


Fig. III. 12 Comparative structures of Phosphatidylcholine and 
Sphingomyel in 
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^D-NMR spectra (left) and best-fitting simulations 
(right) for PC/Cyclohexane/DjO at 263 K; (a) R = 

1 and (b) R = 3. 


Fig. III. 13 
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Fig. III. 14 ^D-NMR spectra (left) and best-fitting simulations 

(right) for PC/ Cyclohexane/D20 at 263 K; (a) R 

= 5 and (b) R = 6. 
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^D-NMR spectra (left) and best-fitting simulations 
(right) for PC/DCC/ Cyclohexane/D20 at 263 K; (a) 

R * 1 and (b) R = 3. 


Fig. III. 16 
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Table III. 6 Parameters used in the 'best- fit' NMR simulation 
for the system PC/DCC/Cyclohexane/D 20 


R 

Tjd) 

T2(2) 

Exchange 
Rate , - 

(lO^s"-') 

Fractional 

of 

population 

'bound water 
quadrupole 
coupling 
constant , 
KHz 


(10 s) 

(10 s) 

site 1 

P(l) 

site 2 
P(2) 

1 

9.8 

3.9 

4.31 

0.99 

0.01 

155 

3 

9.7 

1.5 

4.52 

0.85 

0.15 

162 

5 

9.5 

0.7 

4.40 

0.70 

0.30 

160 

6 

9.4 

0.7 

4.60 

0.63 

0.37 

174 

8 

9.5 

0.6 

4.50 

0.55 

0.45 

190 

10 

9.6 

1.1 

4.90 

0.35 

0.65 

170 

12 

9.6 

0.7 

4.85 

0.15 

0.85 

180 


From the Tables III. 5 and III. 6, in each micellar system 

the T values for both sites are lower, as expected, compared to 
2 

their corresponding values at 296 K. Exchange rates have also 
decreased nearly 10 ^ times . At this temperature, both the systems 
give similar types of spectra at all R values, except in the range 
of R = 5-8. This is a strong indication for the fact that the 
structure or shape of the DCC-containing reverse micelles at R = 5 
and 6 is aspherical, because the lineshape is a typical powder 
pattern with t) 0. Although the same trend should also be 
observable in the DCC-free system, we have not observed it. This 
may also be due to the difference in the degree of organization of 
water molecules extending into the cores of these two systems. 
Measurements of freezing-points of 'bound water' and 'free water' 
by other specific techniques such as differential scanning 
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calorimetry (DSC) , as well as direct electron microscopic 
measurements of the aggregation size and geometry of these reverse 
micelles, are necessary for guantitating our present observations 
further . 

DMR spectra of the R = 10 and R = 12 cases of DCC-added 
system clearly show the superposition of a sharp line on a broader 
line. It is found from the literature that these lineshapes 
resemble those of lecithin-water sy stem ^ ^ ^ with a D 2 O : dipalmitoyl 
lecithin mole ratio of 10.7 at 43°C and of 14.6 at 45.6°C. Such 
composite lineshapes are characteristic of 'free water' exchanging 
slowly with the 'bound water' . 

The appearance of a small hximp in the middle of an 
. . 2 

otherwise typical powder D pattern expected for an asymmetric 
parameter tj « 0 in the R = 5, 6, 8 cases of our DCC-added system 
and in the R = 8 case of DCC-free system clearly shows that two 
physically distinguishable D 2 O sites exist in each of these 
reverse micelles with different T 2 values. The existence of water 
molecules with distinguishable dynamics in two different binding 
sites within the same thermodynamically and macroscopically 
homogeneous system is an interesting result of our studies. 

Even though bilayer structures are different from those 
of reverse micelles, it is at any rate safe to assume that the 
hydration of phosphocholine head group containing natural 
surfactants (lipids) by water occurs in an analogous manner in 
these two structures. One may then conclude with reasonable 
certainty that, in the reverse micellar systems studied by us at 
263 K, only the 'free water' dispersed as small droplets in the 
cyclohexane and present in the core of the 'water pool' is very 
close to freezing point, in correspondence with results of the 
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layer system reported by Wu et 

In both the reverse micellar systems studied, some of 
le experimental spectra observed at 263 K could not be exactly 
.mulated. This is possibly due to (i) the inclusion of only 
jcular terms of in our computations, and (ii) the use of a 
Lngle Lorentzian lineshape convolution function for the two 
jectral components ; it is expected that the use of two 
Drentizian functions, or a Lorentizian shape function for the 
irrower component and a Gausian one for the broader component, 

Duld improve the simulations further. 

II. 5 Conclusions 

2 

. The appearance of the sharp D line in both reverse micellar 

systems at all R values is unambiguously assigned as being 

2 

due to the naturally aboundant D of cyclohexane and 
lecithin. 

. Computer simulated spectra at 296 K show that T 2 values for 

the 'bound water' are nearly 100 times less than corresponding' . 
'free water' values, providing evidence for the higher 
microviscosity of water near the interface. Substantial increase 
in exchange rates of water in both systems above R = 3 supports 
our earlier findings that the structure of reverse micelle changes 
from spherical to a rod-like tubular shape. 

Computer simulated spectra at 263' K show that the exchange 
rates of water decrease nearly 10^ times compared wi'th the same 
system at 296 K. Resemblance of line shape with typical D power 
patterns for tj “ 0 suggest that the shape of DCC-added reverse 

A 

micelle between R valuesof 5 and 8 may be aspherical. Clear 
indication of exchange of 'free water' with 'bound water' has been 
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provided from the experimental spectra corresponding to R=10 and 
R=12 cases of DCC— added system. The appearance of a small central 
hump for some of the R values studied shows the existence of 
water molecules at two different sites with distinct dynamic 
properties. In comparison with the hydration of similar head 
group containing lipids, only the 'free water' dispersed in the 
cyclohexane and present in the core of the 'water pool' is thought 
to be very close to freezing point at 263 K. 
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CHAPTER IV 

PROTON MULTIPLE QUANTUM RELAXATION OF A DEOXYCHOLATE METHYL 

GROUP IN REVERSE MICELLE 

IV. 1. Introduction 

It is now well established that the structure and 
dynamics of macromolecules and molecular aggregates can be 

obtained by extracting cross-relaxation times of coupled nuclear 

• ( 1 2 ) 
spins through two dimensional NMR methods' ' . However, 

interpretation of the resulting data from the molecular 

aggregates is rendered difficult due to the modulating effects of 

motions of a group in a molecule, or of the molecule itself, or 

even of the aggregate as a whole, these various motions often 

occurring at different time scales. For example, the methyl group 

of a surfactant in a reverse micelle may be complicated by the 

following motions : 

(i) motion of a methyl group about its axis, 

(ii) motion of the whole molecule in which the methyl group is 
present and 

(iii) motion of the reverse micelle as a whole. 

Additionally, for systems such as reverse micelles consisting of 

heterogeneous molecules, the inter- and intra- molecular 

interactions can also affect interactions among spins thereby 

adding to the complexity of our interpretation. 

These difficulties were overcome by Kay and 

Prestegard^^^ , Wong et al.^^^ and Wang and Hwang by using the 

results of an one-dimensional analog of n-quantum filtered 2D 

(2 4^ 

experiment (n=0,2,3) combined with transverse' ' ' and 

. . 

longitudinal' ' relaxations. 
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In an effort to study the dynamics of polar head groups 
of surfactants in a reverse micelle and to get dynamic information 
at the interface as well as the interaction of reverse micellar 
'pool water' with the polar head groups, we first attempted to 
study the -N (CH 2)3 of lecithin theoretically by calculating the 
response of this proton-spin system to rf pulses. However, the 
mathematical complexity of handling this Ag case was soon realised 
and we then took up the well studied sodium deoxycholate as the 
surfactant for forming reverse micelles. 

This chapter contains (i) a brief introduction to 

cross-relaxation effects in coupled nuclear spin systems and the 

consequent appearance of forbidden peaks in multiple quantum 

filtered NMR, (ii) our theoretical calculations for an A- spin 

+ ^ 

system (formed by the nine protons of -NCCH^)^), and (iii) 
experimental results of zero- and double- quantum filtered ID 
analog of 2D experiments on deoxycholate reverse micelles. 
Motional parameters have been extracted by way of matching the 
experimental data profile to the theoretical one assuming an 
appropriate model for the reverse micelle, and certain overall 
conclusions have been made based on our results. 

IV. 2. Cross-relaxation Times 

The fast molecular motions of macromolecules and 
molecular aggregates can be probed by C spin— lattice relaxation 
time (Tj^) and nuclear Overhauser effect (nOe) measurements 
performed at different magnetic field strengths. The spectral 
densities of fast molecular motions fall in the frequency range of 
10® - 10® s”^. But the information about slow motions or 
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procGssss such 3s diffusion, cheitiical ©xchangs etc. (correlation 
time (T^) > 10 s) which have negligible spectral density^^^ near 
or above the Larmor frequency are obtained by spin echo 
measurement of spin-spin relaxation time (T 2 ) . The rates for 
these processes are obtained from the pulse spacing dependent 
dephasing rates in 2-pulse or multiple-pulse experiments. 
However, it was believed earlier that the relaxation rate 
dependence on pulse spacing is due to the existence of low 
frequency components of molecular motions, whereas on the contrary 
it is due to the existence of cross-relaxation between components 
of transverse magnetization of spin coupled multiplets. This is 
explained below. 

Two important relaxation mechanisms for spin-1/2 nuclei 
in liquids in the absence of inhomogeneous magnetic field and 
chemical exchange are : inter- and intramolecular dipole-dipole 
and spin-rotation interactions. Much of the earlier NMR 
literature has dealt with intramolecular dipole-dipole interaction 
of two coupled spins only. For the case of three or more coupled 
spins, assuming that the motion of any two nuclei is \incorrelated 
with that of a third (ie. , the case where cross-correlation 
functions go to zero) , it can be readily shown that the total 
relaxation rate is just the sum of all individual pairwise 
dipole-dipole interactions. In reality, however, for a methyl 
group the motion of the third proton is clearly determined by the 
motions of the other two. An elegant pictorial explanation of 
this has been provided by Werbelow and Grant ^ ^ , and this figure 
is reproduced in Fig. IV. 1 
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(a) SINGLE PAIR OF TWO IDENTICAL SPINS 


SPINS 


UNIQUE TIME CONSTANT ^ 

'communication via dipolar 

COUPLING 


LATTICE 

.NONSPIN DEGREES OF 
•FREEDOM ‘ 




(b) CHAIN OF COMMUNICATION IN A SYSTEM OF THREE IDENTICAL SPINS 



Independent 

pairwise 

interactions 

only 

(incomplete 

description) 


Correlated pairwise 
interactions included 
(complete description) 


Longitudinal relaxation process in spin systems 
through dipolar interaction (adapted from ref. (7)) 

(a) for a single pair of identical spin -1/2 (I and J) 
nuclei ; 

(b) for a system consisting of three identical spins 
(I, J and K) . 


Fig. IV. 1 
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For a system of two identical coupled spin -1/2 nuclei, 
I and J (see Fig. IV. 1(a)), the longitudinal relaxation 
processes are characterized completely by one unique time 
constant. When we consider spin systems consisting of three 
identical spins and assume that the system has only three 
independent pairwise interactions, we may, again, always propose 
that a unique time constant exists for this system. Strictly, 
however, it is realised from Fig. IV 1(b) that these three 
pairwise interactions are interlinked and any resultant change at 
any one nucleus due to a coupling with the other two must be 
correlated. Therefore, the correct relaxation equation should 
contain terms corresponding to these correlation processes. The 
correlation processes are of two types : 

(i) auto-correlation, when the orientation of I-J dipole-dipole 
vector at time t is related to the orientation of the same 
vector at time t+r; 

(ii) cross-correlation, when the orientation of I-J dipole-dipole 
vector at time t is related to the orientation of J— K dipolar 
vector at time t+r. 

Void and coworkers and Grant and coworkers have 

calculated these cross-correlation spectral density terms by 
fitting T. and T_ relaxation data; they have shown the usefulness 
of these terms for studying the motional properties of spin 
systems. Later on two other methods for extracting these terms 
have been reported, namely, the use of polarization transfer pulse 
sequence (DEPT and INEPT) and of the heteronuclear multi-pulse 
schemes . 
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The spin systems consideired till to— date are mostly AX 
• 13 19 1 

or having C or F and H. However, the above methods 

invariably require isotopic enrichment, which is difficult in the 
case of biomolecules. To overcome this and other problems related 
to the extraction of cross-correlation terms from T^^ and T 2 data, 
Kay and Prestegard' have proposed an interesting ID analogue of 
double quantum 2D experiment. Their method is explained in detail 
in the next section. 

IV. 3 Extraction of Cross-correlation Spectral Density Terms 

Using Multiple Quantum NMR 

In 1985, Muller et al.^^^^ reported the appearance of 

forbidden cross-peaks, in violation of coherence- transfer 

(12 13 ^ 

selection rules' ' ‘ derived neglecting transverse relaxation, 

in proton multiple quantum filtered correlation spectroscopy and 
multiple quantum NMR. The selection rules governing multiple 
quantum filtered 2D correlated spectroscopy of weakly coupled spin 
systems are: 

(i) The 'diagonal' peaks, representing inequivalent nuclear 
sites, in n-quantum filtered spectra imply that the spin in 
question has resolved scalar couplings to at least (n-1) 
equivalent or nonequivalent spins. 

(ii) The 'cross' peak, representing J-coupling between two spins 
in an n-quantum filtered spectrum indicates that, in addition 
to their direct coupling, there are nonvanishing couplings to 
a common set of at least (n-1) additional spins. 

For a macromolecule having a methyl groi^ attached to a 
tertiary carbon, for instance, these rules predict that a single 
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resonanc© line corresponding to this methyl should not appear in 

twO”Or higher guantum filtered correlation spectroscopy. However ^ 

contrary to this prediction, MQ peaks do arise in similar 

degenerate spin systems ' due to multiexponential (or cross- 

relaxation between pairs of spins) ' transverse relaxation. 

Although the appearance of such forbidden cross-peaks can lead to 

wrong resonance assignments in 2D NMR, these peaks can 

nevertheless give useful motional as well as structural parameters 

( 3 ) 

for a macromolecule' ’ . A brief summary of equations derived for 
A^ spin system for zero- and double quantum filtered ID analogue 
of 2D NMR experiments by Kay and Prestegard ^ ^ ^ and Wong et al. 
is given below, since we have used their derivations as such for 
our work on deoxycholate reverse micelles. 

The energy level diagram and basis functions for an A^ 
spin systems, typified by the methyl (“CH^) group, are given in 
Fig. IV. 2. The pulse sequence and the schemes used for 
transmitter and receiver phase cycling' ' ' in order to observe 
the zero- and double quantum filtered forbidden peaks selectively, 
are given in Fig. IV. 3 and Tables IV. 1 and 2, respectively. 

The first two Tr/2 pulses in the sequence (Fig. IV. 3) 
generate different orders of multiple quantum coherence, and the n 
pulse refocuses the precessional loss due to chemical shift 
offset during the period 2 t. Therefore, only transverse 

relaxation effects are present at the end of 2 t in the case of 
unresolved scalar coupled spin systems. It is also assvimed that 
the relaxation effects originate from proton dipolar interaction 
only^^^^ . The spin-lattice interaction Hamiltonian of this system 
is written as 
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Table IV. 1 Phase cycling scheme used in the selection of zero 

quantum filtered spectra (adapted from ref. (4)). 



U CM 

■e- 

receiver phase 

^2 

receiver phase 

X 

X 

X 

•X 

-X 

-X 

y 

X 

X 

-y 

-X 

-X 

-X 

X 

X 

X 

-X 

-X 

-y 

X 

X 

y 

-X 

-X 

y 

Y 

y 

-y 

-y 

-y 

-X 

y 

y 

X 

-y 

-y 

-y 

y 

y 

y 

-y 

-y 

X 

y 

y 

-X 

-y 

-y 


Table 

IV. 2 

Phase cycling 

scheme 

used in 

the selection of 



double (juantum 
(3)). 

filtered 

spectra 

(adapted from ref. 


^2 

receiver phase 

h 

receiver phase 

X 

X 

X 

X 

-X 

-X 

y 

X 

-X 

y 

-X 

X 

-X 

X 

X 

-X 

-X 

-X 

-y 

X 

-X 

-y 

-X 

X 

X 

y 

-y 

X 

-y 

y 

y 

y 

y 

y 

-y 

-y 

-X 

y 

-y 

“X 

-y 

y 

-y 

y 

y 

-y 

-y 

-y 


^phase of excitation pulses. ^Phase of observation pulses. 
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?{(t) = I e 

i<j 


ID 


I (- 1 ) 

m 


m 






+ E 


I (-1) 
i m 


"" T^(i)B^®(i,t) 


[IV. 1] 


In the above equation. 


= I (671) 


|(67r)'/y5] y^ 




Tj(i,j) 

ij) = 


nin 


ID 


B"”‘(i,t) = 


the dipolar spin tensor operator of rank 2 
connecting spins i and j . 

the spherical harmonic of rank 2 whose value depnds 
on the orientation of the i-j dipolar vector 
relative to the external field, 
spin tensor of rank 1, 

magnetogyric ratio of the methyl proton i, 
distance between protons i and j in the methyl 
group , and 

the lattice random field interaction operator. 


The first term in eqn. [IV. 1] is due to the intramolecular dipolar 

interactions among the three protons of the methyl group and the 

second term takes into account all other interactions. The second 

. ( 7 ) 

term is approximated to be isotropic' , i.e., 

|B°(l,t)|2 . |B+\i,t)|2 _ |B-'(i,t)|^ 

The relaxation processes for the methyl group are described, by 
using Redfield theory , in terms of evolution of the density 

matrix (p) according to the equation of motion given below . 
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5 = T}q 

dt 

where R is the Redfield relaxation matrix and 


P = 


14 

*26 

^35 

’47 

’78 


[IV. 2] 


[IV. 3] 


The subscripted numbers represent the spin energy levels 

shown in Fig. IV. 1. The dipolar part of the relaxation matrix is 
^7) 

given by ' ' 

HI, n 

<clT5(k,l)l<r> - <e |T”(i, j) |P> 

P 

<e|T5(k,l)|<r> - r J?5kl<“pg> <“|T^(i-J)|P> 

</3 |t 5(1 c, 1 ) |v>| (IV. 4) 


where (ca) is a typical spectral density function defined by 

X J JvX 




^ij ^kl 


<Y^(^ij(t)) Y^(^j^i(t+T))> 


cos OTdT 


(IV. 5) 
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The spectral density functions given above are calculated under 
two situations, namely, spin being (i) autocorrelated and (ii) 
cross correlated. These terms have been explained already in 
section IV. 2. A complete listing of all the elements of R in 
terms of auto- and cross-spectral densities can be found in ref. 

( 3 ) . The spectral density function corresponding to the second 

. . <’41 

term of eqn. [IV.l] is given by^ ' 


n 


j (nu) = (-1) r 


11 


<B^^(i,0) 


Br^(j,t)> cos(n£Jt)dt, 


[IV. 6] 


For the specific case of a methyl group in a macromolecule, the 
above spectral densities are evaluated by using an appropriately 
modified "cone model” of Brainard and Szabo^^®^ . Fig. IV. 4 thus 
describes the dynamics of the methyl groups in a spherical 
micelle. Using the motional parameters defined in Fig. IV. 4, eqn. 
[IV. 5] is written as 




= 6 


(-l)^3yV 


lOr® 


1 c2 6D 

^ (6D)^+w^ 


1 

4 


3 

4 


3 

4 


(1-S^) [6D+6Dy (1-S^) ] 
[6D+6Dy (1-S^)]^+W^ 


COS(nij,^l) 

cos(nijj^i) 


S^(6D+4D.) 

2 2 ^ 

{6D+4D^) +w 

(1-S^) [6D+4D^+2Dy (1-S^) ] 
[6D+4D +2D /(l”S^)]^+w^ 


[IV. 7] 


In our model (depicted in Fig. IV. 4), the overall reorientational 
motion of the micelle is characterized by a correlation time - 
the internal rotation of the methyl group about its ^ 
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has a correlation time , and the correlation time 

corresponding to the wobbling motion of the methyl axis is 

^ methyl group is assumed to be 

w 

wobbling within a cone of angle 6 with a frequency of 6D^ and an 
order parameter of S. S is related to e as^^®^ 


^cone = (1+cose) . 


[IV. 8] 


In equation [IV. 7], ^ijkl angle between i-j and k-1 ^H-^H 

vectors. The value of this angle for auto-correlated relaxation 
process is whereas it will be = 60°, 

for cross-correlated relaxation ij kl, as the three protons of a 

methyl group occupy the vertices of an equilateral triangle. 

(7) 

Using standard techniques of matrix algebra' one can 
solve eqn. [IV. 2] to get the relevant density matrix elements for 
the zero-and double quant\im coherences after the second (90) ^ 
pulse. These elements are 
(i) for zero-quantum coherence : 


Pll 

Pll 

P44 

P44 



II 

1 

Pi4(2t“) + 

+ P7j(2x‘)j 

(2X^) 




1 

4 

['^Pi4(2t“) 

- P„(2t') + P^st^i;')] 

(2t'^) 

= - Pyji2x^) 



V [IV. 9] 
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(ii) for double- quantum coherence : 


Pi7(2'c'^) = - P47(2r") + p^g(2T~)j 

[p14<2'c") - ^ P47(2t') + p^g(2T“)j 


[IV. 10] 


In the above equations [IV. 9] and [IV. 10], and in the following 
ones, the + and - superscript signs represent the density matrix 
just after and before, respectively. The multiple quantum 
coherences due to nonuniform relaxation of the transverse 
elements P^jr P 78 ^3 system, created by the 

second rr/2 pulse, evolve during the time A. These nonobservable 
coherences are converted into detectable single-quant\am coherence 
by applying the last (n/2) pulse. The zero- and doxible-quantum 
coherences are selectively observed by employing the 
respective phase-cycling schemes given in Tables IV. 1 and IV. 2. 
The elements (P 2 ^ 4 r ^ 47 ' ^78^ density matrix for 

zero-quantum selective excitation are : 

P^4(2x+A'^) = [p^^(2t+A") + P44(2t+A") + p^^(2x+A^) 

+ Pgg(2T+A”) ] 

P4^(2t+A‘^) = I |^3p^^(2T+A") - p^^(2x+A~) + p^j^(2x+L ) 

-3P44(2t+A") j 

P7g(2T+A'^) = j^p^^(2T+A") + P44(2 t+A") + p^^{2x+A ) 

+ P 88 ( 2 ^+A”) ] 
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excitation are j-ecrive 


P^4(2t+A ) - 1 Lp ^2t+a“) -3p (2x+tr\ 4.-^ 

L 1/ ) +3p^^(2T+A ) 

+ P84 (2t+A~) j 

= !| [-P„(2..A-, . p^3(2.+a-, - p,^(2..A-, 

+ P84(2T+A") J 

P,3(2«A^ = 1 [3p^,(2r«-, + p^3(2„a-, - p^^fj^A*, 

-3Pg4(2T+A") J J 


[IV. 12] 


The 

period t > 2 t 


observed magne t ization 
+ A is given by 




"signal” during the 


\= rh Tr(pi^) 

“ Pl4 2^47 + P,» + P 


14 


+ 2p 


47 Pys) 

[IV. 13] 


Kay and Prestegard^^^ and Wong et al.^^) have shown that, in 

essence, one can calculate theoretically the MQ filtered forbidden 

intensity profile as a function of MQ evolution time by 

using eqns. iv.ii or IV. 12 and IV. 13. They have also demonstrated 

these profiles can be matched with those of the experimental 

uta, whereby useful motional parameters for the spin system 
could be derived. 
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IV. 4. The Ajj Spin System (e.g ; the -N (CHg)^ head group): 

Derivation of Pure State Basis Functions and 
Demonstration of Difficulties in Mathematical 
Tractability 

A system of nine coupled 'like-spin 1/2' nuclei will 

9 

have a total number of 2 = 512 eigenstates or basis functions. A 

convenient spin basis is necessary to formulate these many states 
for further calculations. The most convenient basis in which to 
represent these many states is the one adopted by Bersohn and 
Gutowsky^^®^ / Lalowicz et al.^^^^ and Raghunathan and Ray^^^^ 
based on coupling of angular momenta , for coupled spin 
systems. The representation is given in |IM> basis where 

A = If + I 2 
3=13+14 

D = I 7 + Ig 
E = A + B 
F = C + D 
G = E + F 
I = G + I 9 

and M is the z-component of 1 , the total angular momentum quantum 
number, projected along the direction of B^ such that -I < M < + 
I. The abovementioned coupling pattern is given below in a 


simpler form 
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^1 -"2 

La— I 

I 


^3 L 
LrJ 


^5 ^6 

LcJ 

I 


I? ig 

LiJ 


The complete basis set for representing the case will therefore 
become |ABCDEFGIM>. Under the operation of the Zeeman 
Hamiltonian, these 512 pure states can be arranged according to 
their energy. The energy levels, and degeneracy (no. of states) , 
are given in Fig. IV. 5. 

While considering these many pure states, and giving a 
true representation to each state by taking a linear combination 
among these and going further to give a matrix representation for 
the whole system for calculating its evolution under an rf 
excitation, one soon realises that the problem is not easily 
tractable mathematically. Further theoretical work on the Ag 
system was therefore abandoned, and an easier micellar model was 
pursued instead, choosing the well-studied sodium deoxycholate as 
the surfactant. The structure of the deoxycholate anion is given 
in Fig. IV. 6. we recall here that sodium deoxycholate is the 
well-known 'bile-salt' which is useful for forming large and 
well-defined micellar structure models. 


IV. 5 Studies on the Deoxycholate (DOC) Reverse Micelles 

IV. 5.1 Experimental Details 

The reverse micellar system chosen for our study is 

sodium deoxycholate/benzyl alcohol/H 20 , as simila yp 
soap/alcohol/water systems form reverse micellar so u ions 
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Fig. IV. 5 


Degeneracy present at each M value ofthe Ag spin 
system, under Zeeman Hamiltonian only. 
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Sodium deoxycholate used was procured from S.D. Fine-chem, India, 
and was further purified by dissolving it in ethanol/acetone 
mixture and recrystallizing the salt after filtering the hot 
solution through a charcoal bed. Benzyl alcohol (Fisher) was 
distilled once before use, and the water used was doxible-distilled 
and deionised. 

Proton zero-and double-quantum filtered ID analog's of 
the 2D experiment were done at 297 K on a Bruker WM400 FT NMR 
spectrometer operating at the proton frequency of 400 MHz. The 
pulse sequence and phase cycling schemes used are given in Fig. 
IV. 3 and Tables IV. 1 and IV. 2, respectively. The 90° pulse width 
corresponding to protons in this spectrometer was 16 jis . All 
spectra were recorded with the following conditions : 192 scans, 
sweep width of 3496 Hz, relaxation delay of 5T^ and 8 K data 
points. All spectra were zero-filled to 16 K before Fourier 
transforming, and recorded in magnitude corrected mode, 

VRe(w)^+Im(a))^, where Re(w) and Im(w) are data from real and 
imaginairy channels, respectively. The data were not apodized. 
The multiple quantum evolution time (A) used was 5 ms. 

IV. 5. 2 Results and Discussion 

In the structure of DOC (Fig. IV. 6 (a)), there are three 
methyl groups, labelled as 18, 19, and 21. The proton ID NMR 

spectrum of our DOC reverse micelles is given in Fig. IV. 6(b). The 
peak selected for multiple quantum NMR study is that of methyl 
group-18, as it is well-resolved and can be considered as an 
isolated A^ spin system ^ . Theoretical data were generated by 
using a program developed by Hwang (Dept, of Chemistry, National 
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ppm 


la IV 6 (a) Structure of the Deoxycholate(DOC) anion 

Lg.IV.6 (a) reverse micelles 

Se^^er'ins ot nethyl groups given m (a, . 
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Taiwan University, and Institute of Atomic and Molecular Science, 
Academia Sinica, Taipei, Taiwan, Republic of China) . 

We note that our model for the reverse micelle is simply 
the reverse of the one given by Kay and Prestegard^^^ (for the 
micelle), and is illustrated in Fig. IV. 7 . The symbols used in 
this figure are defined in a manner similar to the ones given for 
Fig. IV. 4. The input data used for getting double-quantum 
filtered forbidden peak intensity profile theoretically for R = 
2.88 and R = 11.52 are given in Table iv.3. The theoretically 
calculated as well as experimentally determined DQ-profiles 
corresponding to both R values are given in Figs. IV. 8 and IV. 9. 


Table IV.3 Input parameters used for theoretical DQ profile 

generation 


R 

S 

T^(10"^S) 

T^(10"^^S) 


j ( 0 ) 

j (W) 

2.88 

0.87 

3.1 

4.6 

46.0 

15.4 

-0.75 

11.52 

0.78 

3.7 

3.9 

31.5 

8.6 

-0.55 


Table IV.3 shows that the isotropoic rotational 

correlation time (r ) of the reverse micelle increases as R 

o 

increases, which indicates that, as expected, the size of the 
aggregate increases at higher R values thereby showing increase in 
However, decreases as R increases. This may be due to a 
less— compact structure at R = 11.52 as compared to that at R — 
2.88 and/or due to an increase in the overall micellar size which 
generates more space for methyl-18 'wobbling^ motion. The order 
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'Inverted cone' model describing motion of 
group-18 present in DOC molecule in our 
micelles. 


methyl 

reverse 


Fig. IV. 7 






2t(10'^s) 


Fig. IV. 8 400 MHz proton double quantum peak, intensities as a 

function of evolution time for methyl-18, 

recorded for the R =2.88 DOC reverse micelles; 
the solid line represents the profile generated 
theoretically and o's are experimental points. 
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2T (W’s) 


400 MHz proton double quantum peak intensities as a 
function of evolution time for methyl-18, 
recorded for the R =11.52 DOC reverse micelles; 
the solid line represents the profile generated 
theoretically and o's are experimental points. 


Fig. IV. 9 
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parameter (i.e., ordering of methyl group-18) shows a marginal 

decrease with increase in R values indicating that as aggregate 

size increases the relative ordering of methyl group-18 is less, 

and decreases as well . The theoretical profile generation is 

very sensitive to the input variables S, r^, j (o) , and j (w) . An 

almost similar theoretical profile can be obtained by using more 

than one data set^^'^^, as there are six input variables for its 

generation. The parameters given in Table IV. 3 is one such input 

variable set. Realistically, these parameters may have to be 

optimized further either by giving more limiting conditions or 

13 

from results of other independent (e.g. light scattering , C T^, 
nOe at several magnetic fields and using the recently reported 
differential line broadening (DLB) effects in degenerate spin 
systems) experiments at each R value. 

Experimental zero-quantum filtered forbidden peak 
intensity profiles as a function of the excitation time are given 
for R = 2.88 and R = 11.52 in Figs. IV. 10 and IV. 11, respectively. 
Unfortunately, the data appear very scattered for R = 11 *52 

system. Due to this, and due to the nonavailability of data 
for both R values, we have fitted the experimental zero quantum 
peak intensities to a simple single exponential function only and 
we have not undertaken any theoretical intensity profile 
generation for getting other meaningful information. The solid 
line (exponential fit-line) in Figs. IV. 10 and IV. 11 shows clearly 
the need for multiexponential function fit and also indicates the 
existence of cross— relaxation among spins. 
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Fig. IV. 10 400 MHz proton zero quantiom peak intensities as a 

function of evolution time for inethyl-18, recorded 
for the R *2.88 DOC reverse micelles; the solid 
line corresponds to the single exponential fit of 
experimental data (o's) . 
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400 MHz proton zero quantum peak 
function of evolution time for 

for the R =11.52 DOC reverse micelles; toe solid 
line corresponds to toe single exponential fit 
experimental data (o's) . 



Fig. IV. 11 
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IV. 6 Conclusions 

Using recently reported techniques, multiple quantum 
filtered ID analog of 2D experiments has been performed on DOC/ 
benzyl alcohol / H 2 O reverse micelles. For the DOC reverse 
micelles, the input data list used to generate a theoritical 
double quantum peak intensity profile which provides the 
'best-fit' to the experimental profile leads us to infer that the 
size and motional parameters offer good support to a general, 
intuitive model of spherical-structure expansion as a function of 
water addition. The nonuniqueness of the input data required to 
generate a particular profile may be eliminated by using 
information/data from additional NMR experiments. Although the 
input data list is in agreement with that of micelles reported 
earlier, we cannot comment further about it because its uniqueness 
has not been verified by us. As expected, the zero quantum peak 
intensities as a function of evolution time does show 
multiexponential functional form. 

We have also shown the mathematical difficulty in 

+ 

treating A^ spin systems (e.g. ,-N(CH 2 ) ^ overcomes this, 
an easy access would be available for probing the polar head group 
of lecithin-based reverse micelles. 

To a reasonable extent, we can conclude that meaningful 
motional and structural parameters can be obtained for a methyl 
group present on the surfactant molecule and aggregate as a whole 
by extracting the cross-relaxation effects through proton multiple 
quantum NMR and analysing them in combination with conventional 
NMR measurements. 
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CONSTRUCTION OF A VERSATILE MAGNETIC FIELD GRADIENT CONTROLLER AND 
STUDY OF SELF-DIFFUSION OF MICELLAR WATER 


V.i Introduction 

A magnetic field gradient may be described. 


by a second 

rank tensor 

. 

It 

is 

written 

in 

coordinates 

as a dyadic , 
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in general, 
Cartesian 


[V.I] 


However, for a small gradient imposed on a high static 

magnetic field (B^) , i.e., » |AB|, where B^ defines the 

z-axis of the system, the components AB i and AB j can be 

X y 

discarded. The field gradient is then equivalent to retaining 
only the last row of 'S in eqn. [V.I], whereby the vector G 
results , with components 


G 


X 



[V.2a] 



[V.2b] 


G 


z 



[V.2c] 


As shown in Fig. V.I, three types of magnetic field gradient have 


been employed in NMR, namely. 
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(i) 

(ii) 

(iii) 


J.O^ 


linear, 

switched or pulsed linear, 
oscillating or time-varying. 

The choice of a particular type of gradient is often 
dictated by the type of NMR parameter to be measured. For 
example, diffusion and flow measurements would require one of 
the first two types of gradient, whereas NMR 'imaging^ variously 
employs all three types. 

In the next (second) section of this chapter, we review 
the methods available in the literature for generating the three 
types of gradients. In section three, detailed discussion is 
made of a particular circuit design for generating all the three 
types of gradients, including our own modifications for 
constructing a versatile magnetic field gradient controller. 
Performance test results of this novel circuit, and measurement 
of self-diffusion coefficient of H 2 O in Lecithin-based reverse 
micelles, are given in the subsequent two sections. The last 
section incorporates the overall conclusions of this chapter. 

V.2 Types of Magnetic Field Gradient Circuits : 

A Survey 

A literature survey of how magnetic field gradients have 
been generated and used in NMR applications is given in the 
following three subsections. 

V,2.1 Linear Magnetic Field Gradients ; 

The idea of using linear* magnetic field gradients in 
NMR is almost as old as NMR itself . The early workers^ 

It would perhaps be more correct to say 'gradients constant in 
space' 
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realized that gradients or inhomogeneities in were often the 
limitation in obtaining narrow NMR absorption lines in liquids. 

A linear gradient is employed in self diffusion 
coefficient measurements and more recently, in the important area 
of NMR imaging. Since linear magnetic field gradient along any 
direction can be easily generated with a pair of current-carrying 
coils, numerous methods are available in the literature. As a 
brief sampling, a number of methods published within the last 
decade may be cited. 

V.2.2 Switched or Pulsed Linear Magnetic Field Gradients 

The homogeneous (uniform) magnetic field pulsing 
techniques were pioneered by Kapitza^^^^ who produced a 50 T 
field by a battery discharge through a small low resistance coil. 

Since Kapitza's work, a number of workers have achieved pulsed 

. . . . fl3) 

high and very high magnetic fields' ' . 

In 1955, Anderson et al. demonstrated how nuclear 

spin echoes in the presence of pulsed field gradient may be used 
for the purpose of information storage. Although they noted the 
effect of diffusion on their experiment, they did not further 
develop the idea for the measurement of self-diffusion 
coefficients. McCall et al. realised the experimental 

potential of this technique for measuring self-dif fusion 
coefficients, but did not undertake a mathematical analysis or 

experimental test of the method. 

Following McCall et al.^^^^, Stejskal and Tanner 
presented the methodology, first experiments with the pulsed 
gradient and a detailed analysis Tanner's*^” pulsed field 
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gradient circuit has two principal parts (i) the coils for 
producing linear magnetic field gradient (ii) transistor 
switches for gating the current through these coils. The 
transistor switch consisted of a 2N1480 and a 2N174 connected 
together in a Darlington (series) mode, with some minor 
modifications of the circuit elements to improve switching 
characteristics ^ ^ . This transistor switch was later further 
refined by Tanner himself . 

Many circuit designs have been published since. 
Different schemes for producing pulsed linear field gradient have 
been discussed in detail by Gross et al.^^®^, Karlicek and 
Lowe^^^^, Hrovat et al.^^®^, and Callaghan et al.^^^^. 

With the advent and rapid progress of NMR imaging, the 

need for realising fast switched or pulsed linear magnetic field 

gradient has increased further. In 1975, Kumar et al. have 

produced pulsed magnetic fields by means of a fast solid state 

. ( 23 ) 

DIP relays for use in NMR Zeugmatography. Hutchison et al.' 
provide perhaps the most detailed design for generating pulsed 
magnetic fields by using SCR's to charge capacitors which are 
then discharged into gradient coils. More recent technical 

advances for generating pulsed linear magnetic field gradient 
have been published^^'^^^^^ . 

V.2.3 Oscillating or Time-varying Gradients 

Oscillating or time— varying gradients are presently 
used extensively only in NMR imaging. Since the early experiments 
of NMR imaging using gradient fields for spatial 
selectivity(25.26) ^^^e variety of novel imaging approaches has 



165 


been used. 

Imaging experiment with oscillating gradients was first 

proposed by Hinshaw ' , who applied sinusoidally oscillating 

field in X, Y and Z directions in such a manner that only one 

point “ the sensitive point “ of the sample experienced a time — 

independent field. This sensitive point is the only element 

whose NMR signal is not modulated by the alternating gradients, 

and therefore it alone passes through a low-pass filter. This 

technique is simple and direct, and requires no calculations or 

computer for producing the image; it did in fact produce the 
. ( 29 ) 

first 3D NMR images' The mathematical solution of the NMR 

signal produced by a series of pulses in the presence of 
time-dependent gradient has also been published^^*^^ . 

Since the oscillating gradient method requires a longer 
signal acquisition time, it is not competitive with other imaging 
modalities. In addition, the impulse response or localization 
function is relatively inflexible, and results in poor 
localization . 

Time-dependent gradients have subsequently been 

employed by Mansfield and Mansfield and Pykett^^^^ in their 

"Echo-Planar Imaging (EPI)" method. Following these workers, a 

(34) 

modified EPI method has been proposed by Ljunggren . In the 
original experiment, only an alternating gradient with a 
rectangular waveform was applied, but the use of cosinusoidal 
gradient waveforms have also been theoretically analysed by Feiner 
and Locher^^^^ and by Tropper^^^^. Experimental results in the 
presence of sinusoidal or cosinusoidal time-varying gradients have 
been given by Frahm and Hanicke^^^^ . 
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Following up on the proposal of Tropper^^®^, Macovski 

/ 3 n 3 8 ) 

et al . ' have introduced newer techniques in MRI using 

time-varying gradients, • the basic innovative feature consisting of 
post processing the FID recorded in the presence of time-varying 
gradients. It is shown by Macovski et al. that (i) it is 
possible to achieve very good (''narrow'') localization functions 
for the imaging process and the processed images have a 

significant insensitivity to inhomogeneity and chemical shift^^^^ 
and (ii) images are possible even in the presence of B 

o 

inhomogeneities, and one can simultaneously acquire data from more 
than one dimension . 

Norton proposed a new imaging scheme employing a 

time-varying gradient that simultaneously oscillates and rotates. 
This method offers the potential of great imaging speed, with the 
advantages of reducing or eliminating switching transistors, 
reducing the system bandwidth, enhancing frequency selectivity and 
reducing power dissipation by utilizing resonant circuits. 

Finally, oscillating gradient methods, compared to methods 

employing stationary gradients, are found to be inherently less 

sensitive to field inhomogeneities when yAB < where AB is the 
maximum field inhomogeneity and is the modulation frequency. 

Further novel approaches to magnetic resonance imaging 

( 41 ) 

using time-varying gradients have been published recently 

V.3 Construction of a Versatile Magnetic Field Gradient 

Controller 

In our laboratory, the need for a versatile magnetic 
field gradient controller arose not only for studying/measuring 
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the translational self-diffusion coefficient of 'pool water^ 
protons in reverse micelles but also because of our continuing 
interest in NMR imaging studies of biosystems ^ ^ . 

The first circuit for the generation of linear and 

oscillating gradients, published in 1979 by Lauterbur et 

was based on the principle that linear field gradients 

corresponded to vectors, and could therefore be generated from 

their orthogonal components. Their microprocessor based 

controller was used in two-dimensional reconstruction 

experiments and for sensitive point scanning ^ ^ * The same 

idea was later adapted to construct a device for three dimensional 

(45) 

imaging purposes' * . 

In 1982, Fitzsimmons described a simple design for 
generating an oscillatory gradient and suggested methods for 
generating linear and pulsed linear magnetic field gradients from 
the same circuit with a little modification/addition of 
components . 

Based on the work of Fitzsimmons ^ ^ , we have designed 
and constructed a field gradient controller for generating all 
the three types of gradients by simply making modifications in 
order to increase the resolution of the oscillatory gradient. In 
our circuit, certain complexities of the original design^ ^ have 
been simplified. 

The block diagram of the complete gradient control 
system used in this work is given in Fig. V. 2. 
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Block diagram of 
sys'tem used in this work 


complete gradient control 


Fig.V.2 
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The hardware may be divided into five subsections : 

(i) digital oscillator and DC reference power supply, 

(ii) computer interface and manual switch for gradient 
controller, 

(iii) circuit for generating complementary outputs, 

(iv) current driver circuit, and 

(v) display circuit. 

The design logic and detailed circuitry for generating 
all the three types of gradients are described in detail below. 

V.3.1 Design Logic 

The generation of linear and oscillating field gradients 
is explained pictorially in Fig. V.3. The different 

configurations of Fig. V.3 lead to the simple idea that with a 
pair of field gradient coils, and by controlling the relative 
amplitude of the direct/alternating current for each coil, one 
can generate easily the linear/oscillating magnetic field 
gradient. Following this basic idea, we have designed and 
constructed a versatile gradient controller. The necessary 
circuits are described one by one below. 

V.3. 2 Digital Oscillator and DC Reference 

As shown in the previous sxabsections , the required 
alternating (sine/cosine) current amplitude waveform for 
oscillatory gradient is generated by digital oscillator, and ix: 
reference for linear magnetic field gradient is derived from DC 
power supply. 
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(a) 


Jo 

2 






Fig.V. 3 Gradient selection using currents of different 

zunplitudes in a pair of coils, (a) -(d) : Different 
direct current (DC) amplitudes in left and right 
coils generate linear magnetic field gradients of 
different magnitude between the coils as shown 
(e) Alternating current (AC) amplitude in the 

right coil and no AC amplitude in the left coil 
places the 'selective plane' on the extreme left 
side, (f) AC amplitude 31^4 in the right, and 

complementary AC amplitude Iq/ 4 in the left coil 

places the 'selective plane' one-fourth the total 
distance from the left coil. 



The oscillator consists essentially of a digital 
square-wave multivibrator, digital serial shift registers with 
inverter in feedback and a resistive network to approximate the 
sine and cosine (sine wave phase shifted by n/2) wave forms. 
This staircase approximation is smoothed using a voltage 
controlled voltage source (VCVS) and a low-pass filter to 
obtain sine/ cosine waveforms. 

The detailed circuit, along with actual component values 
and the function of important components, is described below 
sequentially. As seen from Fig. V.4, a 555 timer chip is used to 
generate the basic clock waveform. This waveform is, however, 
rectangular in shape (i.e., duty cycle is not 50%). in order to 
get a 50% duty cycle waveform (square wave) the rectangular 
waveform is divided by two using a JK flip-flop (7476) . The 
square wave thus obtained has half the oscillating frequency of 
the 555 timer. This square wave is used to clock two pairs of 
digital serial shift-register IC's (74164). Sixteen registers of 
one pair of shift-register along with an inverter feed back IC 
(7404) shift out a sine wave, and sixteen registers of the other 
pair shift out a cosine wave. 

The digital pattern that circulates in the sine 
generating shift registers is set forth in Table V.l. This 
pattern must be converted by a resistive addition network to a 
staircase approximation sine wave. The values of resistors used 
by Fitzsimmons cannot be used here because the length of our 
shift register sequence is sixteen (compared to his twelve) . As 
the sine wave is odd symmetric, we therefore have to obtain only 
8 values of resistances (also known as amplification 
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Fig.V 4 Circuit diagram for the digital oscillator. 





















Table V.l 


Circulating digital pattern in sine wave generating 
shift registers. 
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in binary 


in Hex 


0000 

0000 

0000 

0001 

0001 < — 1 

0000 

0000 

0000 

0011 

0003 

0000 

0000 

0000 

0111 

0007 

0000 

0000 

0000 

1111 

OOOF 

0000 

0000 

0001 

1111 

OOIF 

0000 

0000 

0011 

1111 

003F 

0000 

0000 

0111 

1111 

007F 

0000 

0000 

1111 

1111 

OOFF 

0000 

0001 

1111 

1111 

OlFF 

0000 

0011 

1111 

1111 

03FF 

0000 

0111 

1111 

1111 

07FF 

0000 

1111 

1111 

1111 

OFFF 

0001 

1111 

1111 

1111 

IFFF 

0011 

1111 

1111 

1111 

3FFF 

0111 

1111 

1111 

1111 

7FFF 

1111 

1111 

1111 

1111 

FFFF-wnax 

1111 

1111 

1111 

1110 

FFFE 

1111 

1111 

1111 

1100 

FFFC 

1111 

1111 

1111 

1000 

. FFF8 

1111 

1111 

1111 

0000 

FFFO 

1111 

1111 

1110 

0000 

FFEO 

1111 

1111 

1100 

0000 

FFCO 

1111 

1111 

1000 

0000 

FF80 

1111 

1111 

0000 

0000 

FFOO 

1111 

1110 

0000 

0000 

FEOO 

1111 

1100 

0000 

0000 

FCOO 

1111 

1000 

0000 

0000 

F800 

1111 

0000 

0000 

0000 

FOOO 

1110 

0000 

0000 

0000 

EOOO 

1100 

0000 

0000 

0000 

COOO 

1000 

0000 

0000 

0000 

8000 

0000 

0000 

0000 

0000 

OOOO-wain ' 
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coefficients) associated with each bit. 

By solving the following equation for the nth clock 
pulse (n increases modulo 16) , 


+ »2<V2 * 


®16<Vl6 


= 1 + Sin 


fnrr Tr) 
[16 " 2) 


[V.3] 


where is the amplification coefficient corresponding to ith 

number of output signal of IC^s 74164 and i varies from 1 to 16, 
with the additional constraints that 
(Xj^)i V 1 s i s 16 


(x^)^ = 1 if i s n 
= 0 if i > n 


and a^ = V l s i s 16 


[V.4] 


[V.5] 


we have the following set of simultaneous equations : 

“ 1 + sin [Hg - 2 ) 

, , . f2Tr n ] 
ai+a2 = 1 + sin ^3^ - 2 J 

, . . (311 n ] 

1 [le " 2 J ' 

or, in general, 
n 

Y * 1 + sill “ f ) 1 ^ n s 16 . [V.6] 

i=l 

Using eqn. [V.5] in [V.6], 

I a. = 1 + sin [S| - 2 I for 1 a n s 8 

i=l 


[V.7] 
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Eqn. [V.7] yields a set of eight eguations, which are sufficient 
to determine uniquely the values a^ to ag. These can be 
suitably multiplied by a factor (we simply used a factor of 2) 
to obtain the required sine wave amplification value. We must 
keep one important point in mind while selecting the 
multiplication factor, that is, that the peak— to— peak voltage 
output of the waveform should not exceed +10V, at which voltage 
the output of a 741 Op. Amp. with Vss = + 12V normally saturates. 

The a^^'s obtained from eqn. [V.7] are given below in 

Table V.2. 

Table V.2 Amplification Coefficients for sine-wave Resistive 

Network 


Coefficients of Eqn. V.7 

lOOa^ 

^i 

^1 = 

0.0192 

100 

“2 = 

0.0569 

33.8 


0.0924 

20.8 

^4 = 

0.1244 

15.4 


0.1515 

12.7 

^6 = 

0.1729 

11.1 

^7 = 

0.1876 

10.2 

-8 = 

0.195*1 

09.8 


The numbers listed in the second column of Table V.2 are the 
resistance values in KO to be used for getting a staircase 
approximated sine wave; the resistance values chosen for our 
circuit are the corresponding rounded-off values. Our resistive 
adder network is given in Fig. V.5. The sine wave shiftout 
procedure from the digital serial shift registers (74164) 
followed by the above resistive network is given in Fig. V.6. A 
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Fig.V.5 



Stair case 
approximated 
sine wave 


A 



4 - 


Pirise no.(n) at input of 
^ift registers 


0001-0OFF~FFFF 0000- 


01 


0. — 0 


16 


Hex pattern at output 
of shift registers 


32 


Angle, ©n, where 

0n=('^^/l6-n^2) 


Fig,V.6 


Procedure for shifting out sine wave from digital 
serial shift registers (74164's). 
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little tuning is required for the resistive network to account for 
any round off errors in various steps of the calculation. The Op. 
Amp. 741 (in Fig. V.5) sums up the 16 bit resistive network 
outputs and produces a staircase approximated sine wave as its 
output. There should not be a clipping at either extremity of 
the waveform, with a maximum peak-to-peak excursion of about 9V. 
In case clipping is observed, the feedback resistor (1.88 KQ) in 
the summing operational amplifier may be further reduced till the 
peak-to-peak excursion reduces sufficiently to become 
distortionless . 

At higher frequencies, the waveform may be distorted; 
in that case all the resistors in the network should be increased 
by a factor of 5 to 6 or even more in order to realise undistorted 
waveforms. 

The next stage is a low-pass filter (Fig. V.4) which 
removes the high frequency component of the staircase sine wave 
and gives smoothed sine wave at its output. In order to carry 
this to the other parts of the controller, it is current-buffered 
by a 741 Op. Amp. in 'voltage-follower' configuration because of 
its superb buffering qualities, extremely high input impedance, 
low output impedance and unity voltage gain. 

Cosine wave is generated by coupling its shift register 
to the shift register of sine wave appropriately. Sine and 
cosine outputs are digitally phase locked regardless of the 
frequency of operation. The sine and cosine outputs allow 
gradients to be controlled in quadrature when this is needed. 

The frequency of the sine/cosine wave can be varied 
easily by varying the 100 KQ potentiometer connected to the 555 
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timer (see Fig. V.4) within a fixed range depending on the value 
of the capacitor marked with asterisk in Fig. V.4. The frequency 
range can be increased further by simply replacing it with a 
suitable capacitor; in our circuit, the range is 1-25 Hz with 0.1 
jiF and 12-300 Hz with 100 nF. 

V. 3.3 Computer Interface and Manual Switch to Gradient 

Controller 

We have a twofold objective for designing a computer 
interface and manual switch for gradient controller, namely : 

(i) to select at least 256 (numbered between 0 and 255) 
'selective planes', as shown in design logic subsection, for 
imaging purposes; 

(ii) to vary the magnitude (in T/m) of the continuous/pulsed 
linear magnetic field gradient. 

To achieve these objectives via microprocessor (uP) 
control, we used two 4-bit binary counters (7493) cascaded to form 
an 8-bit serial to parallel converter shown in Fig. V.7. As 
argued by Fitzsimmons , this requires only two pulse lines (one 
for reset and another for outputting plane position number of 
pulses) from microprocessor, whereas a parallel design would 
require as many as eight pulse lines from the microprocessor. 

The nimber of pulses which is equal to the number of 
planes to be selected, or the magnitude of the required linear 
magnetic field gradient serially coming from the microprocessor, 
are converted into a binary form by this 8— bit binary counter 
(serial - to - parallel conversion) . A reset pulse is required to 
initialize the counter. 
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Display Circuits 


ig.V.7 


Circuit diagraai of our computer interface and 
manual switch to gradient controller. 
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The same objectives are also achieved in manual mode by 
setting an 8 -bit DIP switch in either on or off position to 
select the required binary number which is displayed in decimal 
form. 

The selection of computer or manual control of gradients 
is done by a set of multiplexers (74257) which takes either 
computer outputs or manual switch outputs as its input depending 
upon the position of the toggling switch (SWl/1) . 


V. 3 . 4 Circuit Generating Complementary Outputs 

The purpose of this circuit is to generate two DC/AC 
(sine or cosine) signals needed for a pair of gradient coils 
with the condition that the current amplitude of these two 
signals should complement each other. Our circuit generates a 
total of 256 complementary outputs between 0 and 255. The 
resulting outputs, which have the following relationship (see Fig. 
V.3), 


II 


coil 1' 


1 1 


coil 2 


I = |lol 


[V.8] 


are fed to the gradient coils. The corresponding circuit, 
represented in Fig. V.8, consists of a set of Digital-to-Analog 
converters (DAC) and four 741 Op. Amp's. The inputs to the DAC's 
are : 

(i) 8-bit binary data from the multiplexer (see Fig. V.7), 

(ii) sine wave from the digital oscillator or -12V DC from 

power supply, as its negative reference, . 
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Complementary outputs circuits 











182 


The output of DAC's is the same as their input with only 
the amplitudes dependent upon the binary input data (also called 
as attenuation factor) . One DAC receives the 8-bit binary data 
from the multiplexer (Fig. V.7) directly and another receives the 
inverted multiplexer-output so that the attenuation factor of the 
DAC's will have a complementary relationship. 

Output current (^^nalog^ given by 


I 


analog 


^ref(-) f ^ + 
^5 I 2 



[V.9] 


where binary data, A^^ = 1 (A^ at high level) , 
and ^n ~ ^^n level) , 

^ref(-) ~ negative reference voltage (sinewave or 
-12V DC) 

R, _ = resistance at pin no. 15 of DAC IC MC1408. 

15 

The negative is selected here because the digital 

oscillator circuit output is negative due to the current sinking 
property of transistor-transistor-logic (TTL) shift register 
(74164) . 

Since the output of the DAC is in current mode, it is 
converted into voltage mode by a 741 OP. Amp. In the case of sine 
wave, the output of Op. Amp. is DC-blocked by using a 10 4 F 
capacitor and then fed to the voltage gain amplifier (741 Op. 
Amp.), but in the case of DC output it is directly applied to the 
voltage gain amplifier. 



V.3.5 


Switch Circuit for Pulsed 
Gradient Generation 


Linear Magnetic Field 
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The linear magnetic field gradient is achieved by simply 
replacing the sine wave reference given to the DAC (pin number 15 
of MC1408 in Fig. V.8) with a -12V DC. The magnitude of the 
gradient can be varied by varying the 8-bit data input to DAC's. 
It can also be made to rotate by applying 'reset' and 'position' 
pulses continuously from the microprocessor. 

The linear magnetic field gradient thus generated should 
be gated at a suitable point before feeding it to the gradient 
coils in order to achieve a pulsed linear magnetic field 
gradient. We place a simple JFET switch (LF11333) at the output 
of the 741 Op. Amp. in voltage gain configuration (see Fig. V.8). 
This switch can be either microprocessor-controlled, or set 
manually by applying +5V DC. The outputs of this IC are fed to 
current drivers, described below. 

V.3.6 Current Driver Circuit 

We have used a fairly conventional design with a 
single Op. Amp. (741) providing the input to the driver and 
complementary pair of power transistors. The complete circuit 
diagram is shown in Fig. V.9, and consists of a 741 in negative 
feedback mode, transistors (2N 3055, ECP 055) and two diodes 
(1N914) . The negative feedback path has a push-pull cimplifier 
comprising a 2N 3055, ECP 055 and one 2N 1183. The l£2 feedback 
resistor (marked with * in the circuit) in the circuit provides a 
suitable point for measuring the current being delivered to the 
gradient coils since it is in series with them. The gain can be 



Current driver circuit 


Fig.V.9 
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varied with the help of the 3Kn potentiometer in the feedback 
path. 


V. 3.7 Display Circuit 

The binary number selected for generating complementary 
outputs can be displayed conveniently in two ways, namely, (i) in 
a binary form and (ii) in a decimal form. 

(i) Binary Display 

This is done by simply connecting eight Light Emitting 
Diodes (LED's), one for each bit, in the 8-bit data output of the 
multiplexer (see Fig. V. 10) . The common anode of all LED's is 
given +5V DC, which enables them to display the complement of the 
binary number selected, in spite of the fact that the LED's are 
not present in the inverted multiplexer-output line. 

(ii) Decimal Display 

Since it was inconvenient to see the selected number in 
an 8-bit binary form, we decided to display the same in a decimal 
form. 

To accomplish this, the 8-bit binary data from the 
multiplexer is first converted into three nibbles (nibble = a 
group of four bits) of Binary Coded Decimal (BCD) by a set of 
three binary-to-BCD converters (74185), shown in Fig. V.IO. Each 
of the resulting BCD nibbles is given to a BCD-to-7 segment 
decoder (7448) , which is an active 'high' requiring pull-up 
resistors for converting the BCD to a decimal niunber. The seven 
output lines of 7448 are connected to a common cathode 
seven-segment display (LTS543R) . 
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From Multiplexers 



Fig. V. 10 


8-bit binary LED and 3 -digit decimal seven segment 
display circuit. 
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V.3.8 Position of Switches for Mode (^tP/manual) and Sine 

Wave/DC Selection 

There are a total of seven switches in the circuits 
which are to be put in appropriate position to select mode 
(fzP/manual) and reference signal (sine wave/ DC) to DAC. These 
are grouped into two and given in Table V.3 for easy visualization 
and generation of linear/oscillating magnetic field gradient. 


V. 4 Details of the Microprocessor and Its Use in Gradient 

Control 


V. 4.1 Description of the Microprocessor 

The microprocessor used in our work is a low cost Intel 
8085 A based 8 -bit system (Microfriend-I, Dynalog Microsystems 
(DMS) , Bombay, India) . Its functional block diagram is shown in 
Fig. V.ll. 

The 8085 A processor chip contains six general purpose 
8-bit registers B,C,D,E,H and L for both data transfer and data 
storage, and an 8-bit accumulator (register A) for both 
arithmetical and logical operations. With 6.144 MHz crystal, the 
CPU operates at 3.072 MHz i.e., 1-3 ns instruction cycle - 

The microprocessor (M?) lias a powerful system monitor in 
2732 Erasable Programmeible Read Only Memory (EPROM) of 4K bytes. 
The on-board 24 pin zero insertion force (ZIF) EPROM Programmer is 
driven by a Programmable Peripheral Interface (PPI) 8255A to burn 
the developed software program into either 2716 or 2732 EPROM 
from Random Access Memory (RAM) . The powerful on-board system 
monitor program enables us to load any software program into 
on-board RAM through a 21 -key keyboard or through an audio 



188 




Switch 


V.7 SWl/1 


Mode 

(PP/ 

manual) 


V.8 SWl/2 



To 

JFET 
LF 11333 
(Pin 1,8) 


SW2/1 

SW2/2 (MCU08, 
Pin 15) 


1-8 KA 


Manual 


OPEN 


>To MUX 74257 


Manual 


From pp 


• ♦5V DC 


Cine fS'ne wave 

o ..-qJ from digital 
f ^oscillator 

DC 

o — mth-o -12V DC 
6.2 KA 


Sine 

wave/ 


From 

SW2/3 OP-Amp.741^ 

V-8 (Current to 

SW2/4 voltoge 
converter) 


I ©OPEN 

DC 

o To Op- Amp. 741 
10 KA (Voltage <3aln) 


From 

JFET 

V.9 SW2/5 Switch 

(Pin 2/7) 


oOPEN 

DC 

■wvw ■ o To Op. An^. 741 
4.7KA (Driver A/B) 


22pF 


TalDle V.3 


Tabulation of circuit switch positions, 
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STd bus 



/ 

K 


7 

/ 

— 

V 


N 

1 

L 


A 

/ 


8085 A 
CPU J 


8155 Timer used 
by system as a 
single stepping 
operation 

(i) 256x8 scratch 

pad RAM 

(ii) 6bit I/O Port 

(iii) 2x8 bit I A) Port 

(iv) Programmable 
14 bit binary 
counter/Timer 

sid' 


^SOD 


Audio 

casette 

Recorder 

or 

TTY 

or 

CRT 




( 



22 lA) Lines 



24 lA) Lines 


Functional block diagram of Intel 8085A based 
microprocessor kit used in this work. 


Fig. V. 11 
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cassette interface or from ZIP EPROM. An on-board 2K RAM (6116) 
user memory is suitable for developing and executing software 
programs. Intel 8279-5 is a programmable keyboard display 
interface IC which is used for a 21-key keyboard and 7-segment 
display. The display consists of a four byte memory address and 
a two byte data. 

The IC 8155 comprises a 256x8 bit static MOS RAM with 
I/O ports and timer. The system monitor uses up about 100 bytes 
of static RAM as CPU scratch pad and stack area. The remaining 
156 bytes are available for the user. The timer is used by the 
system for single stepping operations which are useful for 
debugging the developed machine code program. 

The address, data and control buses fully buffered (by 
IC's 8212, 74LS244 and 74IiS245) are available as an STD bus, which 
is a 56 pin double-sided edge connector provided for expansion of 
the system. Through this STD bus is hooked up a Real Time 
Controller (RTC) card mounted in card cage. 

The RTC card consists of a 8255A-5 PPI, 8253C-5 
programmable interval timer, 8251A Universal Synchronous 
Asynchronous Receiver and Transmitter (USART) and 8259AC 
Programmable Interrupt Controller (PIC) . 

The PPI 8255 A-5 has three 8-bit I/O ports, port A, Port 
B and Port C. Each Port and control register have distinct 
addresses to set control word in control register and to 
communicate with different devices through each port. This 8255 
A-5 PPI can be set in various modes by writing suitable codes in 
the control register. The available modes are : Bit Set Reset 
(BSR) mode, simple I/O mode, handshake mode and bidirectional 
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mode. We use the 8255 a- 5 in a simple output mode to output reset 
and position/DC-magnitu( 3 e pulses and gate-control through BNC 
connectors to gradient controller box for selecting the plane or 
increasing DC magnitude and gate the output. 

V.4.2 Reset and Pcksition/DC-Magnitude Pulses and Gate Control 
Generation 

V. 4.2.1 Initialization 

the PPI jjorts are set in output mode by writing the 
appropriate word in control register- The number (N) 
corresponding to the first plane/DC-magnitude to be selected is 
stored at memory loca.tion M(X1) . Next, a reset pulse is sent at 
BNC-1, which sets gradient controller at zeroth 
plane/DC-magnitude fof a short computer time. Bringing the 
contents of memory location M(X1) in E-register is the next step. 
Here E-register is used as a 'scratch memory' for count down. 

V. 4.2.2 Software Counter 

As explained above, one pulse is sent at BNC-2 and the 
BNC-3 level is set high for gating. The number of pulses is 
brought out from E-register for countdown and the decremented 
number is restored back in the E-register. If countdown is not 
over, the program loops back to pulse generation at BNC-2 stage. 
During each loop a pulse is sent through BNC-2 to gradient 
controller which in turn increments the plane/DC-magnitude by one. 
After the set number of pulses at M(X1) have been sent, the 
countdown ends;- the desired plane/DC-magnitude is selected, and 
the program j\imps to the system monitor program. 
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V.4.2.3 Pulse Output Method 

There is one-to-one correspondence between the eight 
accumulator bits and 8— bits of each port. We set a particular bit 
of accumulator high (and other bits low) corresponding to a bit of 
a port from where the pulse output is taken. The content of the 
accumulator is output to PPI port, say port A. This makes the 
selected pin of port A high. Next, accumulator is immediately 
reset to zero and the content is output in the same port. This 
brings the selected pin of port A level low. This method 
generates about 5tis pulse at selected pin of port A. In the case 
of gating, level is left high for turning on the JFET switch. 


V. 4.2.4 Successive Increment of Plane/DC -magnitude 

Our fiV system has a user key U1 which, when tapped, 

takes the program from system monitor to scratch pad memory 

location 209C. At this memory location we write a jump 

instruction to jump into our program area. Here we have a program 

to read a number of pulses from M(X1) and then increment it by a 

fixed number (N ) and store the result again at M(X1) . After this, 
o 

the program enters the initialization part of the main program. 
Every time the U1 key is pressed it selects/ increases a 
plane/DC— magnitude higher by the fixed number (N^) . In this 
incrementation process, the 'carry' is discarded (i.e. , the 
number rotates between 0 and 255) . 

The pulse generation details described above are given 
below in the form of a flow chart. Fig. V.12,- a sample program for 
generating reset and position/DC-magnitude pulses and gate control 


is given in Table V.4. 




Set all ppi 8255 ports os 
out put ports 

V 

Store number of pulses (N) at 
memory location M(J<1) 


r 

Send reset pulse at BNC-1 


— Jl 

Store number of pulses in 
E- register 



V 

Send one pulse ot BNC-2 
and level high at BNC-3 


1 

Bring number of pulses from 

E- register to occumulator 
and count down and restore 
it in E- register 



To monitor programme 
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Reset 


_n 


Posit ion /DC 

magnitude pulse 

Gate control — 


1 2 - 





.V.12 


Flow chart for the generation of reset and 
position/DC“inagnitude pulse (s) and gate control. 














Table V. 4 


A sample program for generating reset and 
position/DC-magnitude pulse(s) and Gate Control 
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Machine Code 


1000 

MOVI A, 80 
OUT 2F 

3E 80 \ 
D3 2F J 

1004 

MOVI A, 01 
STA 1012 

3E 01 \ 

32 12 lOj 

1009 

MOVI A, 80 
OUT 2C 

MOVI A, 00 
OUT 2C 

3E 80 

D3 2C 

3E 00 

D3 2C ^ 

1011 

MOVI E, 01 

IE 01 } 

1013 

MOVI A, 20 
OUT 2C 

MOVI A, 10 
OUT 2C 

3E 20 

D3 2C 

3E 10 

D3 2C ^ 

lOlB 

MOV A, E 

SUI 01 

MOV E, A 

7B 

D6 01 

5F 

lOlF 

JNZ 1013 

C2 13 10} 

1022 

RST 0 

Cl > 

1022 

JMP 1009 

C3 09 10 ' 

► 

J 

209C 

JMP 1070 

C3 70 10} 

1070 

LDA 1012 

ADI 10 

STA 1012 

3A 12 10' 
C6 01 

32 12 10, 

1078 

JMP 1009 

C3 09 10! 


Comment/R ema rlc 


PPI in output mode 


Number of pulses (N) to memory 
location 1012 (i.e., M(X1) 

A reset pulse is sent at 
BNC-1 (bit 8 of port A is 
used) 


No. of position/DC-magnitude 
pulses is transferred to 
E-register 

For sending position/DC- 
magnitude pulses at BNC-2 
(bit 6 of port A is used) 
and Gate control at BNC-3 
(bit 5 of port A is used) 


(Position/Dc-magnltude -1) 
in E-register 


Loop back if count down is 
not over 

For monitor, ready for U1 key 

These instructions are for 
displaying reset and position/ 
DC-magnitude pulses and gate 
control on oscilloscope, with 
scope triggred by reset pulse 

When U1 key is pressed, 
system monitor enables a jump 
to 209C. Here we provide a 
jump to program area 1070 

Read number of position/DC- 
Magnitude pulses from 1012 
(i.e. M(X1), increment it by 
a fixed number (N ) stored at 
1074 and store back 

Jump in initialization part 
of program where reset pulse 

is generated. 


pe changed as required^ 


Note : Underlined numbers can 
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V.4.3 


Gradient°?PMFGf Pulsed Magnetic Field 

Gradient. (PMFG) Spm Echo Experiment 


All PPI ports are set in output mode by writing the 
appropriate word in control register. Various delay counts (DCl, 
DC2 and PI) are stored in memory locations M(X1) , M{X2) and M(X3) . 

The pulse output method here is the same as described 
in section V.4.2.3. However, delays are introduced to lengthen 
pulse width (PI) and to introduce intervals between pulses. 

For short pulse width, 8-bit number is brought into the 
accumulator and counted down. Thereafter, the level at BNC-3 is 
changed. This generates pulse widths dependent upon an 8-bit 
number . 

For large pulse interval delays, two 8-bit nxombers are 
used. For each decrement of the second 8-bit number, the first 
8 -bit number goes to complete count down. By this method we get 
delays dependent upon multiplication of two 8-bit numbers. 

The details described above are given below in the form 
of a flow chart. Fig. V. 13, followed by a sample program in Table 
V. 5 


V.5 Highlights of Our Gradient Controller Circuit 

In this section we summarise the technical features of 
our improved gradient controller circuit; improvements over the 
earlier design of Fitzsimmons are also pointed out. 

(1) In digital oscillator circuit (Fig. V.4), 

(i) after the 555 timer, a JK flip-flop is added to get a 
square waveform with 50% duty cycle. 
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START 


Bring count down DC 1 from 
M(xl) to occumulotor 


Set all PPI 8255 ports os 
out put ports 


Store delay count DC 1 at memr 
ory location M(xi) for delay 
T1 = DC1«OF 


— — — ^unt dowfK 
^over’/ 

^ES 

Level high at BNC-3 


Store deloy count DC2atmer 
ory location M(5<2) for delay 
T2 = DC2*OF 


Store delay count PI at memo- 
ry location M(x3) 


Level high at 

BNC-3 


f 

Bring PI from M(x3) to accu- 
mulator and count down 

■ J 



^ount dowr 
\pver?^ 


Bring count down PI from 
M(x3) to accumulator 


-^unt do^ 
x^ver^/' 

^ES 

Level low at BNC-3 


I Bring count down OF to accu- 
I mutoter 


count dowr 
\over?> 


Level low ot BNC-3 




1 j 

Bring count down DC 2 from 




- 


M(x2) to accumulator 

T 


Bring count down OF to accu- 
mulator 


count dowm 
\over?/ 


aount dowm 
^Xpver?/ 


Fig. V. 13 


Flow Chart for the generation of switch control for 
pulsed magnetic field gradient (PMFG) spin echo 
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Table V. 5 


A sample program for generating switch control 
a typical pulsed magnetic field gradient fPMFG) 
echo experiment v 


for 

spin 


Address 

Mnemonic 

Machine Code Coitiment/Reiriiark 

1030 

MOVI A, 80 
OUT 2F 

3E 80 

D3 2F 

\ All ppi ports in output 

/ mode 

1034 

MOVI A, 02 
STA 1091 

3E 02 

32 91 1C 

1 Delay count DCl stored at 

)/ memory location 1091 

1039 

MOVI A, 04 
STA 1092 

3E M 1 Delay count DC2 stored at 

32 92 lOj memory location 1092 

103E 

MOVI A, 04 
STA 1093 

3E 04 

32 92 10 

\ Delay count PI stored at 

/ memory location 1093 

1043 

MOVI A, 10 
OUT 2C 

3E 10 

D3 2C 

\ Level high to turn on JFET 

/ switch (bit 5 of port A) 

1047 

LDA 1093 

SUI 01 

JNZ 104A 

3A 93 10 
D6 01 

C2 4A 10 

1 Count down PI 

104F 

MOVI A, 00 
OUT 2C 

3E 00 

D3 2C 

) Level low to turn off 

1 JFET switch 

1053 

MOVI A, OF 
SUI 01 

JNZ 1055 

LDA 1091 

SUI 01 

STA 1091 

JNZ 1053 

3E OF 

D6 01 

C2 55 10 
3A 91 10 
D6 01 

32 91 10 
C2 53 10 

J 

Count down T1 
where T1 = OF x [DCl] 

1065 

MOVI A, 10 
OUT 2C 

3E 10 

D3 2C 

1 Level high to turn on 

JFET switch 

1069 

LDA 9310 

SUI 01 

JNZ 106C 

3A 10 93' 
D6 01 

C2 6C 10^ 

- Count down PI 

1071 

MOVI A, 00 

OUT 2C 

3E 00 \ 
D3 2C / 

Level low to turn off 

JFET switch 

1075 

MOVI A, OF 

SUI 01 

JNZ 1077 

LDA 1092 

SUI 01 

STA 1092 

JNZ 1075 

3F OF 

D6 01 

C2 77 10 
3A 92 10 ^ 
D6 01 

32 92 10 
C2 75 10 

Count down T2 
where T2 = OF x [DC2] 

1087 

JMP 1034 

C3 34 10} 

Jump back to reinitialize 


delay count and give waveform 


Note : Underlined numbers can be changed as required 
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(ii) Four TTL 8-bit parallel out serial shift registers 
(74164 's) have been used instead of three CMOS dual four 
stage registers (4015 ^s) to shift out sine and cosine 
waves . 

(iii) The sine/cosine wave resolution in our circuit is (32 x 
input period to 74164) , as contrasted from Fitzsimmons 
(24 X input period to 4015). 

(iv) To get a smoothed sine/cosine wave, we have used a very 
simple single pole low-pass filter instead of using an 
MC 1408 current sync. IC. 

(2) In the computer interface and manual switch circuit (Fig. 
V.7), we have used two 74257 (quad 2 input multiplexer 
tristate) instead of two 74367 (tristate hex buffer) which 
makes DIP switch connections easier, and operation in both 
modes simpler. 

(3) In the complementary outputs generating and switching circuit 
(Fig. V.8), 

(i) We feed the smoothed sine-wave/DC to negative reference 
of DAC's; 

(ii) .A JFET switch LF 11333 is used in the input to the 

current driver circuit Op. Amp. (Fig. V. 9) instead of 
the CMOS transmission gates suggested by Fitzsimmons^ ^ . 

(4) We have used 2N 1183, 2N 3055 and ECP 055 transistors instead 
of ECG 185, MJE 3055 T and MJE 2955 T in the current driver 
circuit (Fig. V. 9) . 

We have added a decimal seven segment LED display for showing 
the selected plane/DC-magnitude level (Fig. V. 10) . 


( 5 ) 
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In the photograph, Fig. v. 14, we show the experimental 
setup of our current driver, gradient controller, microprocessor 
and RF probe units. 

Details of Experimental Setup and Performance Result 

The versatile magnetic field gradient controller 
described above has been constructed with the idea of dedicating 
it to our microprocessor-controlled low frequency (15 MHz) home 
made NMR spectrometer . Unfortunately, during the last stages 
of our work, our electromagnet and its power supply became 
non- functional . 

As a result, a commercial spectrometer (BRUKER WP-80) 
had to be used to complete our experiments. Even here, we were 
constrained by the fact that the usual experiments could not be 
done without modifying the existing hardware of the WP-80, which 
thus railed out the possibility of doing diffusion measurements on 
the WP-80. With these restrictions, we limited ourselves to 
checking the performance of our versatile magnetic field gradient 
controller. 

Details of the gradient coils used, and of experimental 
set up, are given below : 

Spectrometer ; BRUKER WP-80, with BNC-28 Computer 

Probe used ' WP-80-^^C/^H - - ^H, 5mm 

Gradient Coil\ : OHFC Copper wire thickness : 0.6 mm 

specif icationj 

Coil outer diameter : 13.8 mm 
Number of turns : 10 
Coil used : Miniductor coil No- 3004 
(From Barker & Williamson, Bristol, Pa., 
U.S.A.) 
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Fig. 14 


Es^erimental setup of the current driver, gradient 
cx>ntroller, microprocessor kit and RF Probe unit. 
(Photograph) 
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The gradient coils were fixed rigidly by supporting them 
with batts of thermocol (styrofoam) in the empty space on each 
side of the sample coil (along y-direction) in such a 
configuration that it generated y-gradient. The y-gradient was 
chosen because there was no empty space in the probe box for 
fixing coils to generate a z-gradient. 

( i ) Linear Magnetic Field Gradient : 

For checking this, we have selected a 5mm NMR tube 
containing D 2 O (this is for field-to-freguency locking) , into 
which two H 2 O capillaries of different diameter were inserted. 

The experimental spectra corresponding to four different 
conditions are given in Fig. V. 15 and Fig. V.16. Due to small 
voltage offset difference between the Op. Amp.'s of the current 
driver circuit, the spectra shown in Fig. V.15(a) and (b) are not 
exactly symmetric. Spectrum in Fig. V. 16 (b) shows that when 
the DC-magnitude 127 is selected, both the coils receive the same 
amount of voltage and this imposes a small additional magnetic 
field in the y-direction, as compared to the situation in Fig. V. 
16(a) where no current is passed through the coils. The 
broadening in Fig. V. 16(b) may be explained due to the field 
inhomogeneity generated by the current flowing through the 
gradient coil pair. 

Fig. V. 17 shows spectra at different orientation of the 
capillary tubes with respect to the magnetic field gradient 
direction. 
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(a) 




LOCK ON 

No spinning 
No Gradient 
No. of FID = 1 


(b) 




Proton NMR spectra of H^O in (a) absence and (b) 

presence of a small additional magnetic field in 
the Y-direction. 


Fig. V. 16 



4 



Fig. V. 17 
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(ii) Oscillating Field Gradient 

For checking this, we have used a sealed D^O tube, a 
small amount of H 2 O present there being sufficient to give a 
strong NMR line. Figs. V.18 and 19 have been recorded at two 
different oscillating frequencies. As expected, multiple lines 
are observed in each spectriam due to magnetic field modulation at 
the oscillating frequency. 

V.7 Self -Diffusion Coefficient Measurements on the Lecithin 

Reverse Micellar *Pool Water* 

V. 7 . 1 Motivation 

Self-diffusion is the net result of thermally induced 

random-walk process experienced by molecules in solution, and 

self-diffusion coefficients provide important information on 

aggregate structure and association behaviour because association 

of molecules in solution greatly affects the translational motion 

of the molecules. Indeed, it has been well established that 

self-diffusion coefficients provide a direct approach for 

elucidating the structural organization in aggregated 

systems . Typical self-diffusion coefficients in liquid 

—9 2 —1 

systems at room temperature range from about 10 ms to 
10 “^^ m^s“^. Elsewhere in the literature, measurements of the 
diffusion coefficient of cell water in biological systems have 
disclosed reductions of as much as 50% in the diffusion 
coefficient of this state of water compared to that of free 

water^^^^ ; likewise, in other 'structured' phases such as water 

containing microemulsions ^ ^ , the diffusion coefficient of water 
has been found to be less than 10% of that of free water. 




Fig. V. 18 Proton NMR spectrum of water at the gradient 

oscillating frequency of 50Hz.(V **6V; plane 

Jb^ 

number selected - l27;no.of FlD's collected = 128. 
The magnetic field lock was on, and sample tube was 
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Self-diffusion coefficient can be measured by using the 
following techniques : 

(i) Tracer techniques'^®' 

(ii) Nuclear spin echo method^®^^, 

(iii) Transient optical grating methods , 

(iv) Neutron spin echo method 

(V) Gouy interferometric techniques . 

Our interest is to measure self-diffusion coefficients 
by the nuclear spin echo method, which is therefore described in 
further detail below. 

V.7.2 The Nuclear Spin Echo Method 

Translational self-diffusion coefficients (D) in single 
as well as multiphase systems are measured by observing the loss 
of phase coherence of a set of processing nuclear spins due to 
their random motion over a gradient in the magnetic field. This 
method is the basis of the measurements presented in the present 
section, and a brief review of the principle of the method shall 
be presented first. 

As is well known^^^'^®^ , a general method for measuring 
transverse relaxation times is by studying the decay of the 
transverse component of magnetization, i.e., the component in the 
XY -plane when the Zeeman field is along the Z-direction. 
Hahn^^^^ , in a classic early paper, has pioneered the use of the 
90 °-t- 180° RF pulse sequence for measuring T 2 . Fig. V. 20 (i) 
illustrates the vector model which is usually employed to 
describe the process. 
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A 90 pulse causes the macroscopic magnetization, 
M, to tip to the y* axis. This magnetization loses phase 
coherence as there is a spread of frequencies in the sample owing 
to the spin-spin relaxation process (time constant T^) and 
inhomogeneities of the magnetic field, some nuclear moments 
precessing more quickly than average, some more slowly. The 
signal detected along the y'axis by a coil set in the xy-plane, 
the free induction decay (FID) , therefore decreases in magnitude. 
A 180° pulse applied at a short time, t, after the initial pulse 
reverses the moments by rotation about the x'-axis. Phase 
coherence is first regained, and subsequently lost. Thus the 
signal appears as a reversed FID followed by a second pulse to 
give the echo. The echo height is less than that of the FID, and 

I 

the echo is slightly asymmetric, owing to the spin-spin relaxation 
and diffusion during the time interval of the experiment, 2r, but 
the additional loss of phase due to the imperfect field, which 
would otherwise have occurred, is eliminated by this Hahn pulse 
sequence. Therefore, measurements of spin echo amplitude for 
different values of t provides a means for determining 

In Hahn's method, after each point is recorded it is 
necessary to wait for the magnetization to return to equilibrium 
(say s 5 Tj^) before the next pair of pulses is triggered. 

A variation of Hahn's spin echo method has been 
developed by Carr and Purcell in their method, a 90° pulse 
is followed by a train of 180° pulses at times r, 3 t, 5t, etc. 
The major drawback of this sequence is that the cumulative effect 
of the out-of-plane refocussing due to imperfect 180° pulses can 


be serious. 
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To overcome this difficulty, Meiboom and 
further modified the Carr-Purcell sequence. Here, the 180° pulses 
are not applied in phase with the 90° pulses, but rather with a 
phase shift of 90° (Fig. V. 20(ii)). This modification of a 90° 
phase shift after the first 90° pulse renders the precise setting 
of the 180° pulse width less critical. The phase shift also 
eliminates the effects caused by inhomogenities in the RF field. 
At present, the common pulsed NMR method for T 2 measurements, and 
perhaps the most accurate, is the one utilizing the Carr-Purcell- 
Meiboom-Gill (CPMG) pulse sequence. 

The same spin echo experiment, done in an inhomogeneous 
magnetic field (generally linearly inhomogeneous in one direction 
only) , provides the basis for NMR self— diffusion experiments. The 
90 °-t- 180° spin echo experiment refocusses spins of any precession 
frequency, but only if the precession frequency is constant during 
the experiment. Therefore, if the nuclei move in an inhomogeneous 
magnetic field, their frequencies vary as well. For coherent 
motion in the direction of the gradient, the spins just refocus 
completely in another direction. But for the incoherent random 
motions such as self -diffusion, a random phase shift of individual 
spin vectors occurs, resulting in incomplete refocussing at the 
time of the echo. This is best demonstrated pictorially by 

Callaghan , as in Fig. V.21. 

In the figure, the positive and negative sign in the 

triangle indicate the successive, random jump (translational 
motion) of molecules with respect to the gradient direction. The 
horizontal sum of these elements gives the phase shift which 
occurs in the mean time between jumps (t^) while the spin is at 



time 
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jump sign ^ 


t =0 




2t=2nT^ 
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480* RF pulse 






Phase shift 
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Fig. V. 21 


Cumulative phase diagram for spin echo 
a steady magnetic field gradient . 


formation in 
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current location. The net phase shift is the sum of those 
occurring before and after the 180^ RF pulse, and corresponds to 
the dotted region. 

The self-diffusion coefficient of the molecule in which 
the nuclei are situated is calculated quantitatively from the echo 
attenuation effect, according to the equation 

M(2t) = M(0) [^exp (- 2 T/T 2 - I r^DG^T^)j [V.IO] 

where y represents the gyromagnetic ratio of the nuclei, G the 
strength of the linear magnetic field gradient in T/m, and D is 
the diffusion coefficient in m /s and t the interval between the 
90° and 180° pulses. T 2 is the (natural) transverse spin 
relaxation time of the nuclei. Since T 2 can be measured rather 
precisely by using CPMG sequence, and since G is known (by 
calibration) , D can be determined directly from a plot of log 
[M(2 t) + 2 T/T 2 ] versus . 

Use of large currents in the gradient coils for 
producing linear magnetic field gradient has the following 
problems: (i) heating of gradient coils and (ii) NMR linewidth 
increases with corresponding decrease in the width of the FID and 
echo. Therefore, the detector bandwidth must be increased, 
thereby decreasing the signal— to— noise ratio. This linear 

I 

magnetic field gradient teclinigue has a diffusion coefficient 
limit of about lO^^^m^s""^, which is set by the maximum field 
gradient that can be used. 

(72) ! 

To overcome this limitation, Stejskal and Tanner i 

have proposed the pulsed gradient spin echo method (PGSE) which j 
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involves transient field gradient pulses of a few milliseconds. 

The amplitude of these pulses must be of the order of a T/m in the 

—12 

case of the small diffusion coefficients encountered in the 10 
- lo”^® m^s”^ range. 

The amplitude of the echo signal detected is given by 
M(2T) = M(0) j^exp [-St/T^ - (rSG)^ D(A- |)jj [V.ll] 

where 5 is the duration of the gradient pulse and A is the delay 
between two gradient pulses. 

By Fourier transforming the second half of the echo 
signal the pulsed field gradient spin echo method can be 
adapted to high resolution experiments, thus allowing the 

measurement of the diffusion coefficients of different components 

, (75,76) 

in a solution ' 

Interest in self-diffusion coefficient measurements has 

grown rapidly during the last fifteen years. The information that 

can be obtained from the self-diffusion behaviour of surfactant 

systems may be classified into ; 

. . (77) 

(i) surfactant micellization 

(78 ) 

(ii) solubilization equilibria 

(79) 

(iii) micellar breakdown by solubilizates 

( 80 ) 

(iv) relations between diffusional rates and phase structure 

(V) micelle size, shape, hydration, and intermicellar 

. ^ ( 81 ) 
interactions 

Excellent reviews by Stilbs on the FT PGSE studies of 
molecular diffusion*®^’ , mesophase structure of membrane lipids 
by Lindblom<«>, micellar solutions by Chachaty , reversed 
micelle by Vos et al.^®^\ and microemulsions by Lindman and 
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Stilbs have been published. 

From the results of previous studies in this 
laboratory there is strong empirical evidence to suggest 
that, in the lecithin/cyclohexane/H^O reverse micelle, at R 
values in excess of 6, a marked change occurs in the arrangement 
of the phospholipid head groups as well as in the dynamic state of 
water trapped in the core structure of the micelle. These changes 
have been rationalised as being due to the possible onset of a 
non-isotropic phase^®^~®^^ . 

In the following subsections, the above mentioned 
structural transition in the lecithin/cyclohexane/H.O is further 
probed and compared with lecithin/CCl^/H^O by measuring the 
self-diffusion coefficient of water over a range of R values. 

V.7.3 Hater Diffusion in Anisotropic Environment 

For water diffusion in an anisotropic environment, Eqn. 
V.IO has to be modified . In the specific instance of a liquid 
crystalline or two dimensional phase, the diffusion of trapped 
water will be mainly parallel to the long axis. An angular 
dependence is then imparted into the diffusion tensor (which will 
now be cylindrically symmetric) according to the formula 

Daniso = < D^> + (< D||> - < D^>) ain^e tv. 12] 

where < 0||> < Dj^> the components of the average diffusion 
tensor parallel and perpendicular to the long cixis of the 
anisotropic phase, and 6 is the angle between the field gradient 
direction and the normal to the preferred axis for diffusion . 
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Since, for our assumed model, < D„> » < Dj^> which has also been 

generally referred to as the 'capillary' model in the 
• ( 91 ) 

literature ' , the latter may be neglected and Eqn. V.12 may then 
be rewritten as 


Daniso “ < D||> = D (i - 4^] [V.13] 

where D = Dy and u = cose 

With the above redefinition, the echo amplitude, Eqn. V. • 
10, may be modified to 


M(2t) = M(0) exp ^ 


2 

3 





[V.14] 


In our model, all values of u = cose are equally 
probable, and therefore we have to perfosrm an integration over the 
angular variable, namely. 


M(2t) = M(0) exp (- 2 T/T 2 ) 


exp y^G^D(l-h^)T^jdjLi 


M(2t) = M(0) exp (- 2 T/T 2 ) exp y^G^DT^f(y)j 


[V.15] 


where Y = [§ r^G^Dr^j ^ and f(y) = ^ 


. 1 


exp (n }dix 


It is thus seen that the tesrm f(y) supplies the 
requisite correction to eqn. V.IO for our case of anisotropic 
diffusion in the randomly oriented long tiibular (or capillas^^) 
stjnicture. 


t 
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V.7.4 Experimental Details 

Using a BRUKER B-KR 322s NMR spectrometer and a 
precalibrated BRUKER B-KR 3002 18 magnetic field gradient unit, 
the rate of decay of spin echo signals has been measured in the 
presence of known magnetic field gradients to estimate the 
self-diffusion coefficient of the 'pool water' of the two reverse 
micellar systems : (i) lecithin/ cyclohexane/H_0 and (ii) 
lecithin/CCl^/H20. Echoes of the water proton signal have been 
observed at time 2 t of a 90° - t - 180°gQOj NMR pulse train, where 
the subscript (90°) indicates that the phase of the RF carrier of 
the second pulse is shifted by 90° relative to the first. has 
been evaluated to an estimated accuracy of ± 5% using the 
Carr-Purcell-Meiboom-Gill pulse train^^®^ corrected for baseline 
drift by the procedure of Hughes 

All signal amplitude measurements have been either read 
out directly from a Tektronix storage oscilloscope or (for weaker 
signals) recorded by an electronic digital read out at the output 
of a high speed signal averager. 

In our diffusion coefficient measurements, it was 
important that during each experiment the NMR probe was 
thermostated to within ±0.5 K. Depending on the choice of field 
gradient strengths, probe temperatures ranged from 296 K to 298 K. 
Larmor frequencies also had to be readjusted slightly from 
experiment to experiment around a central value of 30 MHz. The 
repetition time of the pulse trains was always greater than six 
times the T^ of the sample. 

To minimize errors due to the probe t^aperature and 
other variations to which absolute value of D may be sensitive. 
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3ll results have been converted to the dimensionless ratios D/D 

o 

the diffusion coefficient of the 'pool water' relative to that of 
triply distilled water; here the sample size and diameters of the 
tubes (5 mm OD) have been chosen to be identical. 

V.7.5 Results and Discussion 

We shall first present the value of obtained by CPMG 

method and the T 2 value calculated from the experimental spectrum 

^2 ~ Ti~^ ] • different R values of lecithin/ cyclohexane/H^O 

1/2 

system (Table V.6). 


Table V, 6 Comparison of T 2 values from experimental linewidth 

and from CPMG sequence at different R values for 
the lecithin/cyclohexane/H^O system 


R 

(nearest 

integer 

values) 

T_ from 
linewidth 

in seconds ' 


R 

T 2 from CPMG 

sequence, in 
seconds 

1 



1.0 

0.045 

2 

0.040 


2.0 

0.066 

4 

0.071 


4.3 

0.080 

6 

0.091 


6.0 

0.112 

8 

0.053 


8.6 

0.063 

10 

0.043 


10.8 

0.045 




13.0 

0.041 
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Figs. V.22 and 23 show the linear least square fit 
plots of In jM(2T) + ^ 1 versus through the set of 

experimental data points at different field gradients for 
samples with different R values. We have used values from our 
(more accurate) CPMG experiments, col. 4 of Table V.6, in our 
calculations . 

From the slopes of Figs. V.22 and V.23 the diffusion 
coefficient, D, of water in the reverse micelle is calculated 
using 

slope = j D. 

The standard-errors in the D value measurements range 
between ± 0.02 X 10~^ and ± 0.04 X lo”^^ m^s. 

Samples with R values > 10 became markedly more 

viscous ' and examination under crossed polarisers revealed 

that the system was optically anisotropic. The non-linear trend 

seen in Fig. V. 23(b) clearly corresponds to diffusion in the 

anisotropic phase. To calculate D from this curve, Eqn. V.15 

was fitted to the data points by a two-parameter (D,t) non-linear 

least squares fit where the integral was numerically evaluated by 

a "Simpson's ^rule" sxibroutine^^^^ . Fig. V. 23(b) shows the 

experimental points and the 'best fit' line through them. The 

ratio of this D to the value obtained for triple-distilled 

water under identical experimental conditions turns out to be 

—9 2 -1 

0.60, as shown in Table V.7. The Dq value, 2.35 X 10 ms, 
used throughout in our calculations is an average value of many 
experiments . 




0 200 400 600 800 1000 1200 

X3(xi0'« s3) 



(a) G - 3.6 X 10 T/m , D ^ 21 ± 0 . 01 )xlo“^ i^/s 

(b) G = 3.5 X 10 T/m r D ( n| 9 i±o. 01 )xlO ^ m^/s 




tn IM(2 t) + 
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Table V.7 Variation of self-dif fusion coefficient of water as 

a function of R for lecithin/cyclohexane/H_0 

system. 


R 

D 

( 10 " Vs “^) 

D 

^O 

(D„ = 2.35 X 10 
o 

1.0 

2.23 

0.95 

2.0 

2.21 

0.94 

4.3 

2.16 

0.92 

6.0 

1.91 

0.81 

8.6 

1.80 

0.77 

10.8 

1.61 

0.69 

13.0 

1.42 

0.60 


The variation of relative self-diffusion coefficient 
D/D versus R is pictorially represented in Fig. V.24. From this 
one sees that the water self-diffusion coefficient stays within 
10% of that of free water upto an R value of 6.0, after which it 
dBCiTG^SBci iTSttliBir sharply. This confirms our BarliBr rBsults j 
that in this reverse micellar system around an R value of six 
there is a transition from the isotropic reverse micellar to an 
anisotropic phase which affects the mobility of the water pool. 

For water diffusion in a barrier-free system, the root 
mean square displacement of water molecule in, say, the 
x-direction, <(Ax)^>^/^, occurring during time t is related to the 
self-diffusion coefficient, D, by the well-known Einstein 

relationship 

<(Ax)^>^/^ = 2Dt I 
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If barriers that hinder free motion are present, the measured 
value of D will of course be less than what it would be in the 
unhindered bulk fluid state. In theory, if 2 t can be made short 
enough (roughly less than the root mean square time for a molecule 
to diffuse between barriers or obstructions) , nmr pulse 

techniques should yield a faithful measurement of D for the 
confined fluid. 

Using the above analysis, if diffusion barriers are 

present and sufficiently close together, D becomes a function of 

the applied gradient and the echoes no longer decay as 

For low R value cases (R < 6) , the value of D in each case was 

2 

constant when G was varied by a factor of 20, and all the echo 

3 

decays were linear in x to 2 t - 25 ms for the smallest gradient 
used. Hence, any barrier must be more than [2D(25 ms)] ' apart, 
i.e. 10 microns apart. This is in agreement with the estimates 
of the upper limit of micellar size by our electron microscopy 
experiments published earlier^®^^ . 

The linear least squares fit plots of In jM(2x) + 

versus x^ for lecithin/CCl^/H20 with different R values are shown 

in Figs. V.25 and V.26. The self diffusion coefficients for 

different R values and P/D^ ratio are given in Table V.8 and the 

variation of D/D versus R is shown in Fig. V.27. 

o 

Fig. V.27 shows that initially added water becomes tightly bound 
to the interface, and is then becoming released for more free 
diffusion or exchange with core water thus provides a dynamical 
contrast to the case of the corresponding cyclohexane-containing 
system. 



lntM(2T)4-^l ln[M(2X) + -^ 
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Fig. V. 27 


Variation of relative self-diffusion coefficient 
as a function of added water for lecithin/CCl -/ H-O 
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Table V. 8 


Variation of self— diffusion coefficients of water 
as a function of R for lecithin/CCl^/water system 


R 

D 

(10”\^s“^) 

D 

® —9 2 —1 
(D^ = 2.35 X 10 S ■^) 

2 

1.85 

0.79 

3 

2.00 

0.85 

4 

2.01 

0.86 

5 

2.02 

0.86 

6 

2.01 

0.86 

7 

2.00 

0.85 

8 

1.98 

0.84 

9 

1.97 

0.84 

10 

1.90 

0.81 

12 

1.88 

0.80 


Although the above estimate of the diffusion of trapped 
water is somewhat crude, it should be mentioned here that more 
sophisticated NMR experiments such as the pulsed field gradient 
technique as well as more recent experimental refinements such 
as (i) techniques for diffusion measurement in presence of small 
and large background gradients and (ii) methods for | 
absolute measurement of D without recourse to a standard 
sa»ple<=«) should further improve our understanding of lecithin , 

I 

based reverse micellar systems* j 



Conclusions 


We have demonstrated the successful construction of a 
versatile magnetic field gradient controller capable of 
generating continuous-linear, pulsed-linear and 
oscillating field gradient types. 

Self-diffusion coefficients of water in lecithin based 
reverse micelles could not be measured using our 
magnetic field gradient controller owing to some 
practical difficulties experienced in performing PMFG 
spin echo experiments on the older BRUKER WP-80 FT NMR 
spectrometer. -Nevertheless, the performance results of 
our gradient controller for two types of gradients have 
been illustrated. 

The self -diffusion coefficients of H 2 O in the 
lecithin/cyclohexane/H20 and lecithin/CCl^/H20 systems 
at a number of R values have been measured on a BRUKER 
B-KR 322s FT NMR spectrometer. The trend revealed by 
D/D^ Vs R plots substantiates the earlier findings of 


our laboratory. 
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CHAPTER VI 

OVERALL CONCLUSIONS AND FURTHER PROSPECTS 

We shall now summarize the principal findings of our 

work. 

(1) We have prepared a novel and potentially valuable new reverse 
micellar system, and have characterized it by and NMR 
spectroscopy . 

(2) From our deutron magnetic resonance measurements, we have 
been able to interpret our result in terms of two motionally 
distinct sites for D 2 O molecules (i.e., 'two- site^ model) 
present in DCC-added and DCC-free reverse micelles over a 
range of R-values. We have also observed a thousand fold 
reduction in the exchange rate of water molecules between 
these two sites while going to a lower temperature. 

(3) The manifestation of cross-relaxation terms, and the 
importance of including them in a detailed description of 
motional and structural parameters of aggregated structures, 
have been investigated with special reference to the DOC 
reverse micelles. 

(4) By way of advancing the overall NMR instrumentation 
capabilities of our laboratory, we have accomplished the 
construction of a versatile magnetic field gradient 
controller. Its performance for two types of gradients have 
been illustrated using water samples. 

(5) The measured self-diffusion coefficients of H 2 O in 
egg-lecithin based reverse micelles in different organic 
solvents show that initially added water becomes tightly 
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bound to the interface, and is then becoming gradually 
released for more free diffusion or exchange with core water. 

During the course of our studies, various other 
interesting experimental possibilities suggested themselves. We 
document some of these as further research prospects. 

(1) It would be interesting to study the phase behavior of our 

reverse micelar system lecithin/DCC/Cyclohexane by 

constructing a four component phase diagram at different 
temperatures . 

(2) The indication of structural changes occurring between R 
values of 5 and 6 in this system may be further verified 
visually by recording transmission electron micrographs. For 
example, in dissolving a suitable electron impermeable stain 
in the aqueous core. 

(3) The freezing points of different types of water present in 
the reverse micelles can be measured by performing 
low— temperature differential scanning calorimetry (DSC) 
experiments . 

(4) During the course of our theoretical investigations, we have 
realized the mathematical complexity of handling coupled 
spin systems. Further theoretical efforts at simplifying 
this mathematical problem may perhaps be worthwhile for a 
better understanding of lecithin polar head group dynamics in 


reverse micelles. 



APPENDIX 


c 

C DEUTERIUM LINESHAPE SIMULATION PROGRAM 

C FILE NAME NEW.F 

C INPUTS : OMGAQl , C»!GAQ2 , RK12 , RK2 1 , T2 1 , T22 , TAU , STEP , PINT , VARI 

C OMGAQl =QUADRUPOLAR FREQUENCY-SITE 1 

C 0MGAQ2=QUADRUP0LAR FREQUENCY -SITE 2 

C RJUMP=RATE OF EXCHANGE OF MOLECULES FROM ONE SITE TO ANOTHER 

C T21=T2 FOR THE SITE 1 MOLECULES 

C T22=T2 FOR THE SITE 2 MOLECULES 

C NSTEP=NO. OF OMEGA VALUES CHOOSEN 

C HALFDELONE HALF OF THE DELDA OMEGA 

C NTHETA=NO. OF THETA VALUES CHOOSEN 

C QCC=QUAD. COUPLING CONST. (E**2*Q*Q/H) 

C 

PARAMETER XYD=500 

COMPLEX AA,BB.CC,DETS,J(2),LAMDA(2),R(2.2),RLAMDAC2,2), 

1 S(2,2).S2(2.2),S3(2,2).S4(2,2),S5(2,2),SINVER(2,2),Sl»l.T0TAL 

REAL OMGAIN.OMGAFI, OMGAQl, CWGAQ2, P(2) .PERIND.STEP, TAU, 

1 T21,T22,XAXIS(XYD),YAXIS(XYD) 

INTEGER N,LDA,LDEVEC 
COMMON/AREAl/N , LDA, LDEVEC 
OPEN(UNIT=21,FILE=’XY.O’ ) 

OPEN(UNIT=20,FILE=’DLINS.CXn'’ ) 

OPEN(UNIT=19,FILE=’DLINS. IN’ ) 

READ(19,*) T21.T22,RJUMP 
READ(19,*) TAU 

READ ( 19 , • )NSTEP, HALFDEL, NTHETA, QCC 

C 

DO 2000 VARI=1,NSTEP 
YAXIS( VARI) =0.0 
XAXIS(VARI)=0.0 
2000 CONTINUE 
N=2 
LDA=2 
LDEVEC=2 
T=2.0*TAU 
PI=4.0*ATAN(1.0) 

PIFACT=PI/180.0 

C0UCm=QCC*3. 14*3/4. 0 

DVARI=90. 0/NTHETA 

THETA=0.0 

A=4tt:21+1.0/T21 

B=«K12+1.0/T22 

DO 1000 L00P=1, NTHETA 

THERAD1=THETA*PIFACT 

THERAD2= ( 104 . 0-THETA) •PIFACT 

C0STH1=C0S(THERAD1 ) 

C0HTH2=C0S ( THERAD2 ) 

C0SSQ1=C0STH1*C0STH1 
C0SSQ2=C0STH2*C0STH2 
OMGAQl =COUCON* (3*C0SSQ1-1 ) 

OMGAQ2=COUCON* (3*C0SSQ2-1 ) 

R(1 , 1 )=QffLX(A,OMGAQl ) 

R ( 1 , 2 ) =CMPLX ( -RJUMP , 0 . 0 ) 

R(2, 1)=CMPLX(-RJUMP,0.0) 

R ( 2 , 2 ) =CWPLX ( B , 01«;AQ2 ) 

C 
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c 

c 

c 
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THE ECHO INTENSITY E IS COMPUTED AS: 

« 

E=l. J=1.EXP(-R TAU){EXP{-R TAU)) P 
EXPCR TAU)=S EXP(LAMDA TAU) S( INVERSE). 


CALL EIGEN(R,LAMDA.S.PERIND) 
RLAMDA ( 1 . 1 ) =CEXP ( -LAMDA ( 1 ) *TAU ) 
RLAMDA (2,2) =CEXP ( -LAMDA ( 2 ) ‘TAU ) 
RLAMDA(1,2)=(0. 0,0.0) 
RLAMDA(2,1)=(0. 0,0.0) 

CALL MTMLCCCS, RLAMDA, S2) 

.CALL INVERS(S,SINVER,DETS) 

CALL MTMLCC(S2,SINVER,S3) 

CALL CMPC0N(S3,S4) 

CALL MTMLCC{S3,S4,S5) 

P(l)=l 

P(2)=l 

CALL MTMLCR(S5,P, J) 
0MGAZ0=61.0E6 
OMGAIN=OMGAZO-HALFDEL 
0MGAFI=C»!GA20+HALFDEL 
OMEGA= (OMGAIN-OMGAZO ) 

STEP= (OMGAFI -OMGAIN ) /NSTEP 
DO 200 MM=1, NSTEP 
T0TAL=(0. 0,0.0) 

DO 300 1=1,2 
DO 300 L=l,2 
DO 300 K=1.2 

AA=S ( I . K)*LAMDA (K) *SINVER(K,L) 
BB=OMEGA*OMEGA+LAMDA (K) ♦LAMDA (K) 
CC=(AA/BB)»J(L) 

TOTAL=TOTAL+CC 

CONTINUE 

XAXIS(MM)=C»1EGA 

YAXIS (MM) =REAL (TOTAL ) +YAXIS (MM) 

OMEGA=C»<EGA+STEP 

CONTINUE 

THETA=THETA+DVARI 


CCWTINUE 
WRITE (20, 10 )R 

FOIttlAT(lX,’R-MATRIX’ ,//,lX,4(E11.3,E11.3)) 


WRITE(20,20)LAIffiA 

FC»MAT(1X, 'LAMDA-EIGEN VALUES’ ,//, 1X,2(E11.3,E11.3)) 
WRITE(20,30)S ^ 

FORMAT ( IX, ’S-EIGEN VECTCBS’ ,//,lX,4(E11.3,E11.3)) 
WRITE(20,40)RLAMDA 

FORMAT (IX, 'RLAMDA’ ,//, 1X,4(E11 .3,E11 .3) ) 
WRITE(20,50)S2 

FORMAT ( IX, ’S2’ ,//,lX,4(E11.3,E11.3)) 


WRITE(20,60)SINVER 

F0RMAT(1X,’SINVER’ ,//,lX,4(E11.3,E11.3)) 


WRITE(20,70)DETS 

FORMAT ( IX, ’ DETS’ , // , IX, 2E11 . 3) 

WRITE (20, 80 )S3 

FORMAT (IX, ’S3’ ,//,lX,4(E11.3,E11.3)) 


^ITE(20,90)S4 

F0RMAT(1X,’S4’ ,//,lX,4(E11.3,E11.3)) 
WRITE(20,110)S5 

F0RMAT(1X,’S5’ ,//,lX,4(E11.3.E11.3)) 
HRITE(20,120)P 
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120 FORMATdX/P’ ,//,lX,2F3.1) 

VmiTE(20,130)J 

130 FORMATCIX, ’ J’ ,//, IX. 2(E11. 3,E11 . 3) ) 

WRITE (20, 140)OMGAIN,OMGAFI 

140 FORMAT { IX, ’OMGAIN.OMGAFI’ ,//, IX, 2E15. 4) 

DO 155 I=1,NSTEP 

WRITE(21 , * )XAXIS( I ) , YAXIS( I ) 

WRITE (20 , * )XAXIS ( I ) , YAXIS( I ) 

155 CONTINUE 

C WRITE (21, ’INSTEP 

WRI TE ( 20 , * ) NSTEP , MM 

STOP 

END 

C 

C EIGEN(R,LAMDA,S,PERIND) 

C THIS GIVES EIGEN VALUES IN ’LAMDA’ AND EIGENVECTORS IN ’S’ 

C EIGENVECTORS ARE NORMALISED 

C IN THIS SUBROUTINE ’ EVCCG’ IS AN ’ IMSL’ ROUTINE 

C 

SUBROUTINE EIGEN(R,LAMDA,S,PERIND) 

COMMON/AREAl/N, LDA, LDEVEC 
COMPLEX R(2,2),LAMDA(2),S(2,2) 

INTEGER NOUT 

EXTERNAL EPICG, EVCCG, UMACH,WRCRN 
CALL EVCCG(N,R,LDA,LAMDA,S,LDEVEC) 

PERIND=EPICG(N, N, R, LDA, LAMDA, S, LDEVEC) 

CALL UMACH(2,N0UT) 

CALL WRCRM(’EVAL’ ,1,N, LAMDA, 1,0) 

CALL WRCRNCEVEC’ ,N,N,S,LDEVEC,0) 

WRITE(NOUT,’ (/,R,F6.3)’ ) ’PERFORMANCE INDEX = ’,PERIND 
RETURN 
END 


INVERS(A,AINVER) 

TO FIND INVERSE OF A MATRIX ’A’ 


SUBROUTINE INVERS(A, AINVER,DETS) 
CCOIPLEX A(2,2),AINVER(2.2),DETS 
DETS=A(1,1)*A(2,2)-A(1,2)*A(2,1) 
AINVER ( 1 , 1 )= ( 1/DETS ) ’A (2 , 2 ) 
AINVER(1,2)=-(1/DETS)»A(1.2) 
AINVER(2. 1 )=-(l/rETS)*A(2. 1 ) 
AINVER (2, 2 )= ( 1/DETS ) *A( 1 , 1 ) 
RETURN 
END 


CMPCON(A,B) 

TO FIND CCWPLEX CONJUGATE OF A 


SUBROUTINE CMPCC»I(A,B) 

COMPLEX A(2,2),B(2,2).C0NJG 

DO 100 1=1,2 

DO 100 J=l,2 

B(I,J)=CONJG(A(I,J)) 

CONTINUE 

RETURN 

END 


MTMLCC(A,B,C) 
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C MATRIX MULTIPLICATION COMPLEX. COMPLEX 

C 

SUBROUTINE MTMLCC(A,B,C) 

COMPLEX A(2,2),B(2.2).C(2,2) 

DO 100 1=1,2 
DO 100 J=l,2 
C(I,J)=(0. 0,0.0) 

DO 100 K=l,2 

C(I,J)=C(I,J)+A(I,K)*B(K,J) 

100 CONTINUE 

RETURN 
END 

C 

C MTMLCR(A,B,C) 

C MATRIX MULTIPLICATION CCMdPLEX.REAL 

C 

SUBROUTINE MTMLCR(A,B.C) 

COMPLEX A(2,2),C(2) 

REAL B(2) 

DO 100 1=1,2 
C(I)=(0. 0,0.0) 

DO 100 K=l,2 
C(I)=C(I)+A(I,K)*B(K) 

100 CONTINUE 

RETURN 
END 



